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Software Solutions to Problems in Basic 
Thermodynamics: Part-ll 


Work, Heat and I Law of Thermodynamics 
applied to Closed systems 


4 Work, Heat and I Law of 
Thermodynamics applied to 
Closed systems 


Learning objectives: 

1. Total energy of a system is defined as the sum of internal energy, kinetic energy and 
potential energy i.e. on a unit mass basis : e = u + C 2 /2 + g.z where C is the velocity z is the 
elevation from a datum. 

2. Energy crosses the boundary of a closed system either as Work or Heat, or as both. 

3. Both Work and Heat are path functions 5 , i.e. inexact differentials. 

4. In Thermodynamics, Work is said to be done by a system if the sole effect things external to 
the system can be reduced to the raising of a weight. 

5. ‘Boundary work 5 for a simple compressible system is given by: 


W 12 = 


p dV 


VI 


6. Similarly other types of work, viz. electrical work, shaft work, paddle work, flow work, work 
in stretching a wire, work due to surface tension, magnetization work, free expansion etc. 
have to be considered, if need be. 

7. ‘Heat transfer 5 is energy transfer due to temperature difference only 

8. Conduction, Convection and Radiation are the main modes of heat transfer. Heat transfer 
may occur in one of these modes or, in some cases, one or more modes may be present. 

9. First Law is a statement of conservation of Energy. 

10. First Law for a system undergoing a cycle, and for processes in a closed system are considered. 

11. Different processes for an ideal gas in a closed system (as in a piston- cylinder device) are of 
special interest. 
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Software Solutions to Problems in Basic 
Thermodynamics: Part-ll 

4.1 Formulas used: 

4.1.1 Work: 

Work = Force x distance, N.m (= 1 Joule) 

Work is a path function i.e. an inexact differential. 

4.1 .2 pdV- work or displacement work: 
dW= pdV 

j- V2 

W| 2 = pdV Integration performed on a quasi-static path 

"VI 

4.1 .3 pdV- work in various quasi-static processes: 

(a). Constant pressure (isobaric) process: 

W 12 =p(V2-Y1) 

(b) . Constant volume (isoc boric) process: 

W 12=0 

(c) . For a process in which pV = const.... Isothermal process: 

■' p i 

Wp = pi- VI In £- 

(d) . For a process in which pV = const.... reversible adiabatic or isentropic process: 


Work, Heat and I Law of Thermodynamics 
applied to Closed systems 


7-1 


W 12 = 


pl-Vl - p2 V2 R (T1 - T2) pi VI 


7-1 


7-1 


n- 1 


M3 
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Also, for a perfect gas: 
pv= RT 

And for isentropic process, pv" = const., we have: 
Tv ' = constant 




7-1 


T2 

Tl 



7 


(e). For a process in which pV n = const.. ..polytropic process: 


pi Vl - P 2 V2 pi VI 

■■ ■ | = = 

12 n - 1 n - 1 


n- i 


-er. 


,, w 12 = r -< T1 -^ 

n - 1 


Also: for a polytropic process: 




n-l 
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For a perfect gas: 

pv= RT du=cv-dT 


Then, heat transfer during a polytropic process (for a perfect gas): 
Q = (u2 - ul) + W= c v (T2 - Tl) + R (T1 - T2) 


Simplifying, we get: 

_ y-n T2; 

.. _ i n _t 


Le - Viy=hf' w p°iy 

Polytropic sp. heat: 

Polytr. sp. heat: c n = c v -- — - 

1 — 11 


Mean Effective Pressure (MEP, or p m ): 

_ Are a of Indicator diagram , , 0 

MEP = = — = = — = k ...where k = Spring constant. 

Stroke_lensth 
Indicated Power (IP): 


kW...for a two stroke engine.. where pm is in kPa. A in 
m A 2. L in m. N in RPM.... this is IP for one cylinder 

Note: Put N = N/2 for four stroke engine 
Brake Power (BP): 

2-tc-N-T 

BP= — — — ...where N is RPM. T is Torque 

dO 

Mech. efficiency: 

_ BP 
r -mech ” 


IP = 


Pm ALN 

m 
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applied to Closed systems 


4.1 .4 Other types of Work transfer: 

1. Electrical Power: 


"'dot = EI 


2. Shaft Work: 

"'shaft = Te ...where T is Torque, co is angular velocity 

3. Paddle work or Stirring work: 


W = 


m- s dZ = 


1 


TdS 


i 


4. Flow Work: 

W flow = p-v ...per unit mass 

5. Work done in stretching a wire: 


W= - JdL ..where J is the tension. dL is expansion of wire 

1 

6. Work done in changing the area of a surface film: 


W= - 


g dA ..where a is the surface tension N/m) 


1 


7. Work done in magnetization of a paramagnetic solid: 


W= - EdI ..where H is the field strength and I is the component of magnetization 
^1 field in the direction of the field 


3. Work done in Free expansion: 


"fr ee _ e ,pn=<> ...since there is no resistance to the fluid at boundary 
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4.1.5 Heat Transfer, Q: 

Q is positive while flowing into the system; 

W is positive if work is done by the system. 

Heat Transfer Q 12 : 

Heat transfer is a path function. 


Ql2- | 

j l 


T ds where T is in K : s is entropy 


Specific heat, c: 

It is the amount of heat required to raise a unit mass of substance through a unit rise in temperature. 

C= — J/kg.K 

m-At 
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For a gas, for a constant pressure, reversible non-flow process: 
dQ = m- Cp- dl 

For a gas, for a constant volume, reversible non-flow process: 
dQ = m- e v - dT 


4.1 .6 First Law for a system undergoing a cycle: 

ZW= I IQ ....for a cycle. J = 1 in S.l. Units. i.e. 1 N.m = 1 Joule. 

4.1 .7 First Law for a closed system undergoing a change of state: 

Q - W= AE 

or: 

Q = AE + W 

4.1 .8 First Law is a statement of conservation of Energy. 

Energy is a property of the system; it is therefore, a point function. 

Considering only the kinetic, potential and internal energies, Total energy is: 

E= E k “ E p + U int 


4.1.9 For an Ideal gas: 

Internal Energy U is a function of T only. 

We write: 

dQ = dE + dW 
i.e. dQ=dU-dW 

i.e. dQ = dU - p dV ...when only pdV work is present 
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Enthalpy, h: 

h = u - p v J/kg 
For a perfect gas: 

h = c v T + R T = (c v + r)-T = c -T 

4.1 .1 0 Fist Law for non-flow processes or for Closed systems: 
For reversible, const, volume process: 


Q = (u2 - ul) - W ... where W = p dv 


....But. W - 0, since dV - 0 


Therefore: Q = u2 - ul = <y(T2 - Tl) . . . J/kg 

For reversible, const, pressure process: 


Q = (u2 - ul) + W .. .where W = p dv 


....But. W = p.(v2-v1) 


Therefore: Q = h2 - hi = c p (T2 - Tl) ...J/kg 


For reversible, Isothermal process: 


Q = (u2 - ul) — W ...where W= p dv 


Therefore: Q = c v (T2 - Tl) + W = 0 + W 


i.e. 



For reversible, adiabatic process: 


Q = (u2 - ul) + W ...where W = p dv 


But. Q = 0 for adiabatic process. 


Therefore: 0=(u2-ul) + W 


i.e. w = ul - u2 ....for any adiabatic process 
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And. for rev. adiab. process: p V = const 


And: 


W = 


Pr v l ~~ P2" v 2 R (T1 - T2) 


7 - 1 


7 - l 


4.2 Now, let us work out a few problems with EES: 

“Prob.4.1. A perfect gas is undergoing a process in which T is proportional to V A (-2/5). Calculate the 
work done by the gas in going from state 1 in which pressure is 100 bar and volume is 4 m A 3 to the state 
2 in which volume is 2 m A 3. Also calculate the final pressure. [VTU-BTD-Dec-06- Jan-07]” 

EES Solution: 

“Data:” 

Pl = 100E05[Pa] 

Vl=4[m A 3] 

V2=2[m A 3] 
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“Calculations:” 

“We have: 

T = kl * V A (-2/5); But, PV = RT for perfect gas. 

i.e. P.V / R = kl * V A (-2/5) 

i.e. P.V A (7/5) = k where k = kl * R, a const.” 

PI * Vl A (7/5) = P2 * V2 A (7/5) 
k = PI * Vl A (7/5) 

Work=integral(kW A (-7/5),V,Vl,V2) “...using the built-in integral function of EES” 

“Note: In the above, we calculate Work as Integral of P.dV. So, P is expressed as a function of V. VI 
and V2 are limits of integration, i.e. from VI to V2.” 

Now, hit F2 to calculate. 

Results: 

Unit Settings: SI K Pa J mass deg 

k = 6.9G4E+07 PI = 1 .000E+07 [Pa] 

V - 2 [m 3 ] VI = A [m 3 ] 

Work = -3.1 95E+07 [J] 


P2 = 2.639E+07 [Pa] 
V2 =2 [m 3 ] 


Thus: 

Work = -3.195E07 W ... negative sign indicates that work is done on the system.... Ans. 
Final pressure, P2 = 2.639E07 Pa = 263.9 bar ... Ans. 
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Additionally, plot the variation of Work and P2 as the final volume varies from 0.5 m A 3 to 3 m A 3: 
First, compute the Parametric Table: 


(Note that we have written the absolute value of Work done, i.e. without the negative sign). 



Now, plot the results: 


x 1G 7 



V2 [m 3 ] 
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x 10 7 



V2 [m 3 ] 
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“Prob.4.2. An engine cylinder has a piston of area 0.12 m A 2 and contains gas at a pressure of 1.5 MPa. 
The gas expands according to a process which is represented by a straight line on a p-V diagram. The 
final pressure is 0.15 MPa. Calculate the work done by the gas if the piston stroke is 0.3 m. [VTU-BTD- 
July/Aug.2004-New-Scheme] ” 



Fig. Prob.4.2 


EES Solution: 

“Data:” 

PI = 1.5E03 [kPa] 

P2 = 0.15E03 [kPa] 

A = 0.12 [m A 2] “...piston area” 
L = 0.3 [m] “...stroke” 


“Calculations:” 

DELTAV = A * L “[m A 3]” 

Work = P2 * DELTAV + (PI - P2) * DELTAV/2 “[kj]” 


Solution: 

Unit Settings: SI C kPa kJ mass deg 

A =0.12 [m 2 ] 

P2 =150 [kPa] 


AV = 0.036 [in 2 ] 
Work = 29.7 [kJ] 


L = 0.3 [m] 


PI = 1500 [kPa] 


Thus: 

Work done by the gas = 29.7 kJ....Ans. 
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(b) Plot the variation of Work as the final pressure P2 varies from 100 kPa to 1000 kPa: 
First, compute the Parametric Table: 



Now, plot the results: 



P2 [kPa] 
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a Prob.4.3. A spherical balloon of 1 m dia contains a gas at 200 kPa pressure. The gas inside the balloon 
is heated until the pressure reaches 500 kPa. During the process of heating, the pressure of gas inside 
the balloon is proportional to the cube of the diameter of the balloon. Determine the work done by the 
gas inside the balloon. [VTU-BTD- June- July-08]” 

EES Solution: 

“Data:” 

Pl=200[kPa] 

Dl = l[m A 3] 

k=200 “...since PI = k * D1 A 3” 

P2=500[kPa] 


“Calculations:” 

P2=k*D2 A 3 “...finds D2” 

Vl=(pi/6)*Dl A 3”[m A 3]” 

V2=(pi/6)*D2 A 3”[m A 3]” 

Work = integral(k*6 * V/pi,V,Vl,V2) “[kJ]....Note the use of built-in integral function of EES” 


“Note: In the above, we calculate Work as Integral of P.dV. So, P is expressed as a function of V. VI 


and V2 are limits of integration, i.e. from VI to V2” 
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Results: 

Unit Settings: SI C kPa kJ mass deg 

D1 =1 [m 3 ] D2 = 1.357 [m] k = 200 PI =200 [kPa] 

P2 = 500 [kPa] V = 1.309 VI = 0.5236 [m 3 ] V2 = 1.309 [m 3 ] 

Work = 274.9 [kJ] 


Thus: 

Work done by the gas = 274.9 kJ .... Ans. 


(b) Plot the variation of Work as the final pressure P2 varies from 300 kPa to 750 kPa: 
First, compute the Parametric Table: 
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Now, plot the results: 



“Prob.4.4. A spherical balloon of 1 m dia contains a gas at 1.5 bar pressure. Due to heating, the pressure 
reaches 4.5 bar. During the process of heating, the pressure is proportional to the cube of the diameter 
of the balloon. Determine the work done by the gas inside the balloon. [VTU-BTD-Feb. 2002]” 

EES Solution: 

This is similar to the previous problem. 

“Data:” 

Pl=150[kPa] 

Dl = l[m A 3] 

k=150 “...since PI = k * D1 A 3” 

P2=450[kPa] 

“Calculations:” 

P2=k*D2 A 3 “...finds D2” 

Vl=(pi/6) >f Dl A 3“[m A 3]” 

V2=(pi/6)*D2 A 3“[m A 3]” 

Work = integral(k * 6 * V/pi,V,Vl,V2) “[kJ]....Note the use of built-in integral function of EES” 
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“Note: In the above, we are calculating Work as Integral of P.dV. So, P is expressed as a function 
of V. 

VI and V2 are limits of integration, i.e. from VI to V2.” 

Results: 

Unit Settings: SI C kPa kJ mass deg 

D1 =1 [m 3 ] D2 =1.442 [m] k = 1 50 

P2 = 450 [kPa] V =1.571 VI = 0.5236 [m 3 ] 

Work =314.2 [kJ] 


PI =150 [kPa] 
V2 =1.571 [m 3 ] 


Thus: Work done by the gas = 314.2 kJ .... Ans. 


“Prob.4.5. A spherical balloon of dia 0.5 m is initially at a pressure of 100 kPa. Due to heating, pressure 
increases to 400 kPa during which the inside pressure varies directly proportional to the square of the 
diameter of the balloon. Determine the displacement work during this process. [VTU-BTD-July-2007]” 
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Work, Heat and I Law of Thermodynamics 
applied to Closed systems 

EES Solution: 

“Data:” 

Pl=100[kPa] 

Dl=0.5[m A 3] 

Pl=k*Dl A 2 “...determines k” 

P2=400[kPa] 

P2=k*D2 A 2 “...determines D2” 

“Calculations:” 


Software Solutions to Problems in Basic 
Thermodynamics: Part-ll 


Vl=(pi/6)*Dl A 3“m3” 

V2=(pi/6)*D2 A 3“m3” 

Work = integral(k >f (6 >f V/pi) A (2/3),V,Vl,V2) “kj using the built-in function integral of EES” 

“Note: In the above, we are calculating Work as Integral of P.dV. So, P is expressed as a function 
of V. 

VI and V2 are limits of integration, i.e. from VI to V2” 

Results: 

Unit Settings: SI C kPa kJ mass deg 

D1 = 0.5 [m 3 ] D2 = 1 [m] k = 400 PI =100 [kPa] 

P2 = 400 [kPa] V = 0.5236 [m 3 ] VI = 0.06545 [m 3 ] V2 = 0.5236 [m 3 ] 

Work =121.7 [kJ] 


Thus: Work done by the gas = 121.7 kj .... Ans. 


“Prob.4.6. A quasi-static process occurs such that P = (V A 2 + 8/V), where P is the pressure in bar and V 
is the volume in m A 3. Find the work done when volume changes from 1 m A 3 to 3 m A 3. [VTU-Jan.2004]”’ 

EES Solution: 

“Data:” 

p=(v A 2+8/v) “bar” 
vl = l”m A 3” 
v2=4”m A 3” 
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“Calculation:” 

W=10 A 5*(integral(p,v,vl,v2)) “J....uses the built-in integral function of EES” 

Results: 

Unit Settings: SI K kPa kJ molar deg 

p = 1 8 [bar] v = A [m 3 ] vl = 1 [m 3 ] v2 = A [m 3 ] 

W = 3.209E+06 [J] 


Thus: Work done by the gas = 3.209E06 J .... Ans. 


(b) In addition, plot Work against V2, as V2 changes from 1.5 m A 3 to 5 m A 3: 
First, compute the Parametric Table: 
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Now, plot the graph: 


x 10 3 




A cate -Lucent 


www.alcatel-lucent.com/careers 


What if 
you could 
build your 
future and 
create the 
future? 




One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it’s strange that 
devices ever had to be “plugged in." To obtain that status, there 

needs to be “The Shift". 



Click on the ad to read more 


Download free eBooks at bookboon.com 



Software Solutions to Problems in Basic 
Thermodynamics: Part-ll 


Work, Heat and I Law of Thermodynamics 
applied to Closed systems 


a Prob.4.7. 1 kg of air at 15 C and 100 kN/m A 2 is compressed isentropically to 600 kN/m A 2. Determine 
the final temp and the work done. If the air is cooled to 15 C at constant pressure, calculate the heat 
transferred. Take gamma = 1.4, cp = 1.0213 kJ/kg.K, R = 0.287 kJ/kg.K. [VTU-Jan. 2005]” 

EES Solution: 

“Data:” 

m = 1 “[kg]” 

Tl=15+273 “[K]” 
pl=10 A 5 “[Pa]” 
p2=6*10 A 5 “[Pa]” 

gamma=1.4 “...ratio of sp. heats, (cp/cv) for air” 
cp=1.0213*10 A 3“ [ J/kg.K] . . . .sp. heat” 

R=287 “[J/kg.K].... gas const.” 

“Calculations:” 

T2/T1 = (p2/pl) A ((gamma-l)/gamma)“...temp ratio for an isentropic process.... determines T2” 
W_ad=R*(Tl-T2)/(gamma-l) “Adiabatic work” 

Q=m*cp*(T2-Tl) “[J]....heat transferred, when cooled to T1 from T2, at const, pressure” 


Results: 


Unit Settings: SI K kPa kJ molar deg 

cp = 1 021 [J.kg-K] y = 1 .4 


p2 = 600000 [Pa] 

Q = 196632 [J] 


T2 = 480.5 [K] 


W ad =-138141 [J] 


m = 1 [kg] 

R = 237 [J/k gK ] 


pi = 100000 [Pa] 
T1 = 288 [K] 


Thus: 

Final temp T2 = 480.5 K .... Ans. 

Work done = -138141 J .... Ans. negative sign indicating work done on the system 
Heat transferred in const, pressure cooling, Q = 196632 J ... Ans. 
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(b) In addition, as p2 varies from 100 kPa to 1000 kPa, plot the variation of T2, W and Q 
against p2: 

First, compute the Parametric Table: 


Table 1 


► 

1 .10 

i IZ 

p2 

[Pa] 

2 Z 

T2 

[K] 

3 Z 

w ad 

[J] 

E o 

R 

Run 1 

100000 

283 

0 

0 

Run 2 

200000 

351 1 

45257 

64419 

Run 3 

300000 

394.2 

76196 

108459 

Run 4 

400000 

428 

100426 

142948 

Run 5 

500000 

456 1 

120640 

171721 

Run 6 

500000 

430 5 

133141 

196632 

Run 7 

700000 

502 2 

153666 

218730 

Run 8 

800000 

521.7 

167677 

238675 

Run 9 

900000 

539.6 

180488 

256910 

Run 10 

1000000 

556 

192319 

273750 


Note that in the above Table, we have taken the absolute value of Work, with the understanding that 
it is the work done on the gas during compression. 


Now, plot the results: 



P2 [Pa] 
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xIO 3 



p2 [Pa] xjo! 
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xlG 3 



p2 [Pa] xjo^ 


“Prob.4.8. 5 kg of Nitrogen at 100 C is heated in a reversible, non-flow, constant volume process till 
the pressure becomes three times the initial pressure. Determine: (i) final temp (ii) change in internal 
energy (iii) change in enthalpy, and (iv) heat transfer. Take R = 0.297 kJ/kg.K, cv = 0.7435 kJ/kg.K. 
[VTU-Jan. 2004]” 

EES Solution: 

“Data:” 

m = 5 “kg” 

PressureRatio = 3 “pr. ratio= p2/pl” 

Tl = 100+273 “K” 
cv=743.5 “J/kg.K” 

R=297 “J/kg.K” 
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“Calculations:” 

T2 = PressureRatio * Tl“.. finds T2... since pl/Tl = p2/T2 at const, volume” 
cp - cv = R “..for Ideal gas... finds cp“ 

DELTAU=m * cv * (T2-T1) “J... change in internal energy” 

DELTAH=m * cp * (T2-T1) “J... change in enthalpy” 

W = 0 “..since it is a const, volume process” 

Q=DELTAU +W “J.. from I law for a closed system” 


Results: 


Unit Settings: SI K kPa kJ molar deg 

cp = 1 CM1 [J/kg-K] cv = 743.5 [J/kg-K] 

AU = 2.773E+06 [J]| m = 5 [kg] 


Q = 2.773E+0G [J] 
T2 =1119 [K] 


R = 297 [J/kg-K] 
W = 0 [J] 


AH = 3.881 E+06 [J] 
PressureRatio = 3 
T1 = 373 [K] 


Thus: 

Final temp, T2 = 1119 K, Change in Int. energy, DELTAU = 2.773E06 J, 

Change in enthalpy, DELTAH = 3.881E06 J, Heat transfer, Q = 2.773E06 J .... Ans. 


“Prob.4.9. 1 kg of air contained in a closed system at 100 kPa and 300 K is compressed isothermally till 
the volume halves. During the process, it is also stirred with a Torque of 1 N.m at 400 RPM for 1 hour. 
Calculate the net work done on the system. Assume R = 0.285 kJ/kg.K. [VTU-July 2003]” 

EES Solution: 

“Data:” 

m=l “kg” 
pl=100*10 A 3 “Pa” 

T 1=300 “K” 

pi =0.5*p2 “...since pi. VI = p2.V2 at constant T” 

N=400*60 “Revolutions in one hour” 

T= 1 “N.m .... torque” 

R=285 “J/kg.K” 
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“Calculations:” 

W_iso = R * T1 * ln(pl/p2) “J.... isothermal work on the system” 
Wl=-2 * pi * N * T CC J.... stirring work on the system” 
W_net=W_iso+Wl “J.... net work on the system” 


Results: 

Unit Settings: SI 

m = 1 [kg] 

R = 235 [J/kg-K] 
W iso =-59264 [J] 


kPa kJ molar deg 

N =24000 [rev.] 
T =1 [N.m] 

W ne t = -21 0061 [J] 


pi = 1 00000 [Pa] 
T1 = 300 [K] 


p2 = 200000 [Pa] 
W1 =-1 50796 [J] 


Thus: 

Net work done on the system = -210061 J Ans. Negative sign indicating that work is done on the 

system. 


“Prob.4.10. 1.5 kg of a gas undergoes a quasi-static process, in which the pressure and sp. vol. are related 
by the equation: p = a - b.v, where a and b are constants. The initial and final pressures are 1000 kPa and 
200 kPa respectively. The corresponding volumes are 0.2 m A 3 and 1.2 m A 3. The specific internal energy 
of the gas is given by the relation: u = 1.5 p v - 35, where u is in kj/kg, p is in kPa, and v is in m A 3/kg. 
Find the magnitude and direction of heat transfer and the max. internal energy of the gas during the 
process. [VTU-Jan. 2005]” 

EES Solution: 

“Data:” 

m=1.5 “kg” 

“u=1.5 * p * v - 35 internal energy” 

pl = 1000 “kPa ... initial pressure” 
vl=0.2/m “m3/kg . initial sp. volume” 
p2=200 “kPa . . . final pressure” 
v2= 1.2/m “m3/kg ... final sp. volume” 
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“Calculations:” 

“To find a and b:” 

PI = a - b * initial pressure” 

P2 = a - b * v2 “...final pressure” 

“To find W, Q and DELTAU:” 

DELTAU = U2-U1 CC J change in internal energy” 

p = a- b*v“ rein, between p and v, by data” 

W = m * 10 A 3 * integral(p,v,vl,v2)“J using the built-in function ‘integral’ of EES” 

Q = W + (U2 - Ul) “J ... by I Law for a closed system” 

U1 = m * (1.5 * pi * vl - 35) * 10 A 3 “J internal energy at state 1” 

U2 = m * (1.5 * p2 * v2 - 35) * 10 A 3 “J ... internal energy at state 2” 



L / i 
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Results: 

Unit Settings: SI K kPa kJ molar deg 


a= 1 1 GO 
p = 200 [kPa] 

U1 = £47500 [J] 
v£ = 0.8 [m 3 ] 


b = 1 200 

pi = 1000 [kPa] 
U2 = 307500 [J] 
W = 600000 [J] 


m = 1 .5 [kg] 

Q = 660000 [J] 
v = 0.8 [m 3 ] vl = 0.1333 [m 3 ] 


AU = 60000 [J/kg] 
p2 = 300 [kPa] 


Thus: 

Q = 660000 J .... Ans. It is positive, indicating that heat is transferred to the system. 
U2 = max. int. energy = 307500 J .... Ans. 


(b) Plot Q, W and DELTAU as final pressure p2 varies from 500 to 100 kPa: 


First, compute the Parametric Table: 


l> 

E 

2 E 

3 E 

4 E 

P* 

W 

AU 

Q 

1..9 

[kPa] 

[J] 

[J/kg] 

[J] 

Run 1 

500 

4S8750 

318750 

787500 

Run 2 

450 

438438 

233053 

737500 

Run 3 

400 

525000 

270000 

735000 

Run 4 

350 

543438 

231553 

730000 

Run 5 

300 

558750 

183750 

752500 

Run 6 

250 

585338 

125563 

712500 

Run 7 

200 

500000 

60000 

650000 

Run 3 

150 

510338 

-15333 

535000 

Run 3 

100 

513750 

-101250 

517500 
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Next, plot the results: 


x 10 3 



p2 [kPa] 


xIO 3 



p2 [kPa] 


Note that after p2 = approx. 160 kPa, DELTAU becomes negative. 
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a 


x 10 3 



p2 [kPa] 


Note that after p2 = approx. 425 kPa, Q decreases. 
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“Prob.4.11. A piston-cylinder mechanism containing a fluid has a stirring device as shown. The piston 
is frictionless and held against the fluid by atm. pressure of 101.325 kPa. The stirring device is turned 
10000 revolutions with an average torque against the fluid, of 1.275 N.m. The piston is 0.6 m dia and it 
moves by 0.8 m. Calculate the net work transfer. [VTU-July 2002]” 

0.8 m 

* ► 

P= 101.325 kPa 



EES Solution: 

“Data:” 

N= 1 0000 “revolutions” 

T=1.275[J] 
d = 0.6 [m] 

L=0.8[m] 

p=101.325E03[Pa] 

“Calculations:” 

Wl=-2 * pi * N * T CC J . . . stirring work done on the system” 
W2=F * L “J. . .boundary work done by the system” 

F=p * A CC N... force exerted on the piston by atm.” 

A=(pi/4) * (d) A 2 “m A 2 .... area of piston” 

W_tot = W1+W2 cc J....net work” 


L = 0.8 [m] 

W1 =-80111 [J] 


Results: 

Unit Settings: SI K kPa kJ molar deg 

A =0.2827 [m 2 ] d = 0.G[m] 

N = 10000 [Rev.] p = 101325 [Pa] 

W2 = 2291 0 [J] |w tot = -571 91 [J] 


F = 28649 [N] 
T = 1.275 [J] 


Thus: Net work done = -57191 J, negative sign indicating work done on the system. 
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“Prob.4.12. A closed system undergoes a cycle composed of 4 processes 1-2, 2-3, 3-4 and 4-1. The 
energy transfers are as tabulated: 


Process 

Q(kJ/min.) 

W (kJ/min.) 

AU (kJ/min.) 

1-2 

400 

150 

- 

2-3 

200 

- 

300 

3-4 

-200 

- 

- 

4-1 

0 

75 

- 


(i) complete the Table (ii) determine the rate of work in kW [VTU-Jan. 2004]“ 

EES Solution: 

“Data:” 

“Process 1-2:” 

Q_12=W_12+DELTAU_12 “kj/min” 

Q_12=400 “kj/min” 

W_12=150 “kj/ min ” 

“Process 2-3:” 

Q_23=W_23+DELTAU_23 “kj/min” 

Q_23=200 “kj/min” 

DELTAU_23=300 “kj/min” 

“Process 3-4:” 

Q_34=W_34+DELTAU_34 “kj/min” 

Q_34=-200 “kj/min” 

“Process 4-1:“ 

Q_4 1 = W_4 1 +DELTAU_4 1 “kj/min” 

Q_41=0 “kj/min” 

W_41=75 “kj/min” 

Q_12+Q_23+Q_34+Q_41=W_12-i-W_23-i-W_34-i-W_41 “...First Law for the whole cycle” 

“Net Heat and Work in cycle:” 

Q_net = (W_12 + W_23 + W_34 + W_41) / 60 “[kj/s]” 

W_net = (W_12+W_23+W_34+W_41) / 60 “[kj/s]” 
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Results: 

Unit Settings: SI K kPa kJ molar deg 


AU 12 =250 [kJ/min] 
AU 41 =-75 [kJ/min] 
Q 34 = -200 [kJ/min] 
W 12 =150 [kJ/min] 
W 41 =75 [kJ/min] 


AU 2 3 = 300 [kJ/min] 
Q-| 2 =400 [kJ/min] 
Q 4 1 = 0 [kJ/min] 

W 23 =-100 [kJ/min] 


W net = 6.667 [kW] 


AU 34 = -475 [kJ/min] 
Q 2 3 = 200 [kJ/min] 
Q net = 6.667 [kW] 
W 34 = 275 [kJ/min] 


Thus: 

Following is the completed Table: 


Process 

Q(kJ/min.) 

W (kJ/min.) 

AU (kJ/min.) 

1-2 

400 

150 

250 

2-3 

200 

-100 

300 

3-4 

-200 

275 

-475 

4-1 

0 

75 

-75 


And, W_net = 6.667 kW .... Ans. 
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“Prob.4.13. During a reversible, constant pressure process in a closed system with p = 105 kPa, properties 
of the system change from VI = 0.25 m A 3, tl = 10 C to V2 = 0.45 m A 3, t2 = 240 C. Specific heat at 
const, pressure, cp is given by: cp = (0.4 + 18 / (t + 40)) kJ/kg.C. Assuming the mass of the system as 
1 kg, determine: (i) heat transfer (ii) work transfer (iii) change in internal energy, and (iv) change in 
enthalpy. [VTU-Jan. 2003]” 

EES Solution: 

“Data:” 

p=105 * 10 A 3 “Pa” 

Vl=0.25 “m3” 

V2=0.45 “m3” 
tl= 10 “C” 
t2=240 “C” 

Cp=(0.4+18/(t+40)) * 10 A 3 “J/kg.C” 

“Calculations:” 

Q=integral(Cp,t,tl,t2) “J.... finds heat transfer” 

W = p * (V2-V1) “J finds work transfer” 

Q = W + DELTAU “. . ..by I Law for a closed system” 

DELTAH = Q “J... change in enthalpy for const, pressure process” 


Results: 

Unit Settings: SI K kPa kJ molar deg 


Cp = ABA. 3 [J/kg-C] 

AH = 1 2301 1 [J/kg] 

AU = 1 0201 1 [J/kg] 

p = 1 05000 [Pa] 

Q =123011 [J/kg] 

t = 240 [C] 

tl =10 [C] 

t2 =2-10 [C] 

VI = 0.25 [m 3 ] 

V2 = 0.45 [m 3 ] 

W = 21 000 [J/kg] 



Thus: 

Q = 123011 J/kg, W = 21000 J/kg, AU = 102011 J/kg, AH = 123011 J/kg Ans. 
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(b) As t2 varies from 200 C to 400 C, plot the variation of Q: 
First, compute the Parametric Table: 


Table 1 

► 

1..11 

t2 

[C] 

2 E 

Q 

[J/kg] 

Run 1 

200 

104236 

Run 2 

220 

113677 

Run 3 

240 

123011 

Run 4 

260 

132253 

Run 5 

280 

141414 

Run 6 

300 

150506 

Run 7 

320 

159535 

Run 8 

340 

168508 

Run 9 

360 

177431 

Run 10 

380 

186309 

Run 11 

400 

195147 | 


Now, plot the results: 



O 


xIO 3 



12 [C] 
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a Prob.4.14. A system receives 200 kj of heat at constant volume. Then, it rejects 70 kj of heat at constant 
pressure and work done on the system being 50 kj. If the system is restored to the initial state by an 
adiabatic process, how much work will be done during the adiabatic process? Calculate the change in 
internal energy for the above mentioned processes and draw the p-V diagram. [VTU-Feb. 2002]” 

EES Solution: 

“Let: 

process 1-2 : constant vol. process, 
process 2-3: constant pr. process, 
process 3-1: adiabatic process.” 

P 



“Process 1-2:” 

Q_12 = 200“kJ.... heat transfer, by data” 

W_12 = 0 “kj ... for const, vol. process” 

Q_12 = W_12 + DELTAU_12 “...by First Law for the process” 

“Process 2-3:” 

Q_23 = -70“kJ.... heat rejected, by data” 

W_23 = -50 “kj... work done on the system, by data” 

Q_23 = W_23 + DELTAU_23 “...by First Law for the process” 

“Process 3-1:” 

Q_31 = 0 “kj since adiabatic, by data” 

Q_31 = W_31 + DELTAU_31 “...by First Law for the process” 

DELTAU_12 + DELTAU_23 + DELTAU_31 = 0 “...since it is a closed cylce” 
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Results: 

Unit Settings: SI K kPa kJ molar deg 


AU 12 = 200 [kJ] 

AU 2 3 =-20 [kJ] 

AU 31 =-180 [kJ] 

Q 12 =200 [kJ] 

Q 23 = -70 [kJ] 
W 31 =180 [kJ] 

□31 = 0 [kJ] 

W-| 2 =0 [kJ] 

W 23 =-50 [kJ] 


Thus: 

Work done in adiabatic process, W_31 = 180 kJ ... Ans. 

Change in internal energies: AU 12 = 200 kJ, AU 23 = -20 kJ, AU 31 = -180 kJ .... Ans. 


“Prob. 4.15. Consider the system shown in fig. Initial conditions of the gas are: VI = 0.1 m A 3, pi = 200 
kPa. Ambient pressure: 100 kPa and the spring exerts a force which is proportional to the displacement 
from its equilibrium position. The gas is heated until the volume is doubled, at which point p2 = 600 
kPa. Determine the work done by the gas. [VTU-Aug. 2001]” 
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p = 200 kPa 
V = 0.1 m A 3 


kx 


-WMr- 


p_atm = 
100 kPa 


> x 


EES Solution: 

“Data:” 

pi =200 “kPa” 

VI = 0.1 “m3” 

V2 = 2* VI “m3” 
p2=600 “kPa” 

“Calculations:” 

W_tot=((pl+p2)/2) * (V2-V1) “kj” 

Results: 

Unit Settings: SI K kPa kJ molar deg 

pi = 200 [kPa] p2 = G00 [kPa] VI = 0.1 [m 3 ] V2 = 0.2 [m 3 ] 
Wtot = 40 [kJ] 


Thus: 

Work done by the gas = 40 kj .... Ans. 


w Prob.4.16. A fluid is heated reversibly at a constant pressure of 1.03 bar until it has a specific volume of 
0.1m A 3/kg. It is then compressed reversibly according to the law pv = constant to a pressure of 4.2 bar, 
then allowed to expand reversibly according to the law: pv A 1.2 = constant to the initial conditions. The 
work done in the constant pressure process is 820 J and the mass of the fluid present is 0.2kg. Calculate 
the net work done on or by the fluid in the process and sketch cycle on a p-v diagram.” 
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EES Solution: 

Data: 

pi = 1.03E05 “Pa” 

p2 = pl“.... const, pr” 

v2 = 0.1“m A 3.... sp. volume” 

w_12 = 820“J” 

p3 = 4.2E05 “Pa” 

mass = 0.2“kg” 

“Calculations:” 

p2 * (v2-vl) = w_12 “gives vl” 
p3 * v3 = p2 * v2 “gives v3“ 

w_23 = (p2 * v2) * ln(v3/v2)“J...work done in isothermal process 2-3” 
w_31 = (p3 * v3 - pi* vl)/(1.2 - 1)“J... work done in polytropic process 3-1” 
w_net = w_12 + w_23 + w_31“J...net work done in the cycle, per kg” 
Work_net = w_net * mass“J...net work done in the cycle for mass = 0.2 kg” 


Results: 

Unit Settings: SI C kPa kJ mass deg 


mass = 0.2 [kg] 
p3 = 420000 [Pa] 
v3 = 0.02452 [m 3 ] 
w 2 3 = -1 4477 [J] 


pi = 103000 [Pa] 
vl = 0.09204 [m 3 ] 
Work ne t =-1911 [J] 
w 31 = 41 00 [J] 


p2 = 103000 [pA] 
v2 = 0.1 [m 3 ] 
w-|2 = 320 [J] 
w n et = "3357 [J] 



Thus: 

Net work done = -1911 J .. negative sign indicating work done on the system.... Ans. 
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a Prob.4.17. A fluid system undergoes a non flow frictionless process from VI = 6 m A 3 and V2 = 2 m A 3. 
The pressure and volume relation during the process is given by following relation, P in N/m2 where V is 
in m A 3. Determine the magnitude and direction of work transfer during the process. [VTU-Sept. 2009]” 

EES Solution: 

“Data:” 

P= (15/V) + 2 “...relation between P and V” 

VI = 6 “m A 3” 

V2 = 2 “m A 3” 


“Calculation:” 

W = integral (P, V, VI, V2) “J....work done ...using the built-in integral function of EES” 

Results: 

Unit Settings: SI C kPa kJ mass deg 

P = 9.5 [Pa] V = 2 [m 3 ] VI = 6 [m 3 ] V2 = 2 [m 3 ] 


W = -24.48 [J] 
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Thus: 

Work done = -24.48 J.... negative sign indicating work done on the system Ans. 


4.3 Now, let us solve a few problems with TEST: 

Prob. 4.18. Nitrogen at an initial state of 300 K, 150 kPa, and 0.2 m A 3 is compressed slowly in an 
isothermal process to a final pressure of 800 kPa. Determine the work done during the process. [Ref.l] 

TEST Solution: 

Let us solve this problem with The Expert System for Thermodynamics (TEST): 

Following are the steps: 

1. Start TEST after logging in to www.thermofluids.net . We get the following Greeting screen: 



thermofluids. n 
Mirror 1 
Mirror 2 
Mirror 3 
Mirror 4 


Quad 
W Btu 


Unit Converter 


Tables B. Charts 


Login Message: You are now logged in using your personal license (valid until Thu Dec 17 
2015). If you encounter any problem in accessing a resource, please write to us at 
support@thermof!uids. net. 

As an educator, you can now acquire a group license (from the MyAccount page) and distribute 
TEST accounts to your students who can use TEST from anywhere, not just from within the 
campus. Please take a look at the key problems in the Problems module. If you assign 
supplemental homework from the Problems module (students, of course, cannot see the 
answers), TEST can monitor progress, do the grading, and generate a detailed report. 


Greetings, Dr. Muliya, 


registered user since 2005 


" Thermo dynamics is a fuimj’’ subject. The first time you go through it, you 
don't understand it at all. The second time you go through it, you think you 
understand it, except for one or two small points. The third time you go 
through it, you know you don't understand it, but by that time you are so used 
to it, it doesn't bother you any more.” Arnold Sonunerfield 

Welcome to TEST - a visual platform to analyze thermofluids problems, verify hand 
calculations, pursue what-if scenarios, visualize thermal systems, and engage in a 
life-long journey to master thermodynamics. 

If you are new to TEST, click on the Animations tab and select a slide from Chapter 9 


Daemons RIAs Property Tables Animations Problems Tutorial Forum MyAccount Logout 


Version:! 0.505e; 


must update Jai 


lyour computer to the latest version every month (by going to www.java.com and clicking the Free Java Download button). Visit the 


Tut-nrial 

infnrmat-inn nn rljomnnc 







V 

£) Home f [J 

RIAs Property Tables 

Animations 

Problems 

Tutorial 

Forum 

MyAccount 

Logout 

Version:! 0.505e; C 


x 
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2. Click on Daemons at the menu bar at the bottom to get the following Daemons Map: 


Basic t ools 


r 


System Analysis 


1 


Unit DcsfcCfll 
Converter Cfaarfs 


Closed 


Unsteady Process 



Steady Slate Steady State 
(cycles) I 


States & Properties 


Uniform 

How 


Unsteady Process 


Generic 


Specific 

I 


Generic 


I 1 

Reciprocating Cycles HVAC/Psychrometry Comb, Equilibrium 


I I I 1 

Uniform Non-Mixing Semi-Mixing Mixing 

System Non-Uniform Non-Uniform Non-Uniform 


I 

Single-Flow 


; 


I 

Non-Mixing Multi-Flow 


I 

Mixing Multi- Flow 


Specific 


T 


T 


I 

Vapor and Gas Refrigeration HVAC Gas Combustion 

Power Cycles Cycles Psychrometry Dynamics Equilibrium 


3. We can choose the States and Properties - Uniform system to get the states 1 and 2, and 
then calculate the work for Isothermal process, OR: go to System Analysis - Closed - 
Generic - Uniform system to make the direct analysis of the process. Choosing the States & 
Properties - Uniform System, we get: 



Launches the PG Model 
Generic, Uniform System, Closed process Daemon 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest aas 


model which obeyes the ideal gas equation (pv = RT) and assumes 
)e constant. In the IG (ideal gas) model, specific heats are assumed 
:emperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 


pv = RT 
IG Model 


Examples: Air is compressed in a piston-cylinder device from a beginning-state to a 
final-state with no possibility of a phase change. For specific examples, click on the 
help icon at the bottom margin of the daemon. 
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4. We choose for the Material model: the PG Model, i.e. cp = const. Clicking on it, we get the 
following screen. Now, choose State 1, Enter pi, Tl, Voll in proper units as given in data: 


System State Daemon: Perfect Gas (PG) Model 

S M rU-JnU.im 


thermofluids.nei 


| Home <rf 
I l-* / T E S T 


States > System > PG-Model 


tv STEtffl K-H 


Mixed C si C English | Case-0 v R? Help Messages On Super-lterate J 

State Panel | 


Super-Caiculate 


< I State-1 v > 


I? pi 


fl50 

n si 


■ No '♦ Yes 


W Tl 


r vi 


r ui 


r hi 


m A 3/kg v p kJ/kg v p 

r e i r jl 


kPa v [300 

R? Veil 


«« »i r 


R? zl 


][~ KMg - f 

W Voll MM1 


kJ/kg v 


R1 


kJ/kg.K v o.O 


PI 


P ml 


kJ/kg v j p 


kJfkg v p 


0.2 


c _p1 


c vi 


kJ/kg.K v p 


kJ/kg.K v p 


m*3 


kg/kmol v 


5. Click on Calculate and state 1 is calculated: 


G Mixed 

r si r 

English <| Case-0 vl > 1 

W Help Messages On 

Super-lterate | 

Super-Calculate 

[^3 

| Super Initialize | 



State Panel 


n 



I/O Panel 




1 1 

} < | ©State-1 v > I 

HI Calculate 

I No-Plots 

v Initialize 

J 


■" No Yes 

N2 


R? pi 


W Tl 


r vi 


r ui 


r hi 


150.0 

■ kp > 

0 
d 

1 

> 

IEI 

v [059386 

” m*3/kg 

V (-87.17103 

U/kg 

V ; [1.90754 

U/kg v 

r si 


W Veil 


W zl 


r el 


1” P 

U/kg v 1 

|6.73127 

~ U/kg.K 

v [0.0 

II. 1 " 1 * 

v do 

■EI 

v| [-87.17103 

U/kg 

_vj [1.90754 


phil 


psil 


r ml 


R? Voll 


MM1 


l "" 

kJ/kg 

■m — 1 

kJ/kg 

v] [0.33678 

~ fcg 

V [ 0.2 

■IE 

_vj| 28 .° 

kg/kmol v 

R1 


c jpl 


c_v1 


kl 




10.29693 

~ 1 kJ/kg.K 

V h.0311 

U/kg.K 

v 10.73417 

kJ/kg.K 

V 1.40444 

Unit Less 



■ ' ~ — ■» 

6. Similarly, choose State 2 

, enter p2, T2 and m2. Note that we wrote T2 = Tl, 

m2 = ml. 
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7. Click on Calculate. State2 is calculated: 



8. Draw the p-V diagram: Select Plots-p-V diagram: 


State Panel 


I/O Panel 



jo. 29693 |fcJ/kg.K V j[l.0311 |u/lcg.K v] |p.73417 jldAfl.K v] [1.40444 |un*Let v] 
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And, immediately, following plot with the two states marked, is presented: 


! TEST: p-v Plot (PG-Model) 


[ResetP|ot]i Connect States | Grid | Draw ^3 ^3l E23 ^3 J lines for All States v Data Table 

| J xMin: |0.0894 xMax: [o?6 yMin: |90. 88289 yMax: |l561.0 Re-plot 3 

Use mouse/stylus to write | 

Note: You can use the zoom buttons (+,-) and drag the plot around to resize. To scribble on the plot by mouse drag, click the pen button first. 


p, kPa (log Scale) 


( 0 . 368 , 3079.446 ) 


1561.0 


90.8829 















































O 2 














































0 < 


0.0894 


v, m A 3/kg (log Scale) 


0.6 


You can format it further, connect the states, draw different lines such as p = c, v = c, T = c etc. (see 
the top line in the above screen shot), and change the axes limits too if required. In the following T = c 
is executed: 


p, kPa (log Scale) 


(0.721. 3238.392 ) 



0.0894 


v, m A 3/kg (log Scale) 


0.6 


53 


Download free eBooks at bookboon.com 


Software Solutions to Problems in Basic 
Thermodynamics: Part-ll 


Work, Heat and I Law of Thermodynamics 
applied to Closed systems 


9. Now, that States 1 and 2 are fully known, we can calculate Isothermal work by going to the 
I/O panel: Click on Super Calculate and go to I/O panel: 



Clicking on I/O panel, we get: 

# ******* ANALYST; Dr Muliya; TEST License: Professional******* 


# Solution logged at: Dec 11, 2013 9:26:29 PM 

# start of TEST-code 


States { 

State- 1: N2; 

Given: { pl= 150.0 kPa; Tl= 300.0 K; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.2 m A 3; } 
State-2: N2; 

Given: { p2= 800.0 kPa; T2= “Tl” K; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; } 

} 

# End of TEST-code 

# Evaluated States: 


# State- 1: N2 > PG-Model; 

# Given: pl= 150.0 kPa; Tl= 300.0 K; Vell= 0.0 m/s; 

# zl= 0.0 m; Voll= 0.2 m A 3; 

# Calculated: vl= 0.5939 m A 3/kg; ul= -87.171 kj/kg; hl= 1.9075 kj/kg; 

# sl= 6.7313 kJ/kg.K; el= -87.171 kj/kg; jl= 1.9075 kj/kg; 

# ml = 0.3368 kg; MM1= 28.0 kg/kmol; Rl= 0.2969 kJ/kg.K; 

# c_pl= 1.0311 kJ/kg.K; c_vl= 0.7342 kJ/kg.K; kl = 1.4044 UnitLess; 

# State-2: N2 > PG-Model; 

# Given: p2= 800.0 kPa; T2= £t Tl” K; Vel2= 0.0 m/s; 

# z2= 0.0 m; m2= “ml” kg; 

# Calculated: v2= 0.1114 m A 3/kg; u2= -87.171 kj/kg; h2= 1.9075 kj/kg; 

# s2= 6.2342 kJ/kg.K; e2= -87.171 kj/kg; j2= 1.9075 kj/kg; 

# Vol2= 0.0375 m A 3; MM2= 28.0 kg/kmol; R2= 0.2969 kJ/kg.K; 

# c_p2= 1.0311 kJ/kg.K; c_v2= 0.7342 kJ/kg.K; k2= 1.4044 UnitLess; 

# 

# Property spreadsheet starts: The following property table can be copied onto a spreadsheet (such 


as Excel) for further analysis or plots. 
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# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

150.0 

300.0 

0.5939 

-87.17 

1.91 

6.731 

# 

2 

800.0 

300.0 

0.1113 

-87.17 

1.91 

6.234 


# 

# Property spreadsheet ends 

# 

******* CALCU LATE VARIABLES: Type in an expression starting with an c =’ sign ( c = mdotl*(h2-hl)’ 
£ = sqrt(4*Al/PI)’ etc.) and press the Enter key)********* 

# 

#Isothermal Work done: W = pl*Voll * in (Vol2/Voll) 

=pl*Voll * In (Vol2/Voll) 

pl*Voll * in (Vol2/Voll) = -50.21929300715014 kj = -50.22 kj .... Ans. 

# Alternatively: 

=m 1 *R1 *T 1 *ln(vol2/vol 1 ) 

ml*Rl*Tl*ln(vol2/voll) = -50.21929300715015.... same as above. 
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Alternatively: 

Use the System Analysis - Closed - Generic - Uniform System Daemon for the Process analysis. 

This is direct method, and is preferable since in addition to work and heat, it calculates exergy and Tost 
work 5 too. 


1. Select the appropriate daemon for process analysis as shown below: 


Ba^ic Tools 


I 

System Analysis 

! 


States & Pr operties 


Open 


I 

Unsteady Process 


Closed 

1 

Steady State Steady State 
(cycles) I 


1 

Uniform 

System 


1 

Uniform 

Flow 


1 

Unsteady Process 


Generic 


r 


Specific 
I 


Generic 


1 


Reciprocating Cycles HVAC/Psychromeity Comb , 3 Equilibrium 


T 


Uniform 

System 


I 

Non-Mixing 

Non-Uniform 


I 

Semi-Mixing 

Non-Uniform 


1 

Mixing 

Non-Uniform 


1 

Single-Flow t 


1 

Non-Mixing Multi-Flow 


1 

Mixing Multi-Flow 


I 

Specific 


i 

Vapor and Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 

Psychrometry 


I 

Gas 

Dynamics 


Combustion 
3 Equilibrium 


2. Clicking on Uniform System, and choosing the Perfect Gas Model with cp = const, gives 
following window. Fill up the known parameters viz, pi, Tl, Voll for State 1, and click on 
Calculate. We get: 


Generic, Uniform-System, Closed Process Daemon: PG Model 


thermofluids.net > Daemons > Systems > Closec 

1 > Process > Generic : 

> Unifom 

n > PG-Model 

mm m w of 

1 

i '-I i 


c„ = const 

I 


Move mouse over a variable to display its value 




1 











Mixed 

r si r 

English | ©Case-0 ^ > 

K? Help Messages On 

Super-lterate 

Super-Calcidate 

ESI 

| Super-lmtialize ■ 

State Panel 

1 

Process 



j 

Exergy Panel 

1 



1 












ip | ©State- 

, r-] 


^ ] No-Plots 

3 

Initialize 


Formation Enthalpy: 

No • Yes 

N2 



1 ? pi 


W Tl 


r 

vl 


r Ul 




150.0 

■ >pg 

V 300.0 

mz 

V 10.59386 

[ rTt'S/kg 

v |-87.17103 

kJ/kg 

V jl .90754 

kJ/kg v | 

r si 


I* Veil 


I? 

zl 


r el 


r ji 


|6 73127 

kJ/kg.K 

v go 

■ r* 

Z °o 


■ft 

v |-87.17103 

kJ/kg 

v |l .90754 

kJ/kg v | 

phil 


psil 


r 

ml 


W Voll 


MM1 


1 

~~\kJ/kg j 

- r 

kJ/kg 

V |0. 33678 

"SZ 

V |) 2 

■1 ^ 

_vj | 28.0 

kg/kmol v- 


c _p1 


kJ/kg.K V 11.0311 


c_v1 kl 

|0. 73417 kJ/kg.K v 1 1.40444 UnitLess 
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3. Fill up known parameters for State 2, click on Calculate: 



4. Go to Process Panel. T_B is already checked there; also check W_0 (i.e. other work) as zero. 
Click on Calculate and get the following results: 


lay its value with more precision. 


Mixed r SI English < | ©Case-0 v] > 
State Panel I Process Panel 


17 Help Messages On Super-lterate j 


Super-Calculate 


Process-A [1-2] v > | 

Bf State-1 vj 

|j State-2 v j 

Initialize 

Hp^constant 

Q 


r w_b 


(7 W_0 


17 TJB 


1-50.21929 


v ' |-50.21929 

["H 

o 

o 

> 


v 298.15 

■l?Z ZM 

r S_gen 


r n 


Delta_E 


Delta_S 


10.00104 

kJ/K 

V [To 

UnitLess 

v jo.O 

fkJ 

V 1-0.1674 

U/K v 



Note that we get: W_B = Boundary work for this Isothermal process as -50.22 kj; 

Also, the Heat rejected Q = W_B = -50.22 kj for Isothermal process, as it should be. 

5. To make exergy analysis, we should first choose ‘dead state’. This is State ‘O’, take it as 
pO = 100 kPa, T_0 = 25 C. Go to States Panel and fill up these parameters for State 0, 
and click on Calculate: 


Move mouse over a variable to display its value 


<*■ Mixed 

r si 

C English < f 

©Case-0 v j > | 


!7 Help Messages On 

Super-lterate 

Super-Calculate 

| t nad | 

| Super-lnitialize | 


State Panel 

Process 






L 




[©State-0 v 


|j No-Plots 

V 

■ 

Initialize 



1' No Yes 

N2 

B1 

17 pO 

17 TO 



r 

vO 


r uO 


r ho 


jioo.o 

p kPa v |25.0 

■ |deg-C 

V 

0.88529 

m*3/kg 

v | [-88.52925 

kJ/kg 

o 

© 

kJ/kg v 

r so 


17 VelO 



(7 

zO 


r e0 


r jo 


[6.84529 

kJ/kg 

|.K v |o.O 

II m,s 

V 

W 


|[m_ 

v [-88.52925 

| kJ/kg 

v [o ~o 

kJ/kg v 

phiO 
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mO 


r voio 
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V 
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6. Now, go to Exergy Panel: 



Observe that I = TO * AS = 0.30969 is the ‘lost work’ 

7. Click on SuperCalculate to generate the TEST Code and to get a record of all calculations, 
and go to I/O Panel to get TEST code etc: 


# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 
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States { 

State-0: N2; 

Given: { p0= 100.0 kPa; T0= 25.0 deg-C; VelO= 0.0 m/s; z0= 0.0 m; } 

State- 1: N2; 

Given: { pl= 150.0 kPa; Tl= 300.0 K; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.2 m A 3; } 

State-2: N2; 

Given: { p2= 800.0 kPa; T2= “Tl” K; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; } 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; } 

} 

# p; n d of TEST-code 

******* DETAILED OUTPUT: All the computed properties and variables are displayed on this 
block ************* 

# Evaluated States: 

# 

# State-0: N2 > PG-Model; 

# Given: p0= 100.0 kPa; T0= 25.0 deg-C; Vel0= 0.0 m/s; 

# z0= 0.0 m; 

# Calculated: v0= 0.8853 m A 3/kg; u0= -88.5292 kj/kg; h0= 0.0 kj/kg; 

# s0= 6.8453 kJ/kg.K; e0= -88.5292 kj/kg; j0= 0.0 kj/kg; 

# phi0= 0.0 kj/kg; psi0= 0.0 kj/kg; MM0= 28.0 kg/kmol; 

# R0= 0.2969 kJ/kg.K; c_p0= 1.0311 kJ/kg.K; c_v0= 0.7342 kJ/kg.K; 

# k0= 1.4044 UnitLess; 

# State- 1: N2 > PG-Model; 

# Given: pl= 150.0 kPa; Tl= 300.0 K; Vdl= 0.0 m/s; 

# zl= 0.0 m; Voll= 0.2 m A 3; 

# Calculated: vl= 0.5939 m A 3/kg; ul= -87.171 kj/kg; hl= 1.9075 kj/kg; 

# sl= 6.7313 kJ/kg.K; el= -87.171 kj/kg; jl = 1.9075 kj/kg; 

# phil= 6.2086 kj/kg; psil= 35.9014 kj/kg; ml= 0.3368 kg; 

# MM1= 28.0 kg/kmol; Rl= 0.2969 kJ/kg.K; c_pl= 1.0311 kJ/kg.K; 
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# c_vl= 0.7342 kJ/kg.K; kl= 1.4044 UnitLess; 

# State-2: N2 > PG-Model; 

# Given: p2= 800.0 kPa; T2= “Tl” K; Vel2= 0.0 m/s; 

# z2= 0.0 m; m2= “ml” kg; 

# Calculated: v2= 0.1114 m A 3/kg; u2= -87.171 kj/kg; h2= 1.9075 kj/kg; 

# s2= 6.2342 kJ/kg.K; e2= -87.171 kj/kg; j2= 1.9075 kj/kg; 

# phi2= 106.1536 kj/kg; psi2= 184.0973 kj/kg; Vol2= 0.0375 m A 3; 

# MM2= 28.0 kg/kmol; R2= 0.2969 kJ/kg.K; c_p2= 1.0311 kJ/kg.K; 

# c_v2= 0.7342 kJ/kg.K; k2= 1.4044 UnitLess; 

# Property spreadsheet starts: The following property table can be copied onto a spreadsheet (such 

as Excel) for further analysis or plots. 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 0 100.0 298.2 0.8853 -88.53 0.0 6.845 

# 1 150.0 300.0 0.5939 -87.17 1.91 6.731 

# 2 800.0 300.0 0.1113 -87.17 1.91 6.234 

# 

# Property spreadsheet ends 

# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f-State = State-2; 

# Given: W_0= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= -50.21929 kj; W_B= -50.21929 kj; S_gen= 0.0010386907 kJ/K; n= 1.0 UnitLess; 

# Delta_E= -0.0 kj; Delta_S= -0.16739765 kJ/K; 

# Exergy Analysis Results: 

# Exergy Analysis for Process - A (Dead state: State-0) 

# Given: Q= -50.21929 kj; T_0= 298.15 K; Q_l= 0.0 kj; 

# T_l= 298.15 K; 

# Calculated: Delta_Phi= 33.65961 kj; W_u= -33.96929 kj; 1= 0.30969 kj; 

# S_gen.univ= 0.00104 kJ/K; W_rev= -33.65961 kj; W= -50.21929 kj; 

# W_atm= -16.25 kj; Q_0= -50.21929 kj; 
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Prob.4.19. A mass of 1.2 kg of Air at 150 kPa and 12 C is contained in a gas-tight friction-less piston- 
cylinder device. The air is now compressed to a final pressure of 600 kPa. During the process heat is 
transferred from air such that the temp inside the cylinder remains constant. Calculate the work done 
during this process. [Ref. 1]. 


TEST Solution: 

Use the System Analysis - Closed - Generic - Uniform System Daemon for the Process analysis . 
Following are the steps: 

1. Select the appropriate daemon for process analysis as shown below: 
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2. Clicking on Uniform System, and choosing the Perfect Gas Model with cp = const, gives 
following window. Select the gas as N2, Fill up the known parameters viz, pi, Tl, Voll for 
State 1, and click on Calculate. We get: 


Generic, Uniform-System, Closed Process Daemon: PG Model 
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3. Select State 2, enter the known parameters, i.e. p2, T2, m2. Click on Calculate: 



4. Go to Process Panel, enter W_0 (i.e. works other than boundary works) as zero, 
click on Calculate: 



Note that for this Isothermal process, Boundary work, W_0 is calculated as -140.85 kj. (Ans.) 


Negative work indicates work done on the system. Obviously, heat transfer Q is equal to W_B and is 
negative, i.e. heat is leaving the system in this Isothermal process. 
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5. Plot below shows the States 1 and 2 on a T-s diagram: 
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6. Click on SuperCalculate to produce the TEST code, (with which we can regenerate these 
calculations later by loading this TEST code in the I/O Panel and clicking SuperCalculate) 
Now, go to I/O panel to view the TEST code and other calculated States. Only part of the 
I/O output is shown below: 

# OUTPUT OF SUPER-CALCULATE: ( 

# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: N2; 

Given: { pl= 150.0 kPa; Tl= 12.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; ml= 1.2 kg; } 

State-2: N2; 

Given: { p2= 600.0 kPa; T2= “Tl” deg-C; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; } 

} 


# p n d of TEST-code 

# Property spreadsheet starts:. — 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

150.0 

285.2 

0.5645 

-98.07 

-13.4 

6.679 

# 

2 

600.0 

285.2 

0.1411 

-98.07 

-13.4 

6.267 


Prob. 4.20. A vessel having a volume of 5 m A 3 contains 0.05 m A 3 of sat. liquid water and 4.95 m A 3 of 
sat. water vapour at 0.1 MPa. Heat is transferred until the vessel is filled with sat. vapour. Determine the 
heat transfer for this process. [Ref:2] 

TEST Solution: 

We shall use the States & Properties - Uniform System, with PC (i.e. Phase Change) Material model 
for Water, and then calculate the heat transferred in the I/O panel, using it as a calculator: 
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Following are the steps: 

1. Select the System State daemon: 
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2. Choose the PC model for Material model. Following daemon presents itself. We shall call 
the sat. liq. As State 1, sat. vapour as State 2, and the combined liq + vapour as State 3. So, 
Fill up the known parameters pi, xl and Vollfor State 1: 


System State Daemon: Phase-Change (PC) Model 
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Note that in the above dryness fraction xl is zero since it is sat. liq. state. Now, click on Calculate and 
we get: 
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Observe that ml, T1 etc are immediately calculated for Sat. water at 100 kPa. 

3. Now, enter known parameters, i.e. p2, x2 and Vol2 for State 2: 
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Here, x2 = 1 since we are dealing with sat. vapour. Click on Calculate, and we get: 



Note that m2, u2, h2 etc. are immediately calculated for sat. water vap at 100 kPa. 


4. Now, enter State 3. This is when the entire tank is filled with sat. water vap. i.e. x3 = 1, and of 
course, Vol3 = Voll + Vol2, and total mass m3 = ml+m2. Enter these parameters for State 3: 
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5. Click on Calculate and we get: 



6. Get the T-s plot where States 1, 2 and 3 are shown: 
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7. Click on SuperCalculate to produce the TEST code and other calculated results. Go to I/O 
panel to see them. Part of the output is shown below: 

# OUTPUT OF SUPER-CALCULATE: 


# Daemon Path: States>System>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: H20; 

Given: { pl= 100.0 kPa; xl= 0.0 fraction; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.05 m A 3; } 
State-2: H20; 

Given: { p2= 100.0 kPa; x2= 1.0 fraction; Vel2= 0.0 m/s; z2= 0.0 m; Vol2= 4.95 m A 3; } 
State-3: H20; 

Given: { x3= 1.0 fraction; Vel3= 0.0 m/s; z3= 0.0 m; m3= “ml+m2” kg; Vol3= “voll+vol2” m A 3; } 

} 

# p n d of TEST-code 
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# Property spreadsheet : 


# State 

p(kPa) 

T(K) 

X 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

100.0 

372.8 

0.0 

0.001 

417.34 

417.44 

1.303 

# 02 

100.0 

372.8 

1.0 

1.694 

2506.06 

2675.46 

7.359 

#03 

2026.13 

486.2 

1.0 

0.0983 

2600.47 

2799.74 

6.336 


# 


# >f>f>f>f > f>f CALCULATE VARIABLES: Type in an expression starting with an £ =’ sign ( c = mdotl*(h2-hl)’, 
£ = sqrt(4*Al/PI)’, etc.) and press the Enter key)********* 

# Calculate the heat transferred in the I/O panel, using it as a calculator: 

Q = [(m3 * u3) - (ml * ul + m2 * u2)] heat transferred, since it is at constant volume 

ml*ul + m2*u2 = 27329.407721629643 kj 
m3*u3 = 132261.66395801632 kj 
#Therefore: 

Q = m3*u3 - (ml*ul + m2*u2) = 104932.25623638667 kj = 104932.26 kj .... Ans. 

In addition, note that the masses of sat. liq. and vapour are: 
ml = 47.938637362598115 = 47.94 kg .... Mass of sat. liq. 
m2 = 2.922078117838602 kg = 2.92 kg .... Mass of sat. vap. 

And, total mass m3 is: 

m3 = 50.86071548043672 = 50.86 kg .... Total mass 


Prob.4.21. A cylinder fitted with a piston has a volume of 0.1 m A 3 and contains 0.5 kg of steam at 0.4 
MPa. Heat is transferred to the steam until the temp is 300 C, while the pressure remains constant. 
Determine the heat transfer and work for this process. [Ref: 2] 
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TEST Solution: 

We use the System Analysis - Closed - Generic - Uniform System Daemon for the Process analysis. 


Following are the steps: 

1. Select the appropriate daemon for process analysis as shown below: 
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2. Clicking on Uniform System, choose the PhaseChange (PC) Model for Material Model since 
we are dealing with Steam/Water. 


Generic, Uniform System, Closed Process Daemons: Select a Material Model 




4k Home of 
I.J TEST 






^ y 




mm 


Select a material model to launch the daemon. 


\£L 

PC Model 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (1-120* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

Working fluids such as H20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 
possibility of a phase transition. 


Examples: R-134a vapor is compressed in a piston-cylinder device from a beginning-state to a finai-state. To 
find the work transfer if compression is assumed isentropic. For specific examples, click on the help icon at the 
bottom margin of the daemon. 
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3. Fill up the known parameters viz, pi, ml, Voll for State 1. Click on Calculate. We get: 

Generic, Uniform-System, Closed Process Daemon: PC Model 




4. Select State 2, enter known parameters, i.e. p2, T2, m2, and click on Calculate: 



5. Go to Process Panel, enter b-state and f-state, enter W_0 = 0 (i.e. works other than pdV 
work), and click on Calculate. We get: 
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Thus: Boundary work, W_B = 90.9656 kj and the heat transferred Q = 771.41 kJ....Ans. 

Note that Work is positive, i.e. work done by the system. 
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Heat transfer q is positive, i.e. Heat transferred into the system. 


6. On a T-s diagram, the State points are shown as follows: 
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7. Click on SuperCalculate to produce the TEST code, calculated State properties etc. Go to 
I/O panel to see the code. Part of I/O output is shown below: 

# 

# Daemon Path: Systems>CIosed>Process>Generic>Uniform>PC-Model; v-10.bb06 

# 

# start of TEST-code 

States { 

State- 1: H20; 

Given: { pl= 400.0 kPa; Vell= 0.0 m/s; zl= 0.0 m; ml= 0.5 kg; Voll= 0.1 m A 3; } 
State-2: H20; 

Given: { p2= “pi” kPa; T2= 300.0 deg-C; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; } 

} 

Analysis { 

Process- A: b-State = State- 1; f-State = State-2; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; } 

} 

# p n d of TEST-code 

# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f-State = State-2; 

# Given: W_0= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 771.41223 kj; W_B= 90.9656 kj; S_gen= -0.7952037 kJ/K; 

# Delta_E= 680.44666 kj; Delta_S= 1.7921257 kJ/K; 


Prob.4.22. Air at 1.02 bar, 22 C, initially occupying a cylinder volume of 0.015 m3, is compressed 
reversibly and adiabatically to a pressure of 6.8 bar. Calculate: (i) Final volume (ii) Final temp, and 
(iii) Work done. [Ref: 4] 

(b) In addition: If State 2 is reached by a polytropic process (n = 1.3) instead of by isentropic process, 
find out the values of Work and Heat transfers and their direction. 
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TEST Solution: 

We use the System Analysis - Closed - Generic - Uniform System Daemon for the Process 
analysis. 

Following are the steps: 

1. Select the appropriate daemon for process analysis as shown below: 
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2. Clicking on Uniform System, choose the Permanent Gas (PG) Model for Material Model 
since we are dealing with Air. Enter parameters pi, T1 and Vollfor State 1, click on 
Calculate. We get: 


Generic, Uniform-System, Closed Process Daemon: PG Model 
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3. Select State 2, enter p2, m2 = ml, and s2 = si since it is an isentropic (i.e. reversible, 
adiabatic) process. Click on Calculate. We get: 
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4. Now, go to Process Panel, enter b_state = State 1, f- state = State 2, and Q = 0 since it is 
adiabatic process; click on Calculate. We get: 



Thus: 

Final volume, Vol2 = 0.00387 m A 3, Final temp, T2 = 234.62 C.... Ans. 

W_B = boundary work = -2.752 kj .... Ans. Negative sign means that work is done on the system. 

1. Plot the States 1 and 2 on the T-s diagram: 


t, k 
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558.55 


265.64 

6.18 s, kJ/kg.K 7.56 



} 2 
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(b) If State 2 is reached by a polytropic process (n = 1.3), what are the values of Q and W_B? 

Let the state after the polytropic process be designated as State 3. Note that State 2 and State 3 are identical: 


1. Select State 3 and enter p3, T3 and ml. These are essentially the same as for State 2. 
Click on Calculate. We get: 
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2. Now, go to Process Panel. Enter b-state = State 1, f- state = State 3, n = 1.3 and Other Works, 
W_0 = 0. Click on Calculate. We get: 



Note that Q = -0.922 kj, W_B = -3.674 kj Ans. Negative sign means: Q leaving the system, 

W_B done on the system. 

3. Click on SuperCalculate to generate TEST code and get all calculated results. See them on 
the I/O panel. Part of the output is given below: 


# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { pl= 102.0 kPa; Tl= 22.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.015 m A 3; } 
State-2: Air; 

Given: { p2= 680.0 kPa; s2= “si” kJ/kg.K; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; } 
State-3: Air; 

Given: { p3= <c p2” kPa; T3= CC T2” K; Vel3= 0.0 m/s; z3= 0.0 m; m3= “ml” kg; } 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-3; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; n= 1.3 UnitLess; } 

} 


# p n d of TEST-code 

# 
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#DETAILED OUTPUT: 

# Evaluated States: 

# State- 1: Air > PG-Model; 

# Given: pl= 102.0 kPa; Tl= 22.0 deg-C; Vell= 0.0 m/s; 

# zl= 0.0 m; Voll= 0.015 m A 3; 

# Calculated: vl= 0.8304 m A 3/kg; ul= -87.7146 kj/kg; hl= -3.0105 kj/kg; 

# sl= 6.8709 kJ/kg.K; el= -87.7146 kj/kg; jl= -3.0105 kj/kg; 

# ml = 0.0181 kg; MM1= 28.97 kg/kmol; Rl= 0.287 kJ/kg.K; 

# c_pl= 1.0035 kJ/kg.K; c_vl= 0.7165 kJ/kg.K; kl= 1.4005 UnitLess; 

# State-2: Air > PG-Model; 

# Given: p2= 680.0 kPa; s2= “si” kJ/kg.K; Vel2= 0.0 m/s; 

# z2= 0.0 m; m2= “ml” kg; 

# Calculated: T2= 507.7737 K; v2= 0.2143 m A 3/kg; u2= 64.6319 kj/kg; 

# h2= 210.3562 kj/kg; e2= 64.6319 kj/kg; )2= 210.3562 kj/kg; 

# Vol2= 0.0039 m A 3; MM2= 28.97 kg/kmol; R2= 0.287 kJ/kg.K; 

# c_p2= 1.0035 kJ/kg.K; c_v2= 0.7165 kJ/kg.K; k2= 1.4005 UnitLess; 

# State-3: Air > PG-Model; 

# Given: p3= “p2” kPa; T3= “T2” K; Vel3= 0.0 m/s; 

# z3= 0.0 m; m3= “ml” kg; 

# Calculated: v3= 0.2143 m A 3/kg; u3= 64.6319 kj/kg; h3= 210.3562 kj/kg; 

# s3= 6.8709 kJ/kg.K; e3= 64.6319 kj/kg; j3= 210.3562 kj/kg; 

# Vol3= 0.0039 m A 3; MM3= 28.97 kg/kmol; R3= 0.287 kJ/kg.K; 

# c_p3= 1.0035 kJ/kg.K; c_v3= 0.7165 kJ/kg.K; k3= 1.4005 UnitLess; 

# 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 102.0 295.2 0.8304 -87.71 -3.01 6.871 

# 2 680.0 507.8 0.2143 64.63 210.36 6.871 

# 3 680.0 507.8 0.2143 64.63 210.36 6.871 

# 

# Property spreadsheet ends 

# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f-State = State-3; 

# Given: W_0= 0.0 kj; T_B= 298.15 K; n= 1.3 UnitLess; 

# Calculated: Q= -0.92218256 kj; W_B= -3.6739995 kj; S_gen= 0.0030930154 kJ/K; 
Delta_E= 2.751817kJ; 

# Delta_S= -0.0 kJ/K; 
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“Prob.4.23. 5 kg of Nitrogen at 100 C is heated in a reversible, non-flow, constant volume process till 
the pressure becomes three times the initial pressure. Determine: (i) final temp (ii) change in internal 
energy (iii) change in enthalpy, and (iv) heat transfer. Take R = 0.297 kJ/kg.K, cv = 0.7435 kJ/kg.K. 
[VTU-Jan. 2004]” 

Note that this is the same as Prob. 4.8 which was solved with EES. 

Now, let us solve it with TEST: 

TEST Solution: 

We use the System Analysis - Closed - Generic - Uniform System Daemon for the Process analysis. 
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Following are the steps: 

1. Select the appropriate daemon for process analysis as shown below: 
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2. Clicking on Uniform System, choose the Permanent Gas (PG) Model for Material Model 

since we are dealing with Air. Enter parameters ml, T1 and Voll (= 1 m A 3.... assumed) for 
State 1, click on Calculate. We get: 


Generic, Uniform-System, Closed Process Daemon: PG Model 
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Note that, in calculations, we will be using the built-in properties for R, cp and cv, as seen in the 
above screenshot. 
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3. Select State 2, enter p2 = 3 * pi, Vol2 = Voll, m2 = ml. Click on Calculate. We get: 
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4. Go to Process Panel, enter W_B = 0 since it is const, vol. process, W_0 = o, since there is 
no other work interaction. Click on Calculate. We get: 



5. States 1 and 2 are shown in the p-V diagram: 


p, kPa (log Scale) 
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v, m A 3/kg (log Scale) 
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Thus: 

Final temp, T2 = 846.3 C = 1119.4 K .... Ans., 

Heat transfer, Q = 2739.57 kj .... Ans. Work is done by the system. 

6. Click on SuperCalculate to get TEST code and calculated results. Also, calculate the change 
in internal energy and enthalpy in the I/O panel. Go to I/O panel. Part of the output is:: 

# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: N2; 

Given: { Tl= 100.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; ml= 5.0 kg; Voll= 1.0 m A 3; } 

State-2: N2; 

Given: { p2= “3*pl” kPa; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; Vol2= “Voll” m A 3; } 

} 
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Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_B= 0.0 kj; W_0= 0.0 kj; T_B= 298.15 K; } 

} 


# End of TEST-code 

# Property spreadsheet starts: 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 553.99 373.2 0.2 -33.47 77.33 6.568 

# 2 1661.98 1119.4 0.2 514.45 846.84 7.375 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Process-A: b-State = State- 1; f-State = State-2; 

# Given: W_B= 0.0 kj; W_0= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 2739.57 kj; S_gen= -5.155695 kJ/K; Delta_E= 2739.5671 kj; Delta_S= 4.0328584 kJ/K; 

#******CALCULATE VARIABLES: Type in an expression starting with an “=’ sign (‘= mdotl*(h2-hl)’, 
‘= sqrt(4*Al/PI)’, etc.) and press the Enter key)********* 

# Change in Internal Energy: AU = ml * (u2 - ul) 

i.e. AU = ml * (u2 - ul) = 2739.5670715475585 = 2739.57 kj... Ans. 

# Change in Enthalpy: AH = ml * (h2 - hi) 

i.e. AH = ml * (h2 - hi) = 3847.5560358332727 = 3847.56 kj .... Ans. 


Compare the above results with those obtained with EES: 


Unit Settings; SI K kPa kJ molar deg 

cp = 1 041 [J/kg-K] cv = 743.5 [J/kg-K] 

AU = 2.773E+06 [J] I m = 5 [kg] 


Q = 2.773E+06 [J] 
T2 =1119 [K] 


R = 297 [J/kg-K] 
W = 0 [J] 


AH = 3.831 E+06 [J] 
PressureRatio = 3 
T1 = 373 [K] 


It is observed that results match very well. 
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Prob.4.24.A piston-cylinder device contains 0.8 kg of Nitrogen initially at 100 kPa and 27 C. The nitrogen 
is now compressed slowly in a polytropic process (P.V A 1.3 = const.) until the volume is reduced by one- 
half. Determine the work done and the heat transfer. [Ref: 1] 


TEST Solution: 

We use: System Analysis - Closed - Generic - Uniform System Daemon for the Process analysis. 


Following are the steps: 

1. Select the appropriate daemon for process analysis as shown below: 
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2. Clicking on Uniform System, choose the Permanent Gas (PG) Model for Material Model since 
we are dealing with N2. Enter parameters pi, Tl, ml for State 1, click on Calculate. We get: 


Generic, Uniform-System, Closed Process Daemon: PG Model 
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3. Enter known quantities for State 2. We have: m2 = ml and Vol2 = 0.5 * Voll. Click on 
Calculate, but the entered data is not sufficient to make all calculations: 
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4. Let us proceed to the Process Panel and enter n = 1.3 (i.e, polytropic index), Other Works, 
W_0 = 0. Click on Calculate. We get: 



5. Since iteration has to be done with reference to other states, we have to click on 
SuperCalculate to complete the calculations. Then, we get: 
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Thus: 

Work done, W_B = -54.93 kj.... Ans. (Negative sign means work done on the system) 


Heat transfer, Q = -14.19 kj Ans. (Negative sign means heat rejected by the system). 
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6. Now, go back to State Panel and examine State 2: 



We see that values of p2 and T2 are now posted for State 2. 

Thus: 

P2 = 246.23 kPa, T2 = 96.38 C .... Ans. 

7. T-s diagram showing States 1 and 2 is easily obtained: 
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8. Go to I/O panel to see the TEST Code and the calculated values: 


# OUTPUT OF SUPER- CALCULATE (starts from your inputs) 

# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: N2; 

Given: { pl= 100.0 kPa; Tl= 27.0 deg-C; Vdl= 0.0 m/s; zl= 0.0 m; ml= 0.8 kg; } 
State-2: N2; 

Given: { Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; Vol2= “0.5*Voll” m A 3; } 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; n= 1.3 UnitLess; } 

} 

# End of TEST-code 


# 
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# Property spreadsheet starts: The following property table can be copied onto a spreadsheet (such 

as Excel) for further analysis or plots. 

# 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

100.0 

300.2 

0.8912 

-87.06 

2.06 

6.852 

# 

2 

246.23 

369.5 

0.4456 

-36.13 

73.6 

6.799 


# Property spreadsheet ends 

# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f- State = State-2; 

# Given: W_0= 0.0 kj; T_B= 298.15 K; n= 1.3 UnitLess; 

# Calculated: Q= -14.185789 kj; W_B= -54.93416 kj; 

# S_gen= 0.0050608176 kJ/K; Delta_E= 40.748367 kj; 

# Delta_S= -0.042518552 kJ/K; 


Prob.4.25. A quantity of air at a pressure of 100 kPa, 27 C occupying a volume of 0.5 m A 3 is compressed 
to a pressure of 500 kPa and volume of 0.12 m A 3 according to the law pv A n = const. Find: (i) the value 
of index n (ii) the mass of air (iii) work transfer (iv) heat transferred during the process, and (v) change 
in entropy. [VTU-BTD-July 2007] 
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TEST Solution: 

Following are the steps: 

1. Select System Analysis-Generic-Uniform System: 


liasic Tools 




\ 

System Analysis 

, I . 


States 8t Properties 


Open 


Unsteady Process 


Closed 

1 

Steady State Steady State 
(cycles) I 


I 

Uniform 

System 


1 

Uniform 

Flow 


1 

Unsteady Process 


Generic 


: 


Specific 


Generic 


T 


1 


Reciprocating Cycles HVAC/Psychro metry Comb, S Equilibrium 


T 


Uniform 

System 


I 

Non-Mixing 
Non l In if arm 


I 

Semi-Mixing 

Non-Uniform 


1 

Mixing 

Non-Uniform 


I 

Single-Flow 


I 

Non-Mixing Multi-Flow 


1 

Mixing Multi- Flow 


Specific 


r~ 

Vapor and Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 

Psychrometry 


I 

Gits 

Dynamics 


Combustion 
S Equilibrium 


2. For Material Model, select ‘Perfect Gas’ (PG) Model: 


>s* 

S3 

SL Model 


solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 


Examples: A block of copper is heated from a beginning-state to a final-state. To find the heat transfer 
necessary for the process. For specific examples, click on the help icon at the bottom margin of the daemon. 



RIA: SL Process 
Simulator 


Process Simulation by RIAs: These rich internet applications (RIAs) can be used to interactively explore a 
process where a system goes from a beginning state to a final state under some given constrain (say, the 
temperature remains constant) Unlike the daemons, the RIAs do not require a thorough thermodynamic 
background and can be used to to gain practical insight alongside learning the underlying theory. 

Examples: Watch the temperature rise as a block of copper is heated from a beginning-state to a final-state. 
For specific examples, click on the help icon at the bottom margin of the daemon. 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrog en, and air are examples of a pure gas. The PG (perfect gas) model is 


Cp = cons] Launches the pg Model as model which obeyes the ideal gas equation (pv = RT) and assumes 

I Generic, Uniform System, Closed process Daemon^ be constant> In the IG (jdeal gas ) mo del, specific heats are assumed 
PG Model M 1 1 


pvJtRT 


| to be function of temperature only. The RG (real gas) model uses generalized 
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3. State 1: Enter pi, Tl, Voll. Hit Enter (or click Calculate). 

Generic, Uniform-System, Closed Process Daemon: PG Model 


Move mouse over a varia 


J 


t Home of 
TEST 


um 


i i i 



n > PG-Model 1 


□ 

O 

□ 

□ 

. j |c„ = const j j 


mm 












G Mixed 

r si 

C English <[ 

Case-0 v > | 

f? Help Messages On 

Super-lterate | 

Super-Calculate 

ESI 

Super-lnitialize H 


State Panel 



1 

Exergy Panel 

: I 


- i 

1 1 

j < | ©State-1 v 

Jjcalculate 

1 1 No-Plots 

v | Initialize 


Formation Enthalpy: 

No ♦ Yes 


- 

i? P i 


W Tl 


r Vi 


r ui 


r hi 


100.0 

n kp a v {270 

| deg-C 

vj |o. 86139 

m*3/kg 

V 1-84 13202 

U/kg 

jvJ [2.00699 

kJ/kg v 

1 r si 


W Veil 


W Zl 


r el 


r jf 


[6.8934 

[idrtcg. 

K v |o.0 


v jo ,0 

|l«L_ 

V [-84.13202 

U/kg 

_vj [2.00699 

U/kg v 

phil 


psil 


r ml 


W Voll 


MM1 


1 

> 

f 

kJ/kg 

_vj [0.58046 

kg 

Z3I°5 

■ks' 

_vj [28.97 

kg/kmol v 

p- 


c _p1 


C_v1 


kl 




[0.28699 

U/kg.K v 1 1.00349 

j U/kg.K 

V 0.71651 

U/kg.K 

vj|l,40054 

UnitLess 

3 













4. State 2: Enter p2, Vol2, and m2 = ml. Hit Enter (or click Calculate). 


Move mouse over a variable to display its value with more precision. 


Mixed r SI English jcJ Case-0 v [✓ Help Messages On Super-lterate | 


Super-Calculate 


Super-lnitialize 


V No Yes Ai 


W p2 


r 72 


1 


r u2 


r h 2 



|500.0 

■ 

V [87.03 

deg-C 

V |0.20673 

jirY^g 

r| [41.12008 

U/kg 

v |62.24673 

U/kg 

v | 

r s2 


W Vel2 


C\J 

N 

, 11 


r e2 


T j2 



[6.61447 

U/kg.K 

v 0.0 

||m/s 

V jo.o 

HI* 

v ; [41.12008 

U/kg 

v [62.24673 

U/kg 

J 

phi2 


psi2 


W m2 


W Vol2 


MM2 



r~ 

kJ/kg 

3r 

kJ/kg 

V 1 I=nn1 

■S 

vpO.12 

■*» 

v [28.97 

kg/kmol 

V 


R2 c _p2 cj/2 k2 

|o. 28699 IfcJAflX "vj [l .00349 [ U/kg.K "v] |o.71651 | fcJ /kg.K ~v ] |l .40054 [ UniLcsa v 


5. Go to Process Panel. Enter b-state and f-state. Click W_0 and enter W_0= 0. Click on Calculate: 
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6. Now, click on SuperCalculate. Go to States Panel and see: 




L / i 

Maastricht University 


Join the best at 
the Maastricht University 
School of Business and 
Economics! 

i i: r ji 

• 33 rd place Financial Times worldwide ranking: MSc 
International Business 

• 1 st place: MSc International Business 

• 1 st place: MSc Financial Economics 

• 2 nd place: MSc Management of Learning 

• 2 nd place: MSc Economics 

• 2 nd place: MSc Econometrics and Operations Research 

• 2 nd place: MSc Global Supply Chain Management and 
Change 

Sources: Keuzegids Master ranking 20 1 3; Elsevier 'Beste Studies' ranking 2012; 
Financial Times Global Masters in Management ranking 2012 


Maastricht 
University is 
the best specialist 

university in the 
Netherlands 

(Elsevier) 

Master's Open Day: 22 February 201 4 

www.mastersopenday.nl 
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And State 2: 



Thus: 

Index, n = 1.12776 

Mass of air = ml = m2 =0.58046 kg 

Work transfer = W_B = -78.27 kj (Work done on the system, therefore negative) 
Heat transfer = -53.3081 kj 

Entropy change = (s2 - si) = -0.2789284256081599 = -0.2789 kJ/kg.K 
Total change in entropy of system = Delta_S = -0.16191 kJ/K .... Ans. 


T_s diagram: 


T, K 


C 5.964, 397.687 ) 


396.2 


270.14 

5.95 



s, kJ/kg.K 


7.53 
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And the I/O panel shows: 

# Daemon Path: Systems>Closed>Process>Generic>Uniform>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { pl= 100.0 kPa; Tl= 27.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.5 m A 3; } 
State-2: Air; 

Given: { p2= 500.0 kPa; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; Vol2= 0.12 m A 3; } 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_0= 0.0 kj; T_B= 298.15 K; } 

} 


# End of TEST-code 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 100.0 300.2 0.8614 -84.13 2.01 6.893 

# 2 500.0 360.2 0.2067 -41.12 62.25 6.614 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f-State = State-2; 

# Given: W_0= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= -53.3081 kj; W_B= -78.27469 kj; S_gen= 0.016890267 kJ/K; 

n= 1.1277552 UnitLess; 

# Delta_E= 24.966587 kj; Delta_S= -0.16190599 kJ/K; 


Prob.4.26. Determine the amount of heat which should be supplied to 2 kg of water at 25 C to convert 
it in to steam at 5 bar and 0.9 dry. [VTU-BTD-Dec. 2007-Jan.2008] 
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TEST Solution: 

Following are the steps: 

1 . Select System analysis - Generic - Uniform System from the Daemon tree: 


l 

I I I 

Basic^Touls System Analysis States & Properties 

Unit DeskCai 
Converter 
I 

Unsteady Process Steady State Steady State Unsteady Process 

I , (cycles) I 1 


1 

Tables & 
Cheats 


Closed 


Open 


1 

Uniform 

System 


1 

Uniform 

Flow 

1 


Generic 


: 


Specific 


Generic 


Specific 


T 


; 


Reciprocating Cycles HVAC/Psychrometry Comb. & Equilibrium 


System 

1 


Single-Flow 


I 


Vapor and Gas 
Power Cycles 


I I 1 

Non-Mixing Semi-Mixing Mixing 

Non-Uniform Non-Uniform Non-Uniform 

I I 

Non-Mixing Multi-Flow Mixing Multi-Flow ^ 

I I I f 

Refrigeration HVAC Gas Combustion 

Cycles Psychrometry Dynamics & Equilibrium 


2. Select Phase Change (PC) for Material model, since we are dealing with Water: 


Select a material model to launch the daemon. 




Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (1-120* as opposed to H20), which uses 

rnmnrpggpH lini lirl fahle> for haftar am irarv 

I A> Launches the PC Model 

Generic, Uniform System, Closed Process Daemon _ 12 . N H3. R-134a. N2. C02. etc., should be treated as PC fluids if there is anv 

PC Model 

possibility of a phase transition. 

Examples: R-134a vapor is compressed in a piston-cylinder device from a beginning-state to a final-state. To 
find the work transfer if compression is assumed isentropic. For specific examples, click on the help icon at the 
bottom margin of the daemon. 


SL Model 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 

Examples: A block of copper is heated from a beginning-state to a final-state. To find the heat transfer 
necessary for the process. For specific examples, click on the help icon at the bottom margin of the daemon. 
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3. Enter parameters pi, Tl, ml for State 1; click on Calculate. We get: 

Generic, Uniform-System, Closed Process Daemon: PC Model 


* H - ” □ □ E3 □ [ 


.ijj.y.iMn.u.1,, 




Move mouse over a variable to d 


s value with more pre 


• Mixed c SI English r Case-0 v > | I* Help Messages On Super-lterate | 

State Panel 


SuperCalciiate 


Super-1 nitialize 



| < 1 ©State-1 v 

> 1 

■ 

No-PlotS v 

Initialize 

|subcooled Liquid 

H20 

■ 

pi 


W Tl 


r xl 


r yl 


r vi 


|oo.o 


v 125.0 

| deg-C 

-r 

fraction 

Htz 
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V [0.001 

n^S/kg v | 

r ui 


r hi 


r si 
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1 104.87847 

kJ/kg 

V 1 1 04.97879 

kJ/kg 

V [0.36732 

kJ/kg. K 

v 0.0 

1 m ^ s 


■ v l 

r el 


r ji 


phil 


psil 


i* ml 


1 104.87847 

kJ/kg 

v |l 04.97879 

U/kg 

Z3j 

kJ/kg 


kJ/kg 

v i|^| 

l ( ' 

1 

L 


mi 








]0. 00201 

m*3 


kg/kmol 

E 







4. Similarly for State 2: enter p2, x2, m2 = ml, and click on Calculate. We get: 


Move mouse over a variable to display its value with more precision. 



5. Go o Process Panel. Enter b-state and f-state, and W_B = 0 and W_0 = 0; 
click on Calculate. We get: 


Move mouse over a variable to display its value with more precision 


(* Mixed r SI r English < [©Case-0 v > I 15* Help Messages On 


Process Panel 


| < [process-A[1-2] 

| State-1 v |.p 

1 1 State-2 7) [7 

Initialize | 

[User Defined 

O 

L 

1* WJB 

1? IV_0 

1* T_B 


[4528.298 kJ 

v |0.0 

v ~^^|kJ 

V [298.15 

■|K 


r S_gen 


Delta_E 


Delta_S 



[-3.27107 

kJ/K 

v 14528.298 

1" 

v |l 1 .91692 

kJ/K 



Note that: Delta_E = 4528.298 kj .... Ans. 
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6. Draw the T-s diagram. Constant pressure lines are also shown (in blue): 

T, K 



- 1.69 


s, kJ/kg.K 


11.71 
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7. Click on Super Calculate. TEST code is produced, and see it in I/O panel: 


# Daemon Path: Systems>Closed>Process>Generic>Uniform>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: H20; 

Given: { pl= 100.0 kPa; Tl= 25.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; ml= 2.0 kg; } 
State-2: H20; 

Given: { p2= 500.0 kPa; x2= 0.9 fraction; Vel2= 0.0 m/s; z2= 0.0 m; m2= “ml” kg; } 

} 

Analysis { 

Process- A: b-State = State- 1; f- State = State-2; 

Given: { W_B= 0.0 kj; W_0= 0.0 kj; T_B= 298.15 K; } 

} 

# p; n d of TEST-code 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

100.0 

298.2 

0.001 

104.88 

104.98 

0.367 

#02 

500.0 

425.0 0.9 

0.3385 

2369.03 

2537.75 

6.326 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: b-State = State- 1; f- State = State-2; 

# Given: W_B= 0.0 kj; W_0= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 4528.298 kj; S_gen= -3.2710671 kJ/K; Delta_E= 4528.298 kj; 
Delta_S= 11.916918 kJ/K; 

# Verify: 

ml*(u2-ul) = 4528.298 kj. 
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I Law of Thermodynamics applied to Flow Processes 


5 I Law of Thermodynamics 
applied to Flow Processes 

Learning objectives: 

1. In this chapter, we consider ‘Steady Flow Energy Equation (SFEE)’ and conservation of 
mass for a control volume. 

2. These two principles, viz. Conservation of mass (i.e. continuity equation) and the 
conservation of energy (i.e. First Law) are applied to a number of practically important 
devices such as Nozzles and Diffusers, Turbines and Compressors, Throttling devices, Heat 
Exchangers and Mixing chambers etc. 

3. Transient processes such as filling a tank with a fluid or discharging from a tank are also 
considered. These are known as Uniform State, Uniform Flow (USUF) processes. 

5.1 Formulas used: 

5.1 .1 Steady Flow Energy Equation (SFEE) for a control volume: 

For unit mass flow, i.e. m = 1 kg/s: 

Let: 1 - inlets. 2 - exits 
h = enthalpy kJ/kg V= velocity rn/s 

z= heieht_above_datum m A= area_of_flow n~i A 2 

q = heat_transfer kJ w = work_tr ansfer kJ 


Heat going in to the system is positive, work done by the system is positive. 


Easier way is to remember: Energy going in = Energy going out, in steady state: 


V, 




qj — h| — — g-Z| — W| + h-j + + gz 2 — eqn. (5.1) 

2 2 

f V 2 2 - Vj 2 

i.e. q| - w'i = jh, - - — - +g-(z 2 - z |) ...eqn. (5.2} 

V 2 J 


i.e. qj - wj = Ah + Jite + Ape ..where all terms are for unit mass flow rate.. . eqn. (5.3) 
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When mass flow rate of stream is mj kg/s: 


f -2 } 


„ 2 ^ 


\ 1 


v 2 

Q + ni|- 

h i + — + 3 z i 

K 2 J 

= w + m i 

ll-. 2Z-. 

\ 2 J 


...if there is one stream 
only.... eqn. (5.3-a) 


Note: If there are more than one stream, add additional terms for each stream to take in to account 
respective enthalpies, K.E. and P.E. 

5.1.2 Mass balance: 

Pf A f v i= P2 a 2 -V2 —kg/s eqn.{5_4) 

A. | ■ V | A'2 ■ V2 

or: — - — - = — — - ....kg/s/ where v = sp. volume .... e q 11 . f S . 5 ) 

v t v 2 

5.1 .3 Examples of Steady flow processes: 

1. Nozzle and Diffuser: 

From eqn. (5.1). we have: 


V, 


q_l + h| — + g-z^ = Wj — h.T + + 2 , zt — eqn. (5.1) 


Here, ql = 0, wl - 0. 


Therefore: If change in P.E. is zero and velocity of approach VI = 0. we get: 


v 2 

hi = ...eqn. [5. 6 J 

1 z --5 


0 r: v 2 = J2 \ h t - h 2 ) m/s . . . exit ve 1 0 c ity 
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2. Turbines and compressors: 

For Turbine, it can be taken as insulated, flow velocities small, and K.E. and P.E. terms neglected: 

Then, SFEE for a turbine becomes: 

hi = Iit + w ....for unit mass flow 
i z 

W 

i.e. w= — = hj - h .for unit mass flow rate. . . eqn. (5. 7) 

m ' 


Similarly, SFEE for an adiabatic pump or compressor: 


hi = Iit - w = In 

1 z z 


i.e. w = I 12 _ h| 


W 

m 

...for unit maas flow rate. ..eqn. (5. 3) 
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3. For Throttling device: 


Here.q = 0. w - Q : and changes in P.E. can he neglected. 


Then SFEE reduces to: 



V 2 

h 2 — ...eqn. (5. 9} 


Often, pipe velocities are small. If we neglect changes in K.E. we get: 
h- = h 2 -- eqn. (5.10) 


4. Heat Exchangers: 

When the two streams do not mix, as in a normal HX or condenser: 

Let: c — cold stream, h — hot stream: 1 - inlets, 2 - exits 

Here. Q = 0 r W = 0. and SFEE becomes: 

m c -h cl + m h h hl = m c h c2 + m h h ll2 ...eqn.(5. 1 1) 

When the two streams mix, as in a de-super-heater or cooling tower, or mixing chambers: 

Here, q = 0, w = 0. 

Mass balance: 

m| + m 2 = ixLj ....eqn. (5. 12) 


SFEE: 


m| h| + m ? h 2 = m- h- neglecting changes in P.E. and K.E eqn. (5. 13) 
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5.1 .4 Uniform State, Uniform Flow (USUF) process: 

ex: filling closed tanks with a gas or liquid, discharge from closed vessels etc.: 
Let: 

Q w = heat entering the control volume 
W w = Work leaving the control volume 

i = inlets to control volume 
e = exits from control volume 
m 1 = initial mass in control volume 

m 2 = final mass in control volume 

(m2-m1) = net mass that enters or leaves the control volume 
h = enthalpy 
V = velocity of fluid 
u = internal energy 
z = height from datum 

1 - initial conditions 

2 - final conditions 


Then, for a time period t, the First Law for USUF process is: 


inlets 


f 2 > 

[h + ~r + s- z i 

A J 


f 2 



r v 2 l 


r v 2 

\~ 

=V 

Z-I 

h e + “T + s ' Z e 

^ J 

- 

m-r 

z 

¥2 

Ui + + 2,-Z-i 

J 

- m r 

v i 

L 2 

- sZl 
J A 




+ W,. 


eq.iI5.14) 


As an example of applying the equation (5.14), consider the following: 
Variable flow process: Filling or emptying a tank: 

1. Filling a tank: 


We assume that changes in P.E. and K.E. are negligible. 

Then, we get: 

Let: 

ml = initial mass in tank 
m2 = final mass in tank 

(m2~m1) = mass that enters the tank from the pipe 

Q = heat transfer = 0 

W = Work transfer = 0 

h_p = enthalpy of fluid in pipe 

V_p = velocity of fluid in pipe 

u = internal energy 

1 - initial conditions 

2 - final conditions 
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Making an energy balance: 


TTLyU- - |mi - ml I 


V 


hp 


- 


= m-vu-i 


.„eqn.{5.15) 


If, initially, the tank is empty, then ml = 0 
Then. 


v 


T“ = u 2 — eqn..{5. 16) if ml = 0 


Also, if pipe velocity (i.e. KE) is negligible, then: 
hp = u-? ....eqn. (5.17)... ifK.E. is negligible 

Note: For an Ideal gas. h = cp. T. u = cv.T 
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2. For emptying the tank: 


( m l - m 2)\ h prime + 


V_prime 


- Q = m 


u - m-yu-i 


1 U 1 


...eqn.{5.18) 


where: 

h_prime = sp. enthalpy of leaving fluid 
V_prime = velocity of leaving fluid 


5.2 Problems solved with EES: 

“Prob.5.1. A Nozzle is a device for increasing the velocity of a steadily flowing stream. At the inlet to a 
certain nozzle, the enthalpy of the fluid is 3000 kj/kg and the velocity is 60 m/s. At the discharge end, 
the enthalpy is 2762 kj/kg. The nozzle is horizontal and there is negligible heat loss from it. (i) Find 
the velocity at the exit of the nozzle (ii) If the inlet area is 0.1 m A 2 and the sp. volume at inlet is 
0.187 m A 3/kg, find the mass flow rate (iii) If the sp.vol. at the exit of nozzle is 0.498 m A 3/kg, find the 
diameter of exit section. [VTU-July 2004]” 



Fig.Prob.5.1 


EES Solution: 

“Data:” 

h_l=3000E03 [J/kg]“... enthalpy at inlet“ 
C_l=60 [m/s]“... inlet velocity” 
h_2=2762E03[J/kg] “...enthalpy at exit“ 

A_l=0.1 [m A 2] ££ ... area at inlet” 
v_l=0.187 [m3/kg] ££ ...sp. volume at inlet” 
v_2=0.498 [m3/kg] ££ ...sp. volume at exit” 
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“Calculations:” 

m_l= A_1 * C_l/ v_l “kg/ s... finds mass flow rate” 

m_l= m_2 “...continuity eqn” 

m_2= A_2 * C_2 / v_2 “...finds finds area at exit” 

A_2= pi * D_2 A 2/4 “ finds diameter at exit” 

Q-W = DELTAh + DELTAKE + DELTAPE “...First Law for Open system” 
Q=0 “. . .by data, for nozzle” 

W=0 “...by data, for nozzle” 

DELTAh = h_2-h_l“[J/kg]“ 

DELTAKE=(C_2 A 2/2)-(C_l A 2/2)“[J/kg] “ 

DELTAPE=0“. . .by data” 


Results: 

Unit Settings: SI K kPa kJ molar deg 


Ai = 0.1 [m 2 ] 

C 2 = 692.5 [m/s] 

APE = 0 [J/kg] 
h 2 = 2.762E+06 [J/kg] 
Q =0 [J/kg] 

W = 0 [J/kg] 


A 2 = 0.02307 [m 2 ] 

Ah = -238000 [J/kg] 
D 2 = 0.1 71 A [m] 
mi = 32.09 [kg/s] 
vi = 0.1 37 [m3/kg] 


Ci = 60 [m/s] 

AKE = 238000 [J/kg] 
h-| = 3.000E+06 [J/kg] 
m 2 = 32.09 [kg/s] 
v 2 = 0.-198 [m3/kg] 


Thus: 

Velocity at exit = C2 = 692.5 m/s . . . Ans. 
Mass flow rate = ml = 32.09 kg/s ... Ans. 
Dia at exit = D2 = 0.1714 m ... Ans. 


“Prob.5.2. 12 kg of air per min. is delivered by a centrifugal air compressor. The inlet and outlet 
conditions of air are: VI = 12 m/s, pi = 1 bar, vl = 0.5 m A 3/kg, and V2 = 90 m/s, p2 = 8 bar, v2 = 0.14 
m A 3/kg. The increase in enthalpy of air passing through the compressor is 150 kj/kg and heat loss to 
surroundings is 700 kj/min. Calculate: (i) Power required to drive the compressor, (ii) ratio of inlet to 
outlet pipe diameters. [VTU-Jan. 2003]” 
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Fig.Prob.5.2 



qaiteye 

Challenge the way we run 


EXPERIENCE THE POWER OF 
FULL ENGAGEMENT... 


RUN FASTER. 

RUN LONGER.. 
RUN EASIER... 




110 




Click on the ad to read more 


Download free eBooks at bookboon.com 


Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


EES Solution: 

“Data:” 

m=12 * convert(kg/min, kg/s) “[kg/s]” 

C_l = 12 [m/s]“... inlet velocity” 
pi = 10 A 5 “Pa ... inlet pressure” 
v_l=0.5“m A 3/kg ... sp. volume at inlet” 

C_2 = 90 [m/s] “...exit velocity” 
p2 = 8E05“Pa .... exit pressure” 
v_2 = 0.14“..m A 3/kg ... sp. vol. at exit” 

Q= - (700E03)/(m * 60)“... J/kg” 

DELTAh=150E03“ J/kg.... change in enthalpy“ 
DELTAKE=(C_2 A 2/2)-(C_l A 2/2)“J/kg.... change in K.E.“ 
DELTAPE=0“. . . change in P.E.” 


“Calculations:” 

Q - W = DELTAh + DELTAKE + DELTAPE “...First 
A_l= m * v_l/C_l“..finds area at inlet, m A 2” 

A_2= m * v_2/C_2“... finds area at exit, m A 2” 
DlbyD2= sqrt(A_l/A_2)“... finds dia ratio” 

W_act = W * m . .finds Work required, W” 

Results: 

Unit Settings: SI K kPa kJ molar deg 


A-| = 0.008333 [m 2 ] 
D1 byD2 = 5.1 75 
m = 0.2 [m/s] 
vi = 0.5 [m3/kg] 


A 2 = 0.00031 1 1 [m 2 ] 
Ah = 1 50000 [J/kg] 
pi = 1 00000 [Pa] 
v 2 = 0.1 A [m3/kg] 


Law for open system’ 


Ci = 12 [m/s] 

AKE = 3078 [J/kg] 
p2 = 800000 [Pa] 
W= -21 2311 [J/kg] 


C 2 =00 [m/s] 
APE = 0 [J/kg] 

Q =-58333 [J/kg] 
W ac t = -42462 [vV] 


Thus: 

Power required for compressor = -42462 W . . . Ans. (negative sign indicates power input to system) 
Ratio of inlet to exit diameters = 5.175 .... Ans. 
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(b) Plot the variation of Power required as the heat loss varies from 0 to -8000 J/kg: 
First, compute the Parametric Table: 


Table 1 

► 

1..9 

ra 

Q 

[J/kg] 

2 IT 

w act 

[W] 

Run 1 

0 

-30796 

Run 2 

-1000 

-30996 

Run 3 

-2000 

-31196 

Run 4 

-3000 

-31396 

Run 5 

4000 

-31596 

Run 6 

-5000 

-31796 

Run 7 

-6000 

-31996 

Run 8 

-7000 

-32196 

Run 9 

-8000 

-32396 


Now, plot the results: 


x 10 3 



Q [J/kg] 


Note that as the heat transfer increases, compressor power required also increases. 
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a Prob.5.3. Air flows steadily through a rotary compressor. At entry, the air is at 20 C and 101 kPa. At the 
exit, the air is at 200 C and 600 kPa. Assuming the flow to be adiabatic, (i) evaluate the work done per 
unit mass of air if the velocities at inlet and exit are negligible (ii) what would be the increase in work 
input if the velocities at inlet and exit are 50 m/s and 110 m/s? [VTU-Jan. 2005]” 

EES Solution: 

“Data:” 

DELTAKE=0 
DELTAPE=0 
Q=0 “...since adiabatic” 

Tl=20 “C“ 

T2=200 “C“ 
pl = 101E03 “Pa“ 
p2=600E03 “Pa“ 

R=287 “J/kg.K .... for air” 
gamma=1.4 “for air” 
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“Calculations:” 

“Case 1: Inlet and exit velocities are negligible:” 

cp = R * gamma/ (gamma - 1) sp. heat at const, pressure, J/kg.K” 

DELTAh = cp * (T2-T1) “J/kg ... change in enthalpy 4 

Q - W = DELTAh + DELTAKE + DELTAPE “...by First Law to Open systems” 

“Case 2: Inlet and exit vel. not negligible:” 

Q1=0 “...since adiabatic” 

Ql- Wl= cp * (T2-T1) + (V2 A 2-Vl A 2)/2 “...First Law for Open system, including the change in K.E.” 
V2=110 “m/s ... exit velocity” 

VI =50 “m/s ... inlet velocity” 


Results: 

Unit Settings: SI K 

cp = 1 005 [J/kg-C] 
APE = 0 [J/kg] 
p2 = 600000 [Pa] 

R = 237 [J/kg-C] 

VI =50 [m/s] 

W1 =-135610 [J/kg] 


kJ molar deg 

Ah = 1 3031 0 [J/kg] 
v = 1 A 

i 

Q =0 [J/kg] 

T1 =20 [C] 

V2 =110 [m/s] 


AKE = 0 [J/kg] 
pi = 1 01 000 [Pa] 
Ql =0 [J/kg] 

T2 = 200 [C] 

W = -1 3031 0 [J/kg] 


Thus: 

W = -180.81 kj/kg when K.E. and RE. are neglected... Ans. 

W1 = -185.61 kj/kg when K.E. is not negligible ...i.e. an increase of about 5 kj/kg... Ans. 
(Note: negative sign indicates work input in to the system.) 


“Prob.5.4. In a gas turbine unit, the gases flow through the turbine at 15 kg/s and the power developed 
by the turbine is 12000 kW. The enthalpies of the gases at the inlet and outlet are 1260 kj/kg and 400 
kj/kg respectively, and the velocities of gases at the inlet and outlet are 50 m/s and 110 m/s respectively. 
Calculate (i) rate at which heat is rejected by the turbine (ii) the area of inlet pipe, given the sp. vol. of 
gases at inlet is 0.45 m A 3/kg. [VTU-Jan. 2005]” 
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Gas, In 



Fig.Prob.5.4 


EES Solution: 

“Data:” 

h_l = 1260E03 “J/kg ... enthalpy at inlet' “ 
h_2=400E03 “J/kg . . . enthalpy at outlet” 

C_l=50 [m/s] “...velocity at inlet” 

C_2=110 [m/s] “...velocity at outlet” 
v_l=0.45 [m A 3/kg]“..sp. vol. at inlet“ 
m_l=15 [kg/s] ”... mass flow rate” 

“Calculations:” 

m_l=A_l * C_l/ v_l “...finds area at inlet, m A 2” 

Q - W = DELTAh + DELTAKE + DELTAPE “...First Law for Open System .... finds Q” 

W = 12E06 [ J/s]/15 [kg/s] “J/kg work output of turbine, by data” 

DELTAh = (h_2-h_l) “J/kg... enthalpy change” 

DELTAKE=(C_2 A 2/2)-(C_l A 2/2)“J/kg ... change in K.E.” 

DELTAPE=0“... change in P.E.” 
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Results: 


Unit Settings: SI K kPa kJ molar deg 


Ai-0.135 [m 2 ] 


Ci = 50 [m/s] 


C 2 =110 [m/s] 
APE = 0 [J/kg] 
mi =15 [kg/s] 


Ah =-860000 [J/kg] 
hi = 1 .260E+06 [J/kg] 


AKE = 4800 [J/kg] 
h 2 = 400000 [J/kg] 
vi = 0.45 [m 3 /kg] 


Q = -55200 [J/kg] 


W = 800000 [J/kg] 


Thus: 

Heat rejected by turbine = Q = 55200 J/kg... .negative sign indicates heat going out of the 
system.... Ans. 

Area of inlet pipe = A1 = 0.135 m A 2 ... Ans. 


“Prob. 5.5. In a steady flow system, 50 kJ of work is done per kg of fluid; values of sp. vol., pressure and 
velocity at inlet and exit sections are: 0.4 m A 3/kg, 600 kPa, 15 m/s and 0.6 m A 3/kg, 100 kPa, and 250 
m/s, respectively. The inlet is 30 m above the exit. The heat loss from the system is 8 kj/kg. Calculate 
the change in internal energy per kg of fluid. [VTU-July 2003]” 

EES Solution: 


Q = - 8E03 “J/kg ... heat rej.” 

W=50E03 “J/kg .... work done by fluid” 
pl=600E03 “Pa ...inlet pressure” 
vl=0.4 “m3/kg... inlet sp. vol.” 

Cl =15 “m/s .. inlet velocity” 
p2=100E03 “Pa ... exit pressure” 
v2=0.6 “m3/kg ... exit sp. vol.” 

C2=250 “m/s ... exit velocity” 

Z2=0 “m ... exit datum level” 

Z 1 -JO “m ... inlet datum level” 

“Calculations:” 

Q - W = DELTAh + DELTAKE + DELTAPE “...First Law for Open System .... finds DELTAh” 
DELTAKE=(C2 A 2-Cl A 2)/2 “J/kg ..., change in K.E.” 


“Data:” 
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DELTAPE=g * (Z2 - Zl) “J/kg ..., change in P.E.” 
g=9.81 “m/s2 ... accn. due to gravity” 

DELTAh = DELTAu + (p2 * v2 - pi * vl) “...from h = u + pV finds DELTAu” 


Results: 

Unit Settings: SI K 

Cl =15 [m/s] 

AKE =31138 [J/kg] 
g = 9.81 [m/s 2 ] 

Q = -8000 [J/kg] 

W = 50000 [J/kg] 


kJ molar deg 

C2 = 250 [m/s] 
APE = -29*1.3 [J/kg] 
pi = 600000 [Pa] 
vl = 0 A [m 3 /kg] 

Zl = 30 [m] 


Ah = -88843 [J/kg] 
Au =91157 [J/kg] 
p2 = 100000 [Pa] 
v2 = 0.6 [m 3 /kg] 

Z2 = 0 [m] 


Thus: 

Change in Internal energy = DELTAu = 91157 J/kg .... Ans. 
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(b) Plot the variation of DELTAu as heat rejected Q varies from 0 to -10 kj/kg: 
First, calculate the Parametric Table: 


0 

i iz 

Q 

2 IZ 

AU 

1..11 

[J/krj] 

[J/kg] 

Run 1 

0 

99157J 

Run 2 

-1000 

98157 

Run 3 

-2000 

97157 

Run 4 

-3000 

96157 

Run 5 

4000 

95157 

Run 6 

-5000 

94157 

Run 7 

-6000 

93157 

Run 3 

-7000 

92157 

Run 9 

-8000 

91157 

Run 10 

-9000 

90157 

Run 11 

-10000 

89157 


Now, plot the results: 


x 10 3 



Note: Negative sign for Q only indicates that heat is being rejected. 


118 


Download free eBooks at bookboon.com 


Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


a Prob.5.6. Air flows steadily at a rate of 0.5 kg/s through a compressor, entering at 7 m/s velocity, 100 
kPa pressure and 0.95 m A 3/kg sp. volume, and leaves at 700 kPa, 5 m/s, and 0.19 m A 3/kg. The internal 
energy of air leaving is 90 kj/kg greater than that of air entering. Cooling water in the compressor 
jacket absorbs heat at a rate of 58 kW. Compute the shaft work input and the ratio of inlet to exit pipe 
diameters. [VTU- July 2002]” 


© © 

+ — 



Air, In 


EES Solution: 

“Data:” 

m=0.5 “kg/s .... mass flow rate” 

Q= - 58E03 cc J/s ...heat rejection rate” 
pl = 100E03 cc Pa ... inlet pressure” 
p2=700E03 “Pa .... exit pressure” 

0=7.0 “m/s ... inlet velocity” 

C2=5.0 “m/s .... exit velocity” 

DELTAu=90E03 “J/kg .... change in internal energy” 
vl=0.95 “m A 3/kg ... sp. vol. at inlet” 
v2=0.19 “m A 3/kg .... sp. vol. at exit” 


“Calculations:” 

Q - W = m * (DELTAu+(p2 * v2 - pi * vl)) + m * (C2 A 2-Cl A 2)/2 “...by First Law for Open system... 
finds W” 

Al=m * vl/Cl “m A 2 ... inlet pipe area” 

A2=m * v2/C2 “m A 2 ... exit pipe area” 

DlbyD2=sqrt(Al/A2) “...diameter ratio” 
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Results: 


Unit Settings: SI K kPa kJ molar deg 

A2 = 0.01 9 [m 2 ] 
D1byD2 = 1.89 


A1 = 0.06786 [m 2 ] 
C2 =5 [m/s] 
m = 0.5 [kg/s] 

Q =-53000 [W| 

W = -1 21 OEM [v V] 


pi = 100000 [Pa] 
vl = 0.95 [m 3 /kg] 


Cl =7 [m/s] 

Au = 90000 [J/kg] 
p2 = 700000 [Pa] 
v2 = 0.1 9 [m 3 /kg] 


Thus: 

Work input to compressor, W = -121.994 kW .... Ans. 
Diameter ratio = D1/D2 = 1.89 ... Ans. 
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a Prob.5.7. A centrifugal air compressor compresses 5.7 m A 3/min of air from 85 kPa to 650 kPa. The 
initial sp. vol. is 0.35 m A 3/kg and final sp. vol. is 0.1 m A 3/kg. If the suction line dia is 10 cm and the 
discharge line dia is 6.25 cm, determine: (i) the change in flow work (ii) the mass rate of flow, and (iii) the 
velocity change. [VTU-Aug. 2000]” 



Air 

out 


Fig.Prob.5.7 


EES Solution: 

“Data:” 

Vl=5.7/60“m3/s .... volume flow rate” 
pl=85E03“Pa ... inlet pressure” 
p2=650E03 “Pa .... exit pressure” 
rhol = l/0.35“kg/m3 ... inlet density” 
rho2=l/0.1“kg/m3 ... exit density” 
dl=0.1“m .. inlet dia” 
d2=0.0625“m .. exit dia” 

“Calculations:” 

m = VI * rhol“kg/s .... mass flow rate” 
Al=pi*dl A 2/4 “m2 ... inlet area” 
A2=pi*d2 A 2/4 “m2 .... exit area” 
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“Change in Flow work:” 

DELTAPV= m * (p2/rho2 - pl/rhol) cc J/s .... change in flow work” 


“Velocity change:” 

m = Al*Cl* rhol “... finds inlet velocity, Cl, m/s 3 
m = A2 * C2 * rho2 “. . . finds exit velocity, C2, m/s” 

Results: 

Unit Settings: SI K kPa kJ molar deg 

A1 = 0.00785^ [m 2 ] A2 = 0.003068 [m 2 ] 
dl = 0.1 [m] d2 = 0.0625 [m] 

pi = 85000 [Pa] p2 = 650000 [Pa] 

VI = 0.095 [m3/s] 


Cl =12.1 [m/s] 


C2 = 8.847 [m/s] 

.iPV = 9568 [J/s; 

1 

m = 0.271 4 [kg/s] 


rhol = 2.857 [kg/m 3 ] rho2 =10 [kg/m 3 ] 


Thus: 

Change in flow work = DELTApv = 9568 W . . . Ans. 
Mass flow rate = m = 0.2714 kg/s ... Ans. 

Velocities: Cl = 12.1 m/s, C2 = 8.847 m/s .... Ans. 


“Prob.5.8. A steam turbine receives steam with a flow rate of 900 kg/min. and experiences a heat loss of 
840 kj/min. The exit pipe is 3 m below the level of inlet pipe. Find the power developed by the turbine if 
the pressure decreases from 62 bar to 9.8 kPa, velocity increases from 30.5 m/s to 274.3 m/s, internal energy 
decreases by 938.5 kj/kg, and sp. vol. increases from 0.058 m A 3/kg to 13.36 m A 3/kg. [VTU-Feb. 2002]” 
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EES Solution: 

“Data:” 

m = 900/60 “kg/s” 

Q = - 840E03/60”J/s” 

Zl= 3 “m” 

Z2 = 0 “m“ 
pl= 62E05 “Pa“ 
p2 = 9.86E03 “Pa“ 

Cl= 30.5 “m/s” 

C2 = 274.3 “m/s” 

DELTAu = - 938.5E03 “J/kg” 
vl= 0.058 “m3/kg” 
v2= 13.36 “m3/kg” 
g=9.81“m/s2” 


“Calculations:” 

Q - W = m * (DELTAu + (p2 * v2 - pi * vl)) + m * (C2 A 2-Cl A 2)/2 + m * g * (Z2-Z1) “...First Law 
for Open system” 


Results: 


Unit Settings: SI K kPa kJ molar deg 


Cl = 30.5 [m/s] 
m = 1 5 [kg/s] 
vl = 0.058 [m 3 /kg] 
Z2 = 0 [m] 


C2 = 27-1.3 [m/s] 
pi = 6.200E+0G [Pa] 
v2 = 1 3.36 [m 3 /kg] 


Au =-938500 |W] g= 9.81 [m/s 2 ] 

Q = -14000 [W] 
Z1 = 3 [m] 


p2 = 9860 [Pa] 

W = 1 .692E+07 [vV] 


Thus: 

Power developed by turbine = W = 1.692E07 W .... Ans. 


“Prob. 5.9. A fluid flows through a steady flow system at the rate of 3 kg/s. The inlet and outlet conditions 
are: pi = 5 bar, Cl = 150 m/s, ul = 2000 kj/kg, and p2 = 1.2 bar, C2 = 80 m/s, and u2 = 1300 kj/kg. The 
change in sp. vol. is from 0.4 m A 3/kg to 1.1 m A 3/kg. The fluid loses 25 kj/kg heat during the process. 
Neglecting potential energy, determine the power output of the system. [VTU-Dec. 2006-Jan. 2007]:” 
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EES Solution: 

“Data:” 

m = 3 “kg/s” 

PI = 500 “kPa” 

Cl = 150 “m/s” 
ul = 2000 “kj/kg” 

P2 = 120 “kPa” 

C2 = 80 “m/s” 

u2 = 1300 “kj/kg” 

vl = 0.4 “m3 /kg” 

v2 = 1.1 “m3 /kg” 

q = -25 “kj/kg.... heat loss” 

“Neglecting Potential energy, determine the Power output:” 


“Write SFEE for 1 kg:” 

q - w = (h2 - hi) + ((C2 A 2 - Cl A 2)/2)*10 A (-3) “...First Law for Open system. ..all quantities in kj/kg” 
hi = ul + PI * vl“kj/kg ... inlet enthalpy” 
h2 = u2 + P2 * v2 “kj/kg . . . exit enthalpy” 

“Power output:” 

Work = w * m “kj/s” 


Results: 

Unit Settings: SI 

Cl =150 [m/s] 
m = 3 [kg/s] 
ul = 2000 [kJ/kg] 
w =751.1 [kJ/kg] 


kPa kJ mass deg 

C2 =80 [m/s] 

PI = 500 [kPa] 
u2 = 1 300 [kJ/kg] 
Work = 2253 [kvV] 


hi = 2200 [kJ/kg] 
P2 =120 [kPa] 
vl =0.4 [m^/kg] 


h2 = 1 432 [kJ/kg] 
q = -25 [kJ/kg] 
v2 = 1.1 [m 3 /kg] 


Thus: Work done by the system - Work - 2253 kW . . . Ans. 


“Prob.5. 10. A fluid flows through a steam turbine at a steady rate of 5000 kg/h, while energy is transferred 
as heat at a rate of 6279 kj/h from the turbine. The condition of the fluid at the turbine inlet and exit 
are: hi - 3153 kj/kg, Cl - 60 m/s, Z1 = 6 m, and h2 = 2713 kj/kg, C2 - 185 m/s, Z2 = 4 m. Find the 
power output from the turbine. Comment on K.E. and PE. changes. [VTU-Dec. 08-Jan. 09]” 
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EES Solution: 

“Data:” 

mass_flow = 5000/3600 “kg/s” 

q = (-6279/3600) / mass_flow “kj/kg....heat transf. from turbine” 
hi = 3153 “kj/kg” 
h2 = 2713 “kj/kg” 

Cl = 60 “m/s” 

C2 = 185 “m/s” 

Z1 = 6 “m” 

Z2 = 4 “m” 
g = 9.81“m/s2” 

“Find the Power output from Turbine and comment on K.E. and P.E. changes:” 

q - w = DELTAh + DELTAke + DELTApe “..First Law for a turbine.... all terms are in kj/kg” 
DELTAh = (h2 - hi) “kj/kg” 

DELTAke = ((C2 A 2 - Cl A 2)/2) * 10 A (-3) “kj/kg” 

DELTApe = g * (Z2 - Zl) * 10 A (-3) “kj/kg” 

Work = mass_flow * w “kj/s” 

Results: 

Unit Settings: SI C kPa kJ mass deg 

Cl =60 [m/s] 


massfiow = 1.330 [kg/s] 
w = 423.5 [kJ/kg] 

Zl = 6 [m] 

Thus: 

Work done by turbine = 588.1 kW .... Ans. 

DELTAke = 15.31 kj/kg, DELTApe = -0.01962 kj/kg .. both are negligible compared to the enthalpy 
difference DELTAh = -440 kj/kg .... Ans. 


Ah = -440 [kJ/kg] 

Ape =-0.01962 [kJ/kg] 
hi = 31 53 [kJ/kg] 


C2 =165 [m/s] 

Ake =15.31 [kJ/kg] 
g= 9.81 [m/s 2 ] 
h2 = 271 3 [kJ/kg] 
q = -1.256 [kJ/kg] 
Work = 588.1 [kW] 
Z2 = 4 [m] 
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(b) Plot the variation of Work as heat loss q varies from 0 to -4 kj/kg: 
First, compute the Parametric Table: 



Now, plot the results: 



Note: Negative sign for q only indicates that heat is being rejected from turbine. As the heat rejected 
increases, work output from the turbine decreases. 
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“Prob.5.11. The working fluid in a steady flow process flows at a rate of 220 kg/min. The fluid 
rejects 100 kj/s of heat passing through the system. The conditions of fluid at inlet and outlet are: 
Cl = 220 m/s, pi = 6 bar, ul = 2000 kj/kg, vl = 0.36 m A 3/kg, and C2 = 140 m/s, p2 = 1.2 bar, 
u2 = 1400 kj/kg, v2 = 1.3 m A 3/kg. Suffix 1 indicates inlet, and 2 the outlet. Determine the power capacity 
of the system in MW. [VTU-BTD- June- July 2009:]” 

EES Solution: 

“Data:” 

mass_flow = 220/60 “kg/s” 
q = -100/mass_flow “kj/kg” 

Cl = 220 “m/s” 

PI = 600 “kPa” 
ul = 2000 “kj/kg” 
vl = 0.36 “m3/kg” 

C2 = 140 “m/s” 

P2 = 120 “kPa” 
u2 = 1400 “kj/kg” 
v2 = 1.3 “m3/kg” 
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“Determine the Power capacity of the system in MW:” 

“Apply I Law to Open System: Energy going In = Energy going Out:” 

(ul+ PI * vl) + (Cl A 2/2)/1000 + q = (u2 + P2 * v2) + (C2 A 2/2)/1000 + w “...where w is work done in 
kj/kg” 

Work = w * mass_flow /1000 “MW” 


Results: 

Unit Settings: SI C kPa kJ mass deg 


Cl = 220 [m/s] 
PI = 600 [kPa] 
ul = 2000 [kJ/kg] 
v2 = 1 .3 [m 3 /kg] 


C2 =140 [m/s] 
P2 =120 [kPa] 
u2 = 1 400 [kJ/kg] 
w = 647.1 [kJ/kg] 


massif = 3.667 [kg/s] 
q = -27.27 [kJ/kg] 
vl = 0.36 [m 3 /kg] 

Work = 2.373 [MW] 


Thus: 

Work done by turbine = 2.373 MW . . . Ans. 


(b) Plot the variation of Work as heat rejected q varies from 0 to -50 kj/kg: 
First, compute the Parametric Table: 




Table 1 | 


l> 

1..11 

IZ 

q 

[kJ/kg] 

2 IZ 

Work 

[MW] 

Run 1 

0 

2473 

Run 2 

-5 

2.454 

Run 3 

-10 

2436 

Run 4 

-15 

2.418 

Run 5 

-20 

2.399 

Run 6 

-25 

2.331 

Run 7 

-30 

2.363 

Run 8 

-35 

2.344 

Run 9 

40 

2.326 

Run 10 

45 

2.308 

Run 11 

-50 

2.239 
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Now, plot the results: 



Note: Negative sign for q only indicates that heat is being rejected from turbine. As the heat rejected 
increases, work output from the turbine decreases. 


“Prob.5.12. A turbine operating under steady flow conditions receives steam at the following state: 
Pressure = 13.8 bar, sp. vol. = 0.143 m A 3/kg, sp. int. energy = 2590 kj/kg, Velocity = 30 m/s. The state 
of steam leaving the turbine is: pressure = 0.35 bar, sp. vol. = 4.37 m A 3/kg, sp. int. energy = 2360 kj/kg, 
velocity = 90 m / s. Heat is rejected to surroundings at the rate of 0.25 kW and the rate of steam flow through 
the turbine is 0.38 kg/s. Calculate the power developed by the turbine. [VTU-BTD-June-July 2008:]” 
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EES Solution: 

“Data:” 

mass_flow = 0.38 “kg/s” 
q = - 0.25 [kj/s] / mass_flow “kj/kg” 

Cl = 30 “m/s” 

PI = 1380 “kPa” 
ul = 2590 “kj/kg” 
vl = 0.143 “m3/kg” 

C2 = 90 “m/s” 

P2 = 35 “kPa” 
u2 = 2360 “kj/kg” 
v2 = 4.37 “m3/kg” 

“Determine the Power capacity of the Turbine:” 

“Apply I Law to Open System: Energy going In = Energy going Out:” 

(ul+Pl * vl) + (Cl A 2/2)/1000 + q = (u2 + P2 * v2) + (C2 A 2/2)/1000 + w “..finds w, where w is work 
done in kj/kg” 

Work = w * mass_flow “kW” 


Results: 

Unit Settings: SI C kPa kJ mass deg 


Cl =30 [m/s] 

PI = 1 330 [kPa] 
ul = 2590 [kJ/kg] 
v2 = <1.37 [m 3 /kg] 


C2 =90 [m/s] 

P2 =35 [kPa] 
u2 = 2360 [kJ/kg] 
w = 270.1 [kJ/kg] 


massf| ow = 0.38 [kg/s] 
q = -0.6579 [kJ/kg] 
vl = 0.1 A3 [m 3 /kg] 
Work = 1 02.7 [kW] 


Thus: Power developed by the turbine = 102.7 kW ... Ans. 


“Prob.5.13. Air enters an adiabatic horizontal nozzle at 400 C with a velocity of 50 m/s. The inlet area is 
240 cm A 2. Temp of air at exit is 80 C. Given that the sp. vol. of air at the inlet and exit are respectively 0.2 
m A 3/kg and 1.02 m A 3/kg, find the area of cross-section of the nozzle at the exit. Assume that enthalpy 
of air is a function of temp only and that cp = 1.005 kJ/kg.K. [VTU-BTD-July 2006:]” 
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EES Solution: 

“Data:” 

T1 = 400+273 “K” 

Cl = 50 “m/s ... velocity at inlet” 

A1 = 240*10 A (-4)“m2 ... area at inlet” 

T2 = 80+273 “K” 
vl = 0.2 cc m3/kg” 
v2 = 1.02 cc m3/kg” 
cp = 1.005 “kJ/kg.K” 
q = 0 “...since adiabatic” 

w = 0 “. . .since there is no work output in nozzle” 


“Calculations: 4 * 

DELTAh = cp * (T2 - Tl) “kj/kg... change in enthalpy” 

q - w = DELTAh + ((C2 A 2 - Cl A 2)/2) * 10 A (-3) “... First Law for Nozzle ...Finds Velocity at exit” 
A1 * Cl/ vl = A2 * C2 / v2 “Finds A2, area at exit” 


Results: 

Unit Settings: SI C kPa kJ mass deg 


Tl = 673 [K] 
vl = 0.2 [m 3 /kg] v2 = 1 .02 [m 3 /kg] 


A1 = 0.02-1 [m 2 ] 
C2 = 803.6 [m/s] 
q = 0 [kJ/kg] 


A2 = 0.007016 [m 2 ] 
cp = 1 .005 [kJ/kg-K] 


Cl =50 [m/s] 

Ah =-321.6 [kJ/kg] 
T2 = 353 [K] 
w =0 [kJ/kg] 


Thus: Area of cross-section of nozzle at exit = A2= 76.16 cm A 2 ... Ans. 
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a Prob.5.14. At the inlet to a certain nozzle, the enthalpy of the fluid is 3025 kj/kg and the velocity is 60 
m/s. At the exit of the nozzle, the enthalpy is 2790 kj/kg. The nozzle is horizontal and there is a heat loss 
of 100 kj/kg from it. Calculate the velocity of fluid at nozzle exit. Also find the mass flow rate of fluid if 
inlet area is 0.1 m A 2 and sp. vol. at inlet is 0.19 m A 3/kg. [VTU-BTD-July-2007:]” 



Fig.Prob.5.14 


EES Solution: 

“Data:” 

hi = 3025 “kj/kg” 

Cl = 60 “m/s” 
h2 = 2790 “kj/kg” 

q = -100 “kj/kg.... heat loss from Nozzle” 
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“Calculations:” 

A1 = 0.1 “m2” 

vl = 0.19 cc m3/kg....sp. vol. at inlet” 

q - w = (h2 - hi) + ((C2 A 2 - Cl A 2)/2) * 10 A (-3) “First Law for nozzle, neglecting PE .... finds C2” 
w = 0 “No work done in Nozzle” 

“Mass flow rate:” 

mass_flow = Al*Cl/vl “kg/s” 

Results: 

Unit Settings: SI C kPa 

A1 =0.1 [m 2 ] 
h£ = 2790 [kJ/kg] 
w = 0 [kJ/kg] 


mass deg 

Cl =60 [m/s] 


C2 = 523.1 [m/s] 


massf| ow =31.58 [kg/s] q = -1 00 [kJ/kg] 


hi = 3025 [kJ/kg] 
vl =0.19 [m 3 /kg] 


Thus: 

Exit velocity, C2 = 523.1 m/s ... Ans. 
Mass flow rate = 31.58 kg/s ... Ans. 


“Prob.5.15. An air receiver of volume 6 m A 3 contains air at 15 bar and 40.5 C. A valve is opened and 
some air is allowed to blow out to atmosphere. The pressure of air in the receiver drops rapidly to 12 
bar and then the valve is closed. Calculate the mass of air blown out. [Ref. 4]” 

EES Solution: 

“Data:” 

Voll = 6[m A 3] 

PI = 15E05[Pa] 

T1 = 40.5 + 273 [K] 

P2 = 12E05[Pa] 

R_air = 287[J/kg-K]“...Gas const, for air” 
gamma =1.4 “...ratio of sp. heats for air” 
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“Calculations:” 

T2/T1 = (P2/P1) A ((gamma - 1) /gamma) “..for isentropic expn.... finds T2 (K)” 
m_l = (PI * Voll) / (R_air * Tl) “kg ... initial mass of air in the receiver” 
m_2 = (P2 * Voll) / (R_air * T2) “kg ... final mass of air in the receiver” 
mass_blown = (m_l - m_2) “kg ... mass blown out” 


Results: 

Unit Settings: SI C kPa kJ 

y = 1 A 

r 

m 2 = 85.29 [kg] 

R air = 287 [J/kg-K] 

Voll = 6 [m 3 ] 


mass deg 

massbbwn = 1 4-74 [kg] 
PI = 1 .500E+06 [Pa] 
Tl = 31 3.5 [K] 


mi =100 [kg] 

P2 = 1 .200E+06 [Pa] 
T2 =294.1 [K] 


Thus: 

Final temp of air in the receiver = T2 = 294.1 K ... Ans. 
Mass of air blown out = 14.74 kg . . . Ans. 
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(b) Plot the final temp and mass blown out as the P2 varies from 6 bar to 12 bar: 
First, compute the Parametric Table: 



Now, plot the results: 
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“Prob.5.16. A tank has a volume of 0.4 m A 3 and is evacuated. Steam at a pressure of 1.4 MPa, 300 C is 
flowing in a pipe and is connected to this tank. The valve is opened and the tank is filled with steam until 
the pressure is 1.4 MPa, and then the valve is closed. If the process takes place adiabatically and K.E. 
and P.E. are negligible, determine the final temp of steam in the tank, and the amount filled in. [Ref.2]” 

EES Solution: 

“Data:” 

Voll = 0.4 [m A 3] 

ml = 0 “...initial mass, since the tank is evacuated” 

P_pipe = 1400[kPa] 

T_pipe = 300 [C] 

h_pipe = Enthalpy(Steam_NBS,T=T_pipe,P=P_pipe)“kJ/kg .... note the use of built-in Function for 
enthalpy of Steam in EES” 

“Let the final temp of fluid after filling in in the tank be T2 deg.C” 

“Then: h_pipe = u2 by First Law;” 

u2 = IntEnergy(Steam_NBS,T=T2,P=P_pipe) “kj/kg ... since tank is filled to a pressure of the steam 
in the pipe” 
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“Note the use of built-in Function for Int. energy of Steam in EES” 
h_pipe = u2 “...finds T2” 

v2 = Volume(Steam_NBS,T=T2,P=P_pipe)“m A 3/kg .... sp. vol. of steam in tank” 
mass = Voll/v2 “kg ... mass filled in the tank” 

Results: 

Unit Settings: SI C kPa kJ mass deg 

Ppipe-HOO [kPa] 
v2 = 0.2357 

Voll = 0 A [m 3 ] 

Thus: 

Final temp of steam in tank = T2 = 452 deg.C.... Ans. 

Amount of steam filled in the tank = 1.697 kg ... Ans. 


hpipe = 3040 [kJ/kg] 


T2 =452 [C] 


ml =0 [kg] 
Tpjpe = 300 [C] 


mass = 1 .697 [kg] 


u2 = 3040 [kJ/kg] 
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a Prob.5.17. Let the tank in the previous example contain initially sat. vapour at 350 kPa. Now the valve 
is opened and the tank is filled with steam until the pressure is 1.4 MPa, and then the valve is closed. If 
the process takes place adiabatically and K.E. and P.E. are negligible, determine the final temp of steam 
in the tank, and the amount filled in. [Ref.2]” 

EES Solution: 

“Data:” 

Voll = 0.4 [m A 3] “..volume of tank” 

PI = 350[kPa] “...Initial pressure of steam in tank” 
xl = 1 “...sat. vapour” 

P_pipe = 1400[kPa] “...pressure of steam in pipe” 

T_pipe = 300 [C] “...temp of steam in pipe” 

“Calculations:” 

ul = IntEnergy(Steam_NBS,x=xl,P=Pl) “kj/kg ...Int. energy of steam in the beginning” 
vl = Volume(Steam_NBS,x=xl,P=Pl) “m A 3/kg .... sp. vol. of steam present initially in tank” 
ml = Voll / vl “...initial mass of steam in the tank “ 

h_pipe = Enthalpy(Steam_NBS,T=T_pipe,P=P_pipe) “kj/kg ....enthalpy of steam in the pipe” 

“Let the final mass in tank be m2, temp of fluid after filling in the tank be T2 deg.C” 

“Then: (m2-ml)* h_pipe = (m2 * u2 - ml * ul) ... by First Law;” 

u2 = IntEnergy(Steam_NBS,T=T2,P=P_pipe) “kj/kg ..int. energy of steam after the tank is filled to a 
pressure of the steam in the pipe” 

v2 = Volume(Steam_NBS,T=T2,P=P_pipe)“m A 3/kg .... sp. vol. of steam in tank, after filling up” 
m2 = Voll / v2 “kg.. mass of steam in tank after filling” 

(m2 - ml) * h_pipe = (m2 * u2 - ml * ul) “...By First Law for filling process” 

“Mass of steam flowing in to the tank:” 

mass_to_tank = (m2 - ml) “kg” 
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Results: 

Unit Settings: SI C kPa kJ mass deg 


h pipe = 3CM0 [kJ/kg] 
mass to,tank = 1 -264 [kg] 


T2 = 3*11 .6 [C] 
u2 = 2855 [kJ/kg] 
Voll = 0.4 [m 3 ] 

Thus: 


ml = 0.763 [kg] 

PI = 350 [kPa] 
Tpipe = 300 [C] 
vl = 0.52-13 [m 3 /kg] 
xl =1 


m2 = 2.027 [kg] 
Ppipe-HOO [kPa] 
ul = 25-19 [kJ/kg] 
v2 = 0.1 973 [m 3 /kg] 


Final temp of steam in tank = 341.8 deg. C ... Ans. 
Mass of steam entering the tank = 1.264 kg ... Ans. 


Prob.5.18. A balloon initially contains 65 m A 3 of helium gas at atmospheric conditions of 100 kPa and 
22 C. The balloon is connected by a valve to a large reservoir that supplies helium gas at 150 kPa and 25 
C. Now the valve is opened, and helium is allowed to enter the balloon until pressure equilibrium with 
the supply line is reached. The material of the balloon is such that its volume increases linearly with 
pressure. If no heat transfer takes place during this process, determine the final temp in the balloon. [Ref: 1 ] 

7 ] = 2S*C 



Fig. Prob.5.18 


We use eqn. (5.14): 
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Here, we have: all K.E. and P.E. changes are negligible, Q cv = 0, use £ h’ for gas flowing, and £ u for gas 
confined to the control volume. Work done is positive since the boundary of balloon is expanding. 

EES Solution: 

“Data:” 

PI = 100 “kPa” 

T1 = 22+273 “K” 

Voll = 65 “m A 3” 

P2 = 150 “kPa” 

Vol2 = Voll * (p2/pl) ££ m A 3 ... since volume is proprtional to pressure” 

P3 = 150 ££ kPa” 

T3 = 25+273 ££ K” 
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“Calculations:” 

Q + m_i * h3 = W + (m2*u2 - ml*ul) “First Law for this case of filling a control volume, see Eqn. 5.14” 

h3=Enthalpy(Hehum,T=T3,P=P3)“kJ/kg” 

u 1 =IntEnergy( Helium, T=T 1 ,P=P 1 ) “kj/kg” 

u2=IntEnergy(Helium,T=T2,P=P2)“kJ/kg” 

m_i = m2 - ml“kg.... mass entering the c.v.” 

vl=Volume(Helium,T=Tl,P=Pl)“..m A 3/kg ... sp. vol. in state 1” 

v2=Volume(Helium,T=T2,P=P2)“..m A 3/kg ... sp. vol. in state 2” 

m2 = Vol2/v2 ££ kg...mass in state 2” 

ml = Voll/vl ££ kg...mass in state 1” 

W = ((PI + P2)/2)*(Vol2-Voll)“kJ...work done by the c.v., since Vol is proportional to pressure” 

Q=0 “...heat going into the c.v. is zero” 


Solution: 


Unit Settings: SI K kPa kJ mass deg 

ml =10.6 [kg] 


h3 = -0.2942 [kJ/kg] 
rrij = 1 0.49 
P3 =150 [kPa] 

T2 = 333.6 [K] 
u2 = -503.9 [kJ/kg] 
Voll = 65 [m 3 ] 


PI =100 [kPa] 

Q =0 [kJ] 

T3 = 298 [K] 
vl = 6.1 31 [m 3 /kg] 
Vol2 = 97.5 [m 3 ] 


m2 =21.09 [kg] 

P2 =150 [kPa] 

T1 = 295 [K] 
ul = -629.1 [kJ/kg] 
v2 = <1.622 [m 3 /kg] 
W = “1063 [kJ] 


Thus: 

Final temp, T2 = 333.6 K ...Ans. 
Work done, W = 4063 kJ ...Ans. 
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(b) If the final pressure P2 varies from 120 to 150 kPa, plot the variation of T2 and W against P2: 
First, compute the Parametric Table: 



Now, plot the results: 



P2 (kPa) 
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P2 (kPa) 
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a Prob.5.19. A 100-L rigid tank contains carbon dioxide gas at 1 MPa, 300 K. A valve is cracked open, 
and carbon dioxide escapes slowly until the tank pressure has dropped to 500 kPa. At this point, the 
valve is closed. The gas remaining inside the tank may be assumed to have undergone a polytropic 
expansion, with polytropic exponent n =1.15. Find the final mass inside and the heat transferred to the 
tank during the process. [Ref:2]” 

EES Solution: 

“Data:” 

PI = 1000 “kPa” 

T1 = 300 “K” 

Voll = 0.1 “m A 3” 
n = 1.15 “... polytropic index” 

P2 = 500 “kPa” 

Vol2 = Voll“m A 3” 

T2/T1 = (P2/Pl) A ((n - l)/n)“[K]... finds temp after polytropic expansion, state 2” 

“Calculations:” 

“Note: Enthalpy of C02 exiting goes on varying from state 1 to state 2. 

So, we take the average value of enthalpy:” 

hl=Enthalpy(CarbonDioxide,T=Tl,P=Pl)“kJ/kg .... enthalpy in state 1” 
h2=Enthalpy(CarbonDioxide,T=T2,P=P2) £C kJ/kg .... enthalpy in state 2” 
h_avg = ( hi + h2 ) / 2“kJ/kg ... average enthalpy of exiting C02” 

Q = m_e * h_avg + (m2 * u2 - ml* ul) + W “...First Law for this case of filling a control volume, see 
Eqn. 5.14” 

W = 0 “..no work done, since volume is const.” 
ul=IntEnergy(CarbonDioxide,T=Tl,P=Pl)“kJ/kg. . ..internal energy” 
u2=IntEnergy(CarbonDioxide,T=T2,P=P2)“kJ/kg ... internal energy” 
m_e = (ml - m2)“kg.... mass exiting the c.v.” 

vl=Volume(CarbonDioxide,T=Tl,P=Pl) ££ ..m A 3/kg ... sp. vol. in state 1” 
v2=Volume(CarbonDioxide,T=T2,P=P2) ££ ..m A 3/kg ... sp. vol. in state 2” 
m2 = Vol2/v2 C£ kg...mass in state 2” 
ml = Voll/vl ££ kg...mass in state 1” 
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Results: 

Unit Settings: SI K kPa kJ mass deg 


hi = -7.942 [kJ/kg] 
ml = 1 .853 [kg] 
n = 1 .1 5 
Q = 24.56 [kJ] 
ul = -61 .76 [kJ/kg] 
v2 = 0.1 001 [m 3 /kg] 
W = 0 [kJ] 


h2 = -25.75 [kJ/kg] 
m2 = 0.9993 [kg] 
PI = 1000 [kPa] 
T1 = 300 [K] 
u2 = -75.79 [kJ/kg] 
Voll = 0.1 [m 3 ] 


h avg =-16.35 [kJ/kg] 
m e = 0.8587 
P2 = 500 [kPa] 

T2 =274.1 [K] 
vl = 0.05382 [m 3 /kg] 
Vol2 = 0.1 [m 3 ] 


Thus: 

Q = 24.56 kJ .... Heat transferred .... Ans. 

m2 = 0.9993 kg .... Final mass inside the tank ... Ans. 

T2 = 274.1 K ...Final temp of gas inside the tank.... Ans. 


(b) As the final pressure (P2) varies from 300 kPa to 600 kPa, plot the variation of T2, m2, and Q 
with P2: 

First, compute the Parametric Table: 
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Now, plot the results: 



P2 [kPa] 
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P2 [kPa] 



P2 [kPa] 
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5.3 Problems solved with The Expert System for Thermodynamics (TEST): 

Nozzles and Diffusers: 

Prob.5.20. Superheated vapour Ammonia enters an insulated nozzle at 20 C, 800 kPa, shown in Fig. 
below, with a low velocity and at a steady rate of 0.01 kg/s. The Ammonia exits at 300 kPa with a velocity 
of 450 m/s. Determine the temperature (or quality, if saturated) and the exit area of the nozzle. [Ref. 2] 



TEST Solution: 

1. Start TEST by going to www. ther mofluids . net , enter the required e-mail address and 

password and, we get the greeting screen. Locate the Daemons tab at the bottom of screen: 



•ipzSi r r j / Thbir/n vdyrjiiwh^-J**^ 


Greetings Dr. Muliya, registered user since 2005 


■339 


Login Message: You are now logged in using your personal license (valid until Thu Dec 17 
2015). If you encounter any problem in accessing a resource , please write to us at 
support&thermofiuids. net. 

As an educator, you can now obtain a oroup license and distribute TEST accounts to your 
students, who can use TEST from anywhere, not just from within the campus. Please take a 
look at the key problems in the Problems module. If you assign some of those as homework 
(students, of course, cannot see the answers), TEST can monitor progress and generate 
homework report. 
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2. Click on Daemons tab to get the Daemon tree to select the required Daemon: 


Daemons (Thermodynamic Calculators) 


I 

Basic Tools 


r 


1 


System Analysis 

, I , 


States & Properties 


Unit DeskCal Tables & 
Convener Charts 


Open 


Unsteady Process 


Closed 

1 

Steady Slate Steady State 

(cycles) I 


I 

Uniform 

System 


1 

Uniform 

Flow 


1 

Unsteady Process 


Generic 


Specific 


Generic 


T 


1 


Reciprocating Cycles HVAC/Psychrometry Comb- & Equilibrium 


I 

Uniform 

System 


I 

Non-Mixing 

Non-Uniform 


Semi- Mixing 
Non-Uniform 


Mixing 

Non-Uniform 



Single- Flo 




I 

Non-Mixing Multi-Flow 


I 

Mixing Multi Flow 


I 

Specific 


Vapor and Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


HVAC 

Psychrometry 


I 

Gas 

Dynamics 


? 


Combustion 
3 Equilibrium 


ick on a node to jump to a particular class of daemons, 
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3. Hover the mouse pointer over System Analysis - Open - Generic - Single Flow shown 
above. We get the following description screen: 



Click to goto pager TEST> Da emorks> System s>Open> Steady >Generic> Single-Flow System! 


Single- Flow Steady Systems: Analyze 
an open steady system with a single inlet 
and a single exit. Examples include 
turbines, compressors, pumps, nozzles, 
diffusers, throttling valves, etc. 

Chapters 4 and 6 deal with generic open 
steady systems. 


We see from the description that this is the daemon to be used for calculations regarding turbines, 
compressors, pumps, nozzles, diffusers, throttling valves etc. i.e. for most of this chapter, we will be 
using this daemon. 


4. Now, click on Single-Flow. A window appears where we have to choose the required 

Material model. In the present case, we deal with vapour/liquid Ammonia; so, we choose PC 
(i.e. Phase Change) model as shown below: 



Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20 3C as opposed to H20), which uses 
compressed liquid table for better accuracy. 

\JK_ 


Launches the PC Model 20 , R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 

Generic, Single-Flow, Open-Steady Daemon 

PC Mode 


jjusaiuiiiLy ui a jjiiasc u ai ibiuui i. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


Pure Solid and Pure Liquid: Constant density and constant specific heats ( c p = c v = c) characterize the 
solid/liquid ( SL ) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 

Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 

Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation (pv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 


c p = const 
PG Model 


pvysRT 
RG Model 


not 

SL Model 
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5. On clicking PC model, we get the following screen, where we have chosen Ammonia 
as the fluid. We have also filled up the data for State 1 as pi = 800 kPa, T1 = 20 C and 
mdotl = 0.01 kg/s. Veil = 0 by default. Click on Calculate and rest of the calculations for 
state 1 are completed in the screen shot shown below: 



6. Similarly, choose State 2 and fill in the given data, i.e. p2 = 300 kPa, Vel2 = 450 m/s and 
mdot2 = mdotl. Click on Calculate: 



Note that no new calculations are made since data is not enough; however, we will return to this State 
after entering other data: 
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7. Go to Device Panel and fill in the known data, i.e. Q = 0, Wdot_ext = 0 for a Nozzle, and 
click on Calculate. We get: 



8. Now, the important step: click on SuperCalculate to update all related States calculations. 
We get: 
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9. Now, go to State Panel and see States 1 and 2: 


State 1: 



State 2: 



Thus: A2 = 8.5668935E-6 m A 2 = 8.567 mm A 2, T2 = -9.259 C, x2 = 0.94729 Ans. 
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10. Draw the indicative T-s diagram after choosing the required plot from Plots tab: 

T,K 

446.0 


145.96 

- 1.72 s, kJ/kg.K 8.38 


(- 1 . 146 , 515.036 ) 



11. Now, go to the I/O panel and see the TEST code which can be used to regenerate the 
calculations at a later date, and also other calculations such as property values at different 
States etc: 


# OUTPUT OF SUPER- CALCULATE (starts from your inputs) 

OPERATION (for further iteration use SUPER-ITERATE) 

#*******ANALYST: Dr. Muliya; TEST License: Professional******* 

# Solution logged at: Dec 20, 2013 9:43:40 PM 

#*******TEST- CO de: To save the solution, copy the codes generated below into a text file. To reproduce 
the solution at a later time, launch the daemon (see path name below), paste the saved TEST-code at 
the bottom of this I/O panel, and click the Load button. 
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# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PC-Model; v-10.bb06 

# Start of TEST-code 


States { 

State- 1: Ammonia(NH3); 

Given: { pl= 800.0 kPa; Tl= 20.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; mdotl= 0.01 kg/s; } 
State-2: Ammonia(NH3); 

Given: { p2= 300.0 kPa; Vel2= 450.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# End 0 f TEST-code 
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# ****** DEXAILED OUTPUT: 

# Evaluated States: 

# State- 1: Ammonia(NH3) > Superheated Vapor; 

# Given: pl= 800.0 kPa; Tl= 20.0 deg-C; Vell= 0.0 m/s; 

# zl= 0.0 m; mdotl= 0.01 kg/s; 

# Calculated: vl= 0.1614 m A 3/kg; ul= 1335.6714 kj/kg; hl= 1464.7263 kj/kg; 

# sl= 5.1321 kJ/kg.K; el= 1335.6714 kj/kg; jl = 1464.7263 kj/kg; 

# Voldotl= 0.0016 m A 3/s; Al= 161.4154 m A 2; MM1= 17.031 kg/kmol; 

# State-2: Ammonia(NH3) > Saturated Mixture; 

# Given: p2= 300.0 kPa; Vel2= 450.0 m/s; z2= 0.0 m; 

# mdot2= “mdotl” kg/s; 

# Calculated: T2= -9.2592 deg-C; x2= 0.9473 fraction; y2= 0.9998 fraction; 

# v2= 0.3855 m A 3/kg; u2= 1248.0526 kj/kg; h2= 1363.4763 kj/kg; 

# s2= 5.1985 kJ/kg.K; e2= 1349.3026 kj/kg; j2= 1464.7263 kj/kg; 

# Voldot2= 0.0039 m A 3/s; A2= 0.0 m A 2; MM2= 17.031 kg/kmol; 

# Property spreadsheet starts: The following property table can be copied onto a spreadsheet (such 

as Excel) for further analysis or plots. 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

800.0 

293.2 

0.1614 

1335.67 

1464.73 

5.132 

# 02 

300.0 

263.9 0.9 

0.3855 

1248.05 

1363.48 

5.199 


# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= “6.64382E-4” kW/K; Jdot_net= 0.0 kW; Sdot_net= “-6.64382E- 4” kW/K; 

Note: In the above calculations, j = h + V A 2/2 + g.Z, and e = u + V A 2/2 + g.Z 
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(b) If the exit pressure varies from 100 to 500 kPa, mass flow rate remaining constant at 0.01 kg/s, 
plot the variation of T2 and A2 with p2: 

The procedure is quite simple: 

Go to State 2, change the pressure p2 to desired value and click Calculate, then click SuperCalculate. 
All values are updated. 

Do this for all desired values of p2 and separately tabulate p2, T2 and A2: 


P2 (kPa) 

T2 (K) 

A2 (mm A 2) 

100 

239.5 

24.74 

200 

254.3 

12.68 

300 

263.9 

8.567 

400 

271.2 

6.472 

500 

277.3 

5.2045 


Now, we can copy these values to EXCEL and draw the graphs: 


T2 vs P2 for Nozzle 



P2 (kPa) 
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Prob.5.21. A Diffuser shown in fig. has air entering at 100 kPa, 300 K, with a velocity of 200 m/s. The 
inlet cross-sectional area of the diffuser is 100 mm A 2. At the exit, the area is 860 mm A 2, and the velocity 
is 20 m/s. Determine the exit pressure and temp of air. [Ref. 2]: 



Fig. Prob.5.21 


TEST Solution: 

1. Go to the Daemon tree and locate System Analysis - Open - Generic - Single Flow: 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


2. Select the Perfect Gas (PG) Model (cp = const.) for Material model, since air is the 
working substance: 


PC Model 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20 sc as opposed to H20), which uses 
compressed liquid table for better accuracy. 

Working fluids such as H20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 
possibility of a phase transition. 


'W' 

SL Model 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simple st gas model which obeyes the ideal gas equation (pv = RT) and assumes 
_ Launches the PG Model eats to be constant. In the IG (ideal gas) model, specific heats are assumed 

f ro d- Generic, s m gie-Flcw, Open-Steady Daemon ^jpn 0 f temperature only. The RG (real gas) model uses generalized 


| | compressiblity charts and is useful for qases near the critical or super-critical 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


3. We get the following screen after clicking on PG model. Now, choose Air as the Working 
substance from the drop down menu. Then, enter known values of PI, Tl, Veil and A1 for 
State 1. Click on Calculate. We get: 

Generic, Open Steady, Single-Flow, Daemon: PG Model 


^j ^te£»« ^MthermofIuids.net > Daemons > Systems > Open > Steady > Generic ; bingleFlow > PG-Model 

*, m 

Move mouse over a variable to display its value with more precision. 



Note that mass flow rate is calculated as mdotl = 0.02323 kg/s. 

4. Now, go to State 2, and enter A2, Vel2 and mdot2 (= mdotl). Click on Calculate. We get: 



5. Go to Device Panel, enter b-state and f-state as State 1 and State 2 respectively as shown. 
Also, enter Qdot = 0 and Wdot_ext = 0, since heat transfer and work transfer for diffuser are 
zero. Click on Calculate. We get: 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


6. Now, click on SuperCalculte. We get: 


'"♦Mixed r SI r English < [©Case-0 v > | I* Help Messages On Super-lterate | 

State Panel I Device Panel I Exergy Par 


Super-Caiculate 


< [Device-A[1-2] 


J 


R/ Qdot 


& Wdot ext 


l P«I 


r Sdot_gen 


0.0 


3 


|oF 


|kW 


Sdot_net 

"v^fs.OE-S 


7 . Now, go back to States panel: 


And, observe States 1 and 2: 
State 1: 



State 2: 



Thus: Exit pressure, p2 = 123.93 kPa, exit temp, T2 = 319.73 K ... Ans. 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


8. Draw the indicative T-s diagram from the Plot tag, after choosing T-s plot: 


T, K 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


9. SuperCalculate also produces TEST code and the detailed property output etc. in the I/O 
panel. Part of the output is shown below: 

# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { pl= 100.0 kPa; Tl= 300.0 K; Vell= 200.0 m/s; zl= 0.0 m; Al= 100.0 mm A 2; } 
State-2: Air; 

Given: { Vel2= 20.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; A2= 860.0 mm A 2; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


-End of TEST-code 


# ****** DETAILED OUTPUT: 

# Evaluated States: 

# State- 1: Air > PG-Model; 

# Given: pl= 100.0 kPa; Tl= 300.0 K; Vell= 200.0 m/s; 

# zl= 0.0 m; Al= 100.0 mm A 2; 

# Calculated: vl= 0.861 m A 3/kg; ul= -84.2395 kj/kg; hl= 1.8565 kj/kg; 

# sl= 6.8929 kJ/kg.K; el= -64.2395 kj/kg; jl = 21.8565 kj/kg; 

# mdotl= 0.0232 kg/s; Voldotl= 0.02 m A 3/s; MM1= 28.97 kg/kmol; 

# Rl= 0.287 kJ/kg.K; c_pl= 1.0035 kJ/kg.K; c_vl= 0.7165 kJ/kg.K; 

# kl= 1.4005 UnitLess; 

# State-2: Air > PG-Model; 

# Given: Vel2= 20.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; 

# A2= 860.0 mm A 2; 

# Calculated: p2= 123.9268 kPa; T2= 319.7311 K; v2= 0.7404 m A 3/kg; 

# u2= -70.102 kj/kg; h2= 21.6565 kj/kg; s2= 6.8952 kJ/kg.K; 

# e2= -69.902 kj/kg; j2= 21.8565 kj/kg; Voldot2= 0.0172 m A 3/s; 

# MM2= 28.97 kg/kmol; R2= 0.287 kJ/kg.K; c_p2= 1.0035 kJ/kg.K; 

# c_v2= 0.7165 kJ/kg.K; k2= 1.4005 UnitLess; 
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Basic Thermodynamics: Software Solutions 


Part II I Law of Thermodynamics applied to Flow Processes 

# Property spreadsheet: 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 100.0 300.0 0.861 -84.24 1.86 6.893 

# 2 123.93 319.7 0.7404 -70.1 21.66 6.895 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= “5.472404E-5” kW/K; Jdot_net= 0.0 kW; Sdot_net= “-5.472404E- 
5” kW/K; 


(b) As A1 varies from 50 to 300 mm A 2, plot the variation of mdot, p2 and T2, other quantities 
remaining unchanged: 

The procedure is as follows 

Go to State 1, change the pressure A1 to desired value and click Calculate (or, hit Enter), then click 
SuperCalculate. All values are updated. 

Do this for all desired values of Aland separately tabulate Al, mdot, p2, and T2: 


A1 (mm A 2) 

mdot (kg/s) 

P2 (kPa) 

T2 (K) 

50 

0.01161 

61.96 

319.73 

100 

0.02323 

123.93 

319.73 

150 

0.03484 

185.89 

319.73 

200 

0.04646 

247.85 

319.73 

250 

0.05807 

309.82 

319.73 

300 

0.06969 

371.78 

319.73 


Note that T2 does not change; but, mdot and P2 vary with A2. 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


Now, plot these results in EXCEL: 


mdot vs A1 for Diffuser 



A1 (mm A 2) 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


P2 vs A1 for Diffuser 



A1 (mm A 2) 


Prob.5.22. Carbon dioxide enters an adiabatic nozzle steadily at 1 MPa, 500 C with a mass flow rate of 
6000 kg/h and leaves at 100 kPa and 450 m/s. The inlet area of the nozzle is 40 cm A 2. Determine (a) 
the inlet velocity, and (b) the exit temperature. [Ref. 1]: 



Fig. Prob.5.22 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


TEST Solution: 

1. Go to Daemon tree, choose System Analysis - Open - Generic - Single Flow as shown below: 


Basic Tools 


System Analysis 

t I , 


1 

States & Pr operties 
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Flow 
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Unsteady Process 
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(cycles) I 
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Non-Mixing Multi-Flow 


I 
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I 
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I I I 

Vapor and Gas Refrigeration HVAC Gas Combustion 

Power Cycles Cycles Psychrometry Dynamics 3 Equilibrium 


2. On clicking Single Flow, we are led to Material model: 


Pure Solid and Pure Liquid: Constant density and constant specific heats ( c p = c v = c) characterize the 
solid/liquid ( SL ) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 



SL Model 


Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed 
to maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be 
analyzed with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no 
possibility of a phase change. For specific examples, click on the help icon at the bottom margin of the 
daemon. 



Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model 
is the simplest gas model which obeyes the ideal gas equation (pv = RT) and 
assumes specific heats to be constant. In the IG (ideal gas) model, specific heats 
are assumed to be function of temperature only. The RG (real gas) model uses 
generalized compressiblity charts and is useful for gases near the critical or 
super-critical conditions for which PC-model data are not avaiable. 


Helium expands steadily in a nozzle from an inlet-state to an exit-state 
possibility of phase change. For specific examples, click on the help icon at 


the bottom margin of the daemon. 
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I Law of Thermodynamics applied to Flow Processes 


3. Choose Ideal Gas (IG) model and select C02 as the working fluid. In IG Model, sp. heat 
is taken as a function of temp. Enter the data given for State 1, i.e. PI, Tl, A1 and mdotl; 
click on Calculate (or, hit Enter). We get: 


Generic, Open Steady, Single-Flow, Daemon: 7G Model 


thermoftuids.net > Daemons > Systems > Open > Steady > Generic > SingleFIow > IG-Model 


3 H T ' 


Home of | P ^ 


i rnr 


nil m 


Move mouse over a variable to display its value with more precision. 



Note that Veil is calculated as Veil = 60.86 m/s.... Ans. 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


4. Enter data known for State 2, i.e. P2, Vel2 and mdot2 = mdotl. Hit Enter. We get: 


Move mouse over a variable to display its value with more precision. 
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5. Go to Device Panel. Enter State 1 for b-state, State 2 for f-state, and Qdot = 0, Wdot_ext = 0. 
Click on Calculate. We get: 



6. Now, click on SuperCalculate. We get: 
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Part II 


I Law of Thermodynamics applied to Flow Processes 


7 . Go back to States Panel: We get: 


State 1: 



And, State 2: 



Thus: Inlet velocity, Veil = 60.86 m/s, exit temp, T2 = 412.76 C ... ans. 
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Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


8. Plot the indicative T-s diagram: 
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Basic Thermodynamics: Software Solutions 
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I Law of Thermodynamics applied to Flow Processes 


9. The I/O panel gives TEST code and other details. Part of it is shown below: 


# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: 002; 

Given: { pl= 1000.0 kPa; Tl= 500.0 deg-C; zl= 0.0 m; mdotl= “6000/3600” kg/s; Al= 40.0 
cm A 2; } 

State-2: C02; 

Given: { p2= 100.0 kPa; Vel2= 450.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# p; n d of TEST-code 

# Property spreadsheet starts: # 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 1000.0 773.2 0.1461 -8600.73 -8454.67 5.379 

# 2 100.0 685.9 1.2958 -8683.64 -8554.07 5.677 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Device-A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= 0.49769834 kW/K; Jdot_net= “-2.1262167E-5” kW; Sdot_net= 
-0.49769834 kW/K; 
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I Law of Thermodynamics applied to Flow Processes 


Prob.5.23. Refrigerant 134a at 700 kPa and 120 C enters an adiabatic nozzle with a velocity of 20 
m/s and leaves at 400 kPa and 30 C. Determine (a) the exit velocity, and (b) ratio of inlet to exit area, 
A1/A2. [Ref. 1] 


R1 34a 



TEST Solution: 

1. In the Daemons tree, select System Analysis - Open - Generic - Single Flow as 
shown below: 
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I Law of Thermodynamics applied to Flow Processes 


2. For Material model, choose Phase Change (PC) model, since R134a is the working fluid: 



Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

L&: 


Launches the PC Model ?o, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 

Generic, Single-Flow, Open-Steady Daemon 

PC Mode 
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Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 

Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 

Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation (jjv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 
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3. In the Window that appears, select the R134a as the substance and enter the data for State 1, 
i.e. PL T1 and Veil, and click Calculate (or, hit Enter). We get: 


Generic, Open Steady, Single-Flow, Daemon: PC Model 


thermofluids.net > Daemons > Systems > Open > Steady > Generic > SingleFlow > PC-Model 

xvrmm *, m 

Move mouse over a variable to display its value with more precision. 



4. Similarly, fill in data for State 2, press Enter: 



5. Go to Device Panel, enter Qdot = 0, Wdot_ext = 0; press Enter: 
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I Law of Thermodynamics applied to Flow Processes 


6. Click on SuperCalculate. Go to States Panel: 


State 2: 



Thus: Exit velocity = Vel2 = 409.53 m/s ... Ans. 

7. Use the I/O panel to calculate A1/A2: 

We have: rhol * A1 * Veil = rho2 * A2 * Vel2.... By mass balance 
i.e. (A1/A2) = ( 1 / v2 ) * Vel2/ ( Vel 1 * (1/vl)) 

i.e. (A1/A2) = 15.578 .... Ans. 


8. Also, from I/O panel, copy the TEST code etc: 


# OUTPUT OF SUPER-CALCULATE 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: R-134a; 

Given: { pl= 700.0 kPa; Tl= 120.0 deg-C; Vell= 20.0 m/s; zl= 0.0 m; } 
State-2: R-134a; 

Given: { p2= 400.0 kPa; T2= 30.0 deg-C; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# p; n d of TEST-code 
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# Property spreadsheet starts: 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

700.0 

393.2 

0.0433 

327.57 

357.89 

1.189 

# 02 

400.0 

303.2 

0.0569 

251.46 

274.23 

0.99 


# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

#******CalCULATE VARIABLES: Type in an expression starting with an sign ( c = mdotl*(h2-hl) 5 , 
£ = sqrt(4*Al/PI)’ etc.) and press the Enter key)********* 

#(A1/A2) = ( 1 /v2) * Vel2/ ( Vel 1 * (1/vl)) 

( l/v2)*Vel2/ (Veil * (1/vl)) = 15.578052777777797 
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Basic Thermodynamics: Software Solutions 
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I Law of Thermodynamics applied to Flow Processes 


Prob.5.24. Air at 80 kPa, 27 C, and 220 m/s enters a Diffuser at a rate of 2.5 kg/s and leaves at 42 C. 
The exit area of the diffuser is 400 cm A 2. The air is estimated to lose heat at a rate of 18 kj/s during this 
process. Determine: (a) the exit velocity, and (b) the exit pressure. [Ref. 1] 



Fig. Prob.5.24 


TEST Solution: 

1. Go to Daemon tree, choose System Analysis - Open - Generic - Single Flow as shown below: 


Basic fools 


System Analysis 

, I , 
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Unit DeskCal Tables 
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I Law of Thermodynamics applied to Flow Processes 


2. For Material model, choose Ideal Gas (IG) model, where cp is taken as a function of temp. 
(PG model also will give almost the same results): 


PC Model 


Working fluids such as H20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is 
any possibility of a phase transition. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the 
help icon at the bottom margin of the daemon. 



SL Model 


Pure Solid and Pure Liquid: Constant density and constant specific heats ( c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed 
to maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be 
analyzed with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no 
possibility of a phase change. For specific examples, click on the help icon at the bottom margin of the 
daemon. 


PG Model 


pv 

IG I 


pv-zRT 
RG Model 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model 
is the simplest gas model which obeyes the ideal gas equation (pv = RT) and 
assumes specific heats to be constant. In the IG (ideal gas) model, specific heats 
are assumed to be function of temperature only. The RG (real gas) model uses 
generalized compressiblity charts and is useful for gases near the critical or 
super-critical conditions for which PC-model data are not avaiable. 


^^ C ri*,Vingle-Ftow! e Open-Steady Daemon |’P les: Helium ex P ands steadil V in a nozzle from a " ******* to an exit-state 

no possibility of phase change. For specific examples, click on the help icon at 


rwTcrfr 

I 1-hp hnl+nm marnin nf 1 Tip Hapmnn 


3. Choose Air as the working substance. Enter given data of PI, Tl, Veil, and mdotl for State 
1. Click Enter: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 


^^^^ ^^ ^ ^ ^^Bthermofluids.net > Daemons > Systems > Open > Steady > Generic > SingleFlow > IG-Model 

■"■lEMPi *, m 

Move mouse over a variable to display its value with more precision. 



Note that in addition to properties such as ul, hi, si, volume flow rate Voldotl and inlet area A1 are 
also calculated. 
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I Law of Thermodynamics applied to Flow Processes 


4. Similarly, enter data for State 2, i.e. T2, A2 and mdot2, and click Enter: 



5. Go to Devices Panel, enter Qdot = -18 kW (negative sign since heat is leaving the system), 
Wdot_ext = 0. And for i-State = State 1, b-state = State 2. Click Calculate (or, hit Enter): 
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I Law of Thermodynamics applied to Flow Processes 


6. Now, click on SuperCalculate: 


Go to State- 1 and 2: 
State 1: 



State 2: 


rho2 = 1 .0069549 kg/m A 3 [Density] 


1 <*■ Mixed 

rs\r 

English ^<J 

©Case-0 v > 

I 17 Help P/lessages On 

Super-lterate | 

Super Calculate 1 oad | 

| Super-lrwtialize | 

State Panel 
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- 

i i 

■ < | ©State-2 v]>| 


B No-Plots 

v Initialize 

J 1 


Fno Yes Air 


P2 


17 72 


r rho2 


r v2 

r u 2 


[91.07284 

[kPa 

v 142.0 


3 [1.00695 

kg/m*3 

v| [0.99309 

m^g v] [-73.34679 

kJ/kg ^ 

r h2 


r s 2 


r Vel2 


17 z2 

r e2 


[17.09702 

kJ/kg 

v| |6. 96924 

kJ/kg.K 

V |62.06832 

nVs 

51 oo 

■ m v 1-71.42055 

j kJ/kg v 

17 j 2 


phi2 


psi2 


17 mdot2 

r Voldot2 

19.023258 

_ kJ/kg 

Jill 

kJ/kg 


kJ/kg 

v =mdot1 

■ kg/s v |2.48273 

(1^3/8 v 

17 >42 


MM2 


R2 


r c _p2 



400.0 

■i”« 

v| |28.97 

kg/kmol 

v] jo. 28699 

kJ/kg.K 

vj|l. 00605 

kJ/kg.K v 



Thus: Vel2 = 62.07 m/s, p2 = 91.07 kPa Ans. 
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I Law of Thermodynamics applied to Flow Processes 


7. Draw the indicative T-s diagram: 

T,K 

346.67 


270.14 

6.26 s, kJ/kg.K 7.67 

8. From the I/O panel, copy the TEST code etc: 

# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { pl= 80.0 kPa; Tl= 27.0 deg-C; Vell= 220.0 m/s; zl= 0.0 m; mdotl= 2.5 kg/s; } 
State-2: Air; 

Given: { T2= 42.0 deg-C; z2= 0.0 m; mdot2= “mdotl” kg/s; A2= 400.0 cm A 2; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= -18.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


( 6 . 379 , 363.601 ) 



# E n d of TEST-code 
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# Property spreadsheet starts: 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

80.0 

300.2 

1.0767 

-84.12 

2.02 

6.957 

# 

2 

91.07 

315.2 

0.9931 

-73.35 

17.1 

6.969 


Prob.5.25. Argon gas enters steadily an adiabatic turbine at 900 kPa and 450 C with a velocity of 80 m/s 
and leaves at 150 kPa with a velocity of 150 m/s. The inlet area of the turbine is 60 cm A 2. If the power 
output of the turbine is 250 kW, determine the exit temp of argon. [Ref. 1] 


PI, T1, Veil, A1 
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I Law of Thermodynamics applied to Flow Processes 


TEST Solution: 

1. Go to System Analysis ---Single Flow daemon as shown below: 


Basic Tools 

■ — i — 


System Analysis 

1 , 


S tates 8t Pro perties 


Closed 
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Uniform 

System 


1 

Uniform 

Flow 
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Unsteady Process 
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I I I 1 
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Single-Flo 




I 

Non-Mixing Multi- Flow 


I 

Mixing Mult if low 


Specific 


T 


T 


I 

Vapor mid Gas Refrigeration HVAC Gas Combustion 

Power Cycles Cycles Psychro merry Dynamics & Equilibrium 


2. For Material Model, choose Perfect Gas (PG) model, and select Argon for the working 
substance. Enter the data, viz. PI, Tl, Veil and A1 for State 1, and click on Calculate 
(or, hit Enter). We get: 


Generic, Open Steady, Single-Flow, Daemon: PG Model 


^ > Systems > Open > Steady > Generic > SingleFlow > PG-Model 

■■■■ mi m & mm *, as 

Move mouse over a variable to display its value with more precision. 



Note that additional properties at State 1 and mass flow rate, mdotl are calculated. 
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I Law of Thermodynamics applied to Flow Processes 


3. Enter given data for State 2, i.e. P2, Vel2 and modot2 = mdotl; hit Enter: 



4. Go to Devices Panel. Enter Qdot = 0 and Wdot_ext = 250 kW. Enter State 1 and State 2 for 
b-state and f-state respectively. Click Calculate: 


Move mouse over a variable to display its value with more precision. 


(• Mixed C SI ^English < [ Case-0 C] > | P Help Messages On 
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State-Null: 
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5. Click SuperCalculate: 
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I Law of Thermodynamics applied to Flow Processes 


6. Now, go back to States Panel, see State 2: 



Observe that T2 = 267.14 deg. C Ans. 
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I Law of Thermodynamics applied to Flow Processes 


7. Indicative T-s diagram, obtained from Plots tab, is as follows: 



3.49 


s, kJ/kg.K 


4.51 


8. The I/O panel gives the TEST code, and other details: 

# OUTPUT OF SUPER-CALCULATE : 

# start of TEST-code 

States { 

State- 1: Ar; 

Given: { pl= 900.0 kPa; Tl= 450.0 deg-C; Vell= 80.0 m/s; zl= 0.0 m; Al= 60.0 cm A 2; } 
State-2: Ar; 

Given: { p2= 150.0 kPa; Vel2= 150.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 250.0 kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 
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# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

900.0 

723.2 

0.1672 

70.63 

221.13 

3.882 

# 

2 

150.0 

540.3 

0.7496 

13.55 

125.99 

4.104 


Prob.5.26. Air flows steadily through an adiabatic turbine, entering at 1 MPa, 500 C and 120 m/s and 
leaving at 150 kPa, 150 C and 250 m/s. The inlet area of turbine is 80 cm A 2. Determine (a) the mass 
flow rate of air, and (b) the power output of turbine. [Ref. 1] 


TEST Solution: 

1. Choose System Analysis Single Flow daemon as in previous cases: 
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I Law of Thermodynamics applied to Flow Processes 


2. For Material model, choose Ideal Gas (IG) model, where cp is a function of temp. (We can 
choose PG model also; results will not be much different): 


I I 

PC Model 


Working fluids such as H20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is 
any possibility of a phase transition. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the 
help icon at the bottom margin of the daemon. 


SL Model 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed 
to maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be 
analyzed with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no 
possibility of a phase change. For specific examples, click on the help icon at the bottom margin of the 
daemon. 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model 
is the simplest gas model which obeyes the ideal gas equation {pv = RT) and 
assumes specific heats to be constant. In the IG (ideal gas) model, specific heats 
are assumed to be function of temperature only. The RG (real gas) model uses 
generalized compressiblity charts and is useful for gases near the critical or 
super-critical conditions for which PC-model data are not avaiable. 


c p = const 
PG Model 


pv- zRT 
RG Model 


pv- 

IG f 


Launches the IG Model 
Generic, Single-Flow, Open-Steady Daemon 

" I II 
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no possibility of phase change. For specific examples, click on the help icon at 
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I Law of Thermodynamics applied to Flow Processes 


3. Choose Air for material, enter data, i.e. PI, Tl, Al, Veil for State 1; press Enter: 

Generic, Open Steady, Single-Flow, Daemon: IG Model 


> Daemons > Systems > Open > Steady > Generic > SingleFlow > IG-Model 

I eesi m as *r m ■ *, se 

Move mouse over a variable to display its value with more precision. 



Note that mdotl is calculated as 4.327 kg/s .... Ans. 

4. Now, enter data for State 2, i.e. P2, T2, Vel2, and mdot2 = mdotl. Press Enter: 



Note that A2 is calculated as 140.11 cm A 2. 

5. Now, go to Device Panel. Enter State 1 and State 2 for b-state and f-state respectively. Also 
Qdot = 0. Press Calculate: 
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I Law of Thermodynamics applied to Flow Processes 


Note that work output is calculated as: Wdot_ext = 1495.3 kW .... Ans. 

6. Indicative T-s diagram is as follows: 



6.41 s, kJ/kg.K 7.94 


7. Clicking on SuperCalculate gives TEST code etc. in the I/O panel: 


# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: Air; 

Given: { pl= 1000.0 kPa; Tl= 500.0 deg-C; Vell= 120.0 m/s; zl= 0.0 m; Al= 80.0 cm A 2; } 
State-2: Air; 

Given: { p2= 150.0 kPa; T2= 150.0 deg-C; Vel2= 250.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 

# 
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# Property spreadsheet starts:. 

# 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

1000.0 

773.2 

0.2219 

274.45 

496.33 

7.215 

# 

2 

150.0 

423.2 

0.8096 

5.24 

126.68 

7.125 


# 

# Property spreadsheet ends 

# 

# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; T_B= 298.15 K; 

# Calculated: Wdot_ext= 1495.2935 kW; Sdot_gen= -0.3899593 kW/K; Jdot_net= 
1495.2935 kW; Sdot_net= 0.3899593 kW/K; 


Prob.5.27. Steam flows steadily through an adiabatic turbine. The inlet conditions of steam are: 6 MPa, 
400 C and 80 m/s and the exit conditions are: 40 kPa, 92% quality and 50 m/s. The mass flow rate of 
steam is 20 kg/s. Determine: (a) the change in K.E. (b) the power output, and (c) the turbine inlet area. 
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(b) Plot the Power output and exit temp against the exit pressure as exit pressure varies from 10 to 200 
kPa. [Ref. 1] 


V] 



JC-y 

vl 


Fig.Prob.5.27 


TEST Solution: 

1. Go to Daemons tree and choose System Analysis Single Flow as shown below: 
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2. Choose PC model for Material Model, since Steam is the working substance: 


Select a material model to launch the single-flow open-steady daemon. 




Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 


L&. 



compressed-liquid sub-model, except for species with an asterisk (H20 ac as opposed to H20), which uses 
compressed liquid table for better accuracy. 


PC Mode 



\A/nrlrinn fli lirlc ci irh ac H7H R-1 7 NH3. R-1 Ida N? CH7 ptr should hp trpstpri as PT fluids if thprp is snv 


Launches the PC Model 

L • • 


Generic, Single-Flow, Open-Steady Daemon 

isition. 




Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


3. Choose H20 for working substance, enter data for State 1, i.e. enter PI, Tl, Veil and mdotl, 
and press Enter (or, Calculate): 

We get: 


Generic, Open Steady, Single-Flow, Daemon: PC Model 


i:v . • -’H rmofluids.net systems > Open > Steady > Generic > SingleFlow > PC-Model 

*, m 

Move mouse over a variable to display its value with more precision. 



Note that Turbine inlet area A1 is calculated as: A1 = 0.01185 m A 2 ... Ans. 

4. Enter data for State 2, i.e. P2, x2, Vel2 and mdotl = mdot2. Press Enter (or, Calculate). We get: 
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5. Go to Device Panel. Enter State 1 and State 2 for b-state and f-state respectively. Also enter 
Qdot = 0, and press Calculate. We get: 



Note that Wdot_ext is calculated as: 14558.3 kW = 14.558 MW = Work output of turbine... Ans. 


Also: Vel2 = 50 m/s, Veil = 80 m/s, and therefore, change in K.E. = (Vel2 A 2 - Vell A 2) / 2. 
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i.e. 


YelT - Veil 



6. Click on Super Calculate and go to I/O panel to get TEST code and other details: 


#• 


OUTPUT OF SUPER-CALCULATE: 


# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PC-Model; v-10.ca08 


#■ 


Start of TEST-code 


States { 

State- 1: H20; 

Given: { pl= 6000.0 kPa; Tl= 400.0 deg-C; Vell= 80.0 m/s; zl= 0.0 m; mdotl= 20.0 kg/s; } 
State-2: H20; 

Given: { p2= 40.0 kPa; x2= 0.92 fraction; Vel2= 50.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; T_B= 298.15 K; } 

} 

# p; n d of TEST-code 

# Property spreadsheet starts: 

# State p(kPa) T(K) x v(m3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 01 6000.0 673.2 0.0474 2892.81 3177.14 6.541 

# 02 40.0 349.0 0.9 3.682 2304.23 2451.18 7.139 

# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; T_B= 298.15 K; 

# Calculated: Wdot_ext= 14558.326 kW; Sdot_gen= 11.966095 kW/K; Jdot_net= 
14558.326 kW; Sdot_net= -11.966095 kW/K; 
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(b) Plot Power output and T2 as P2 varies from 10 to 200 kPa: 

The procedure is quite simple: 

Go to State 2, enter the desired value of P2 and press Enter. Then, press SuperCalculate. Read the value 
of T2 and Wdot_ext and tabulate the values. Results are shown below: 


P2(kPa) 

Wdot_ext (MW) 

T2 (deg.C) 

10 

15.717 

45.81 

30 

14.815 

69.08 

50 

14.354 

81.31 

80 

13.907 

93.48 

110 

13.59 

102.3 

140 

1 3.345 

109.29 

170 

13.143 

115.17 

200 

12.972 

120.23 


Plot the above results in EXCEL: 
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Prob.5.28. Air enters the compressor of a gas turbine plant at 100 kPa, 25 C with a low velocity and 
exits at 1 MPa and 347 C with a velocity of 90 m/s. The compressor is cooled at a rate of 1500 kj/min 
and the power input to the compressor is 250 kW. Determine the mass flow rate of air through the 
compressor. [Ref. 1] 


‘tan 



Fig. Prob.5.28 


TEST Solution: 

1 . Go to System Analysis Single Flow daemon, as in the case of previous problems: 


Basic f ools 


System Analysis 

. I . 


States & Properties 


Unit DeskCal Tables <6 
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Closed 


Open 

I 


I 

Uniform 

System 


1 

Uniform 

Flow 


Unsteady Process 


Steady State Steady State 
(cycles) I 


1 

Unsteady Process 


Generic 


r 


Specific 


Generic 


T 


1 


Reciprocating Cycles HVAC/Psychrometry Comb. & Equilibrium 


I I I 1 

Uniform Non-Mixing Semi-Mixing Mixing 

System Non-Uniform Non-Uniform Non-Uniform 



Single-Flo 


I 

Non-Mixing Multi-Flow 


I 

Mixing Multi-Flow 


Specific 


T 


T 


I 
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2. Choose the Ideal Gas (IG) model, since we are going to use Air as working substance: 



SL Model 


Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which 
can be assumed to maintain their condensed (solid or liquid) phase when a system 
undergoes other changes, can be analyzed with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state 
with no possibility of a phase change. For specific examples, click on the help icon at the 
bottom margin of the daemon. 


c„ = const 
PG Model 


pv-zRT 
RG Model 


pv=RL_ 

Launches the IG Model 
IG Mo Generic, Single-Flow, Open-Steady Daemon 


Pure Gas: A pure gas has a fixed chemical composition across 
space and time. Oxygen, nitrogen, and air are examples of a pure 
gas. The PG (perfect gas) model is the simplest gas model which 
obeyes the ideal gas equation (pv = RT) and assumes specific 
heats to be constant. In the IG (ideal gas) model, specific heats 
are assumed to be function of temperature only. The RG (real 
gas) model uses generalized compressiblity charts and is useful 
for gases near the critical or super-critical conditions for which PC- 
model data are not avaiable. 

Examples: Helium expands steadily in a nozzle from an inlet-state 
to an exit-state with no possibility of phase change. For specific 
examples, click on the help icon at the bottom margin of the 
daemon. 


3. Choose Air for working substance, enter data i.e. PI, Tl, Veil for State 1 and press Enter: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 


thermofluids.net > Daemons > Systems > Open > Steady > Generic > SinqleFlow > IG-Model 

msssbI osm m a? *r m m *. m 

Move mouse over a variable to display its value with more precision. 



4. Enter P2, T2 and Vel2 and mdot2 = mdotl for State 2; press Enter: 
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5. Go to Device Panel, enter b-state and f-state, and also Qdot = 0, Wdot_ext = -250 (negative 
sign since work is input to compressor), press Enter: 



6. Now, click on SuperCalculate to up-date all calculations: 


Go to State Panel, State 1: 



And, State 2: 



Thus: mdotl = 0.6699 kg/s .... Ans. 
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7. Indicative T-s diagram is as follows: 
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8. From the I/O panel, get the TEST code etc.: 

# OUTPUT OF SUPER-CALCULATE: 


# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: Air; 

Given: { pl= 100.0 kPa; Tl= 25.0 deg-C; Veil- 0.0 m/s; zl= 0.0 m; } 

State-2: Air; 

Given: { p2= 1000.0 kPa; T2= 347.0 deg-C; Vel2= 90.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= “-1500/60” kW; Wdot_ext= -250.0 kW; T_B= 298.15 K; } 

} 


# p; n d of TEST-code 

# Property spreadsheet starts: 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 100.0 298.2 0.8557 -85.55 0.02 6.887 

# 2 1000.0 620.2 0.178 153.86 331.83 6.978 

# 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= “-1500/60” kW; Wdot_ext= -250.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen=0. 145071 16kW/K;Jdot_net= -225.0 kW;Sdot_net= -0.061220754 
kW/K; 
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Prob.5.29. A compressor operating at steady state takes in 45 kg/min of methane gas (CH4) at 1 bar, 25 
C, 15 m/s, and compresses it with negligible heat transfer to 2 bar, 50 m/s at exit. The power input to 
the compressor is 1 10 kW. Using the ideal gas model, determine the temp of the gas at the exit. [Ref. 5] 


TEST Solution: 

1. Go to System Analysis ... Single Flow daemon as shown: 


Basic fouls System Analysis 


Closed Open 


Unsteady Process Steady State Steady State 

I . (cycles) I 


States 8t Properties 

I 1 1 

Uniform Uniform 

System Flow 

1 

Unsteady Process 


Generic 

r 


Specific 


Generic - 


Specific 


T 


T 


Reciprocating Cycles HVAC/Psychrometry Comb. 3 Equilibrium 


I 

Uniform 

System 



Vapor and Gas 
Power Cycles 


I I 1 

Non-Mixing Semi-Mixing Mixing 

Non-Uniform Non-Uniform Non-Uniform 

; i 

Non-Mixing Multi-Flow Mixing Multi-Flow ^ 

I I I T 

Refrigeration HVAC Gas Combustion 

Cycles Psychrometry Dynamics & Equilibrium 


2. Select IG model for Material model: 



SL Model 


assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 


pv = RT 

IG Mode 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 

Pure Gas: A pure gas has a fixed chemical composition across space and time. 

Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation (pv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 

Launches the IG Model I 51 Helium expands steadily in a nozzle from an inlet-state to an exit-state with 

Generic, Single-Flow, Open-Steady Daemon 3 j|jty of phase change. For specific examples, click on the help icon at the 



bottom margin of the daemon. 
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3. Choose Methane (CH4) for working substance, enter data for State 1 (i.e. PI, Tl, Veil and 
mdotl), press Enter: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 




> Steady > Generic . SinqleFlow := IG-Model 



mm *,1 

Move mouse over a variable to display its value with more precision. 
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4. Enter data for State 2, (i.e. P2, Vel2 and mdot2), press Enter: 


Move mouse over a variable to display its value with more precision. 


Mixed C si English < [{Case-0 v > | [✓ Help Messages On Super-lterate | 
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5. Go to Device Panel, enter for b-state and f-state, and also Qdot = 0, Wdot_ext = -110 kW 
(negative sign since work is done on the system in compressor). Press Enter: 



6. Now click on SuperCalculate. Go to State Panel, State 2. We get: 
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Thus: T2 = 87.02 deg. C ... Ans. 

7. Indicative T-s diagram is as follows: 

T, K 

396.19 


263.34 

10.46 s, kJ/kg.K 12.36 



8. I/O panel gives TEST code etc: 


# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.ca08 

# start of TEST-code 


States { 

State- 1: Methane(CH4); 

Given: { pl= 100.0 kPa; Tl= 25.0 deg-C; Vell= 15.0 m/s; zl= 0.0 m; mdotl= 45.0 kg/min; } 
State-2: Methane(CH4); 

Given: { p2= 200.0 kPa; Vel2= 90.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/min; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= -110.0 kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 
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# Property spreadsheet starts: 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

100.0 

298.2 

1.5454 

-4822.36 

-4667.82 

11.619 

# 

2 

200.0 

360.2 

0.9334 

-4711.78 

-4525.09 

11.694 


# Property spreadsheet ends 

# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= 0.0 kW; Wdot_ext= -110.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= 0.056342352 kW/K; Jdot_net= -110.00039 kW; Sdot_net= 
-0.056342352 kW/K; 


Prob.5.30. Helium is to be compressed from 120 kPa, 310 K to 700 kPa, 430 K. A heat loss of 20 kj/kg 
occurs during compression. Neglecting K.E. changes, determine the power input required for a mass 
flow rate of 90 kg/min. [Ref. 1] 
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TEST Solution: 

1. Go to System Analysis ... Single Flow daemon as shown: 
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2. Choose the Ideal Gas (IG) model for Material model, since Helium is the 
working substance: 


SL Model 


to maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be 
analyzed with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no 
possibility of a phase change. For specific examples, click on the help icon at the bottom margin of the 
daemon. 



pv=zRT 
RG Model 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model 
is the simplest gas model which obeyes the ideal gas equation (pv = RT) and 
assumes specific heats to be constant. In the IG (ideal gas) model, specific heats 
are assumed to be function of temperature only. The RG (real gas) model uses 
generalized compressiblity charts and is useful for gases near the critical or 
super-critical conditions for which PC-model data are not avaiable. 


jq Mri( Launches the IG Model 
Generic, Single-Flow, Open-Steady Daemon 


;s: Helium expands steadily in a nozzle from an inlet-state to an exit-state 
possibility of phase change. For specific examples, click on the help icon at 


tne Dotfcom margin of the daemon. 


Binary Mixture: The mixture of two gases, A and B, is expressed in terms of the 
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3. Choose He for working substance, enter data for State 1 (i.e. PI, Tl, mdotl = 1.5 kg/s), press 
Enter. We get: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 


ttiermofluids.net > Daemons > Systems > Open > Steady > Genenc > StnqleFlow > IG-Model 

R9|miH kktbh *, m 

Move mouse over a variable to display its value with more precision. 



4. Enter data for State 2, i.e. P2, T2, mdot2 = mdotl. Press Enter. We get: 



5. Go to Device Panel. Enter Qdot = - 20 * mdotl and click on Calculate, and SuperCalculate. 
We get: 
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Thus: W = - 965.37 kW ... Ans. (negative sign, since work is done on the system in compressor) 

6. Indicative T-s diagram from Plots tab: 
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7. I/O panel gives the TEST code and other details: 


# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: He; 

Given: { pl= 120.0 kPa; Tl= 310.0 K; Vell= 0.0 m/s; zl= 0.0 m; mdotl= 1.5 kg/s; } 
State-2: He; 

Given: { p2= 700.0 kPa; T2= 430.0 K; Vel2= 0.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= “-20*mdotr kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 

# Property spreadsheet starts 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 120.0 310.0 5.3695 -582.75 61.58 31.389 

# 2 700.0 430.0 1.2768 -208.59 685.16 29.424 

# Property spreadsheet ends 


# Mass, Energy, and Entropy Analysis Results: 

# Device-A: i-State = State- 1; e-State = State-2; 

# Given: Qdot= “-20*mdotl” kW; T_B= 298.15 K; 

# Calculated: Wdot_ext= -965.37134 kW; Sdot_gen= -2.8472614 kW/K; 

Jdot_net= -935.37134 kW; Sdot_net= 2.9478817 kW/K; 
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Prob.5.31. Refrigerant- 134a is throttled from the sat. liquid state at 800 kPa to a temp of -20 C. Determine 
the pressure of the refrigerant at the final state. [Ref. 1] 

TEST Solution: 

Note that this is a problem on throttling. The daemon to be used is still the same as used earlier, viz. 
Systems>Open>SteadyState>Generic>SingleFlow>IG-Model: 

1. Go to System ... Single Flow daemon: 
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2. Choose PC model for material model since R134a is the working substance: 

Generic, Single-Flow, Open-Steady Daemons: Selects Material Model 



| a m TMll 1 \-&m 2 1 k =° 


Select a material model to launch the single-flow open-steady daemon. 


EH 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub- 
cooled (compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. 
Based on the saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is 
modeled with the compressed-liquid sub-model, except for species with an asterisk (H20 3C as opposed 
to H20), which uses compressed liquid table for better accuracy. 

PC Model Launches the PC Model R ' 12, NH3, R * 134a ' N2 ' c ° 2 > etc -' should be treated as PC fluids if there is 

Generic, Single-Flow, Open-Steady Daemon transition. 



Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the 
help icon at the bottom margin of the daemon. 


SL Model 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed 
to maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be 
analyzed with the SL model. 

Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no 
possibility of a phase change. For specific examples, click on the help icon at the bottom margin of the 
daemon. 
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3. Choose R134a for working substance and enter data for State 1, i.e. PI, xl and press Enter. 
We get: 

Generic, Open Steady, Single-Flow, Daemon: PC Model 


^ > Systems > Open > Steady > Generic > SingleFlow > PC-Model 

I EBt m w *r m m *, m 

Move mouse over a variable to display its value with more precision. 



4. Enter data for State 2, i.e. T2, and h2 = hi since it is throttling process. Click on Calculate 
and SuperCalculate. We get: 



Thus: p2 = 133.7 kPa ... Ans. 
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5. Indicative P-h diagram is easily obtained from the Plots tab: 


p, kPa (log Scale) 


( 200.945. 87050.213 ) 


9329.3 


1.5506 



- 34.42 


h, kJ/kg 


303.52 


6. I/O panel gives the TEST code etc: 


# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: R-134a; 

Given: { pl= 800.0 kPa; xl= 0.0 fraction; Vell= 0.0 m/s; zl= 0.0 m; } 
State-2: R-134a; 

Given: { T2= -20.0 deg-C; h2= “hi” kJ/kg; Vel2= 0.0 m/s; z2= 0.0 m; } 

} 


# p; n d of TEST-code 
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# Property spreadsheet: 


# State 

p(kPa) 

T(K) 

X 

v(m3/kg) 

u(kj/kg) h(kj/kg) 

s(kj/kg) 

# 01 

800.0 

304.4 

0.0 

8.0E-4 

94.0 94.68 

0.351 

#02 

133.7 

253.2 

0.3 

0.0487 

88.16 94.68 

0.378 


Prob.5.32. A hot water stream at 80 C enters a mixing chamber with a mass flow rate of 0.5 kg/s where 
it is mixed with a stream of cold water at 20 C. If it is desired that the mixture leave the chamber at 42 
C, determine the mass flow rate of the cold water stream. Assume that all the streams are at a pressure 
of 250 kPa. [Ref. 1] 


TEST Solution: 

1. This is a problem on mixing chambers. So, choose the appropriate daemon as shown below: 




Basic Tools 


?skCt 


System Analysis 
, I . 


1 

States & Properties 


Closed 


Open 


I 

Uniform 

System 


1 

Uniform 

Flow 


Unsteady Process 


Steady State Steady State 
( cycles ) I 


l 

Unsteady Process 


Generic , 

r 


Specific 


Generic - 


I I 

Reciprocating Cycles HVAC/Psyclvometry Comb. & Equilibrium 


1 1 1 1 

Uniform Non-Mixing Semi-Mixing Mixing 

System Non-Uniform Non-Uniform Non-Uniform 


I 

Single-Flow « 


1 

Non-Mixing Multi-Flow 


Specific 


Mixing Multi-Flow 


I 1 1 1 

Vapor and Gas Refrigeration UVAC Gas Combustion 

Power Cycles Cycles Psychrometry Dynamics & Equilibrium 
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2. Hovering the mouse pointer on Mixing Multi-Flow brings up the following: 


Click to go to page: TEST :> Daemons^ Systems^ Open^ Steady:* Generic^ Multi-Flow Mixing Systems 


Multi-Flow Mixing Systems: Analyze a 
mixing open steady system with two 
inlets and a single exit. Examples include 
a mixing chamber where two non- 
reacting gases are mixed or two different 
phases of a fluid are mixed at steady 
state. 

Mixing chambers are covered in chapters 
4, 6, and 11. 


CP 



3. Choose Phase Change (PC) model, and choose H20 as working substance: 


Two Identical Fluids 



The system has two separate flows 
consisting of two identical fluids, say, 
H20 and H20, which can not mix. Two 
inlet and two exit states - il, i2, el 
and e2-states - are necessary to 




>- Launches the PC Model 
PC Mod Generic, Multi-How, Non-Mixing, upen-bteady LJaemon 


describe the non-mixing device. The 
two fluids can be chemically identical. 





SL Model 


IG Model 


fn— sRT 
RG Model 


Examples; Cold liquid water and 
steam are mixed in a steady-flow 
mixing chamber to produce hot water. 
A steady flow of refrigerant is 
separated into flows of saturated 
vapor and saturated liquid in a 
separation chamber. For specific 
examples, click on the help icon at the 
bottom margin of the daemon. 


This e-book 
is made with 

SetaPDF 



SETASIGN 






£1 


lopers 


www.setasign.com 
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4. Enter data for State 1, i.e. PI, Tl, mdotl; click on Calculate (or, press Enter). We get: 

Generic, Open Steady, Multi-Flow, Mixing Daemon: PC Model 



5. Enter data for State 2 (i.e. cold stream entering), i.e, P2 and T2, press Enter: 



6. Now, enter data for State 3 (i.e. state after mixing), i. P3, T3, mdot3 (= mdot2 + mdotl), 
press Enter: 
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7. Go to Device Panel, enter State 1, State 2 and State 3 for il-state, i2-state and el-state 
respectively. e2-state is maintained as Null-state since there is only one exit. 

Press Enter, and also SuperCalculate: 



8. Now, go to State 2: 



Thus: mdot2 = 0.864 kg/s ... Ans. 

9. Go to I/O panel to see TEST code etc: 

# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>MultiFIowMixed>PC-Modei; v- 10. 
bb06 

# start of TEST-code 

States { 

State- 1: H20; 

Given: { pl= 250.0 kPa; Tl= 80.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; mdotl= 0.5 kg/s; } 
State-2: H20; 

Given: { p2= “pi” kPa; T2= 20.0 deg-C; Vel2= 0.0 m/s; z2= 0.0 m; } 
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State-3: H20; 

Given: { p3= 250.0 kPa; T3= 42.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; mdot3= “mdotl+mdot2” 
kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1, State-2; e-State = State-3; Mixing: true; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# E n d of TEST-code 


# Property spreadsheet : 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

250.0 

353.2 

0.001 

334.86 

335.12 

1.075 

#02 

250.0 

293.2 

0.001 

83.96 

84.21 

0.297 

#03 

250.0 

315.2 

0.001 

175.92 

176.17 

0.599 
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# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1, State-2; e-State = State-3; Mixing: true; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= 0.023273543 kW/K; Jdot_net= “-2.842171E-14” kW; Sdot_net= 
-0.023273543 kW/K; 

Verify: 

******* CALCULATE VARIABLES: Type in an expression starting with an c =’ sign ( c = mdotl*(h2-hl)’, 
£ = sqrt(4*Al/PI)’, etc.) and press the Enter key)********* 

(mdotl*hl+mdot2*h2) = 240.3336599692044 
mdot3*h3 = 240.33365996920443 

# i.e. Energy balance is verified. 


Prob. 5.33. Steam enters the condenser of a steam power plant at 20 kPa as sat. vapour with a mass flow 
rate of 20000 kg/h. It is to be cooled by water from a nearby river, circulating the water through the 
tubes within the condenser. The river water is not allowed to experience a temp rise above 10 C. If the 
steam is to leave the condenser as sat. liquid at 20 kPa, determine the mass flow rate of cooling water 
required. [Ref. 1] 

TEST Solution: 

This is a Non-mixing multi-flow type problem, i.e. the steam and cooling water do not mix. 
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1. Choose the daemon suitable for Non-mixing, multi-flow problem, as shown below: 




r 

Basic Tools 




"1 


System Analysis 

r I , 


States & Proper ties 


Un if Desk Cal Tab les 
Converter Chans 


Open 


Unsteady Process 


Closed 

1 

Steady Sfrrte Steady State 
(cycles) I 


I 

Uniform 

System 


Uniform 
Row 

1 

Unsteady Process 


Generic 


Generic ^ Specific 

i 1 i i 

Reciprocating Cycles HVAC/Psychrometry Comb. & Equilibrium 


I I I 1 

Uniform Non-Mixing Semi-Mixing Mixing 

System Non-Uniform Non-Uniform Non-Uniform 


I 

Single-Flow 


Specific 


Non-Mixing Multi-Flow 


Mixing Multi-Flow 


T 


T 


I 

Vapor and Gas Refrigeration HVAC Gas Combustion 

Power Cycles Cycles Psychrometry Dynamics & Equilibrium 


2. Hovering the mouse pointer over “Non-mixing Multi-Flow” gives following window: 


CEick to go to page: TEST> Daemons:* Systems>Open> Steady >Generic> Multi-Flow Non-Mixing Systems 


Multi-Flow Non-Mixing Systems: 

Analyze a non-mixing open steady 
system with two inlets and two exits. A 
co- flow or counter- flow heat exchanger 
is an example of such a system. The 
working substances can be different for 
the two flows. 

Heat exchangers are covered in 
chapters 4 and 6. 



Ftuw 

DirtttivH 



□ Ca Flc-w FSe-* 
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Note that this is the daemon required to solve parallel flow and counter- flow heat exchangers: 

3. Choose PC model under ‘Two Identical Fluids’ as shown below, since water/ steam is the 
working substance: 


Select a material model to launch the non-mixing mufti-flow daemon. 


Two Identical Fluids 


p 



~ j 


C Mnrl 

Launches the PC Model 

Generic, Multi-Flow, Non-Mixing, Open-Steady Daemon 


1 1 


c„ = const 
PG Model 


pv = R1 
IG Model 


pv = zRT 
RG Model 


The system has two separate flows consisting of two identical fluids, say, H20 and H20, which can not mix. Two inlet and 
two exit states - il, i2, el and e2-states - are necessary to describe the non-mixing device. The two fluids can be 
chemically identical. 

Examples: Heat is exchanged in a counter-flow heat exchanger between a flow of cold liquid water and hot steam. For 
specific examples, click on the help icon at the bottom margin of the daemon. 


Two Different Fluids 


IaA 

PC Model 
+ 

PC Model 


The system has two separate flows consisting of two phase-change (PC) fluids, say, H20 and NH3, which can not mix. 
Two inlet and two exit states - il, i2, el and e2-states - are necessary to describe the non-mixing device. The two fluids 
can be chemically identical. 

Examples: R-134a and H20 are the two fluids in a heat exchanger. Suppose both the inlet states, state-1 (il) and state- 
2 (i2), and one of the exit states, state-3 (el), are completely given. For state-4 (e2 state), set mdot4=mdot2, set up 
the device panel with the known value of Wdot_ext(=0) and Qdot (=0, if adiabatic), and click Super-Calculate to 
evaluate State-4. If T3 and T4 are both unknown, but related, iterative solution is necessary in which Qdot is left as an 
unknown, T3 is guessed until Qdot approaches the known value. For specific examples, click on the help icon at the 
bottom margin of the daemon. 



Deloitte 

Discover the truth at www.deloitte.ca/careers © Deloitte & Touche LLP and affiliated entities. 




228 




Click on the ad to read more 


Download free eBooks at bookboon.com 








Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


4. After choosing H20 as the working substance, enter data for State 1, i.e. PI, xl 
(= 1, since sat. vap. is entering the condenser), and hit Enter: 


Generic, Open Steady, Multi-Flow, Non-Mixing Daemon: PC Model 




5. Enter data for State 2 (i.e. sat. liq. leaving the condenser); i.e. enter P2, x2 (=0.0), mdot2 = 
mdotl. Hit Enter: 



6. For State 3, enter data for river water entering the condenser; hit Enter: 
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7. State 4 is river water exiting the condenser; enter the data, i.e. P4, T4, mdot4 = mdot3, 
and hit Enter: 



Note that exit temp of cooling (river) water is 10 C above the inlet temp. 


8. Now, go to Device Panel, enter i-1 state, i-2 state, e-1 state and e-2 state as shown. Also, 
Qdot = 0, and Wdot_ext = 0. Press Calculate, and SuperCalculate: 
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9. Go to State Panel. 


See State 3: 



We see that: mdot3 = 313.49 kg/s ...flow rate of cooling (river) water required... Ans. 
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10. To see the TEST code etc go to I/O panel: 

# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>SteadyState>Generic>MultiFlowUnmixed>PC-Model; v-10. 
bb06 

# start of TEST-code 

States { 

State- 1: H20; 

Given: { pl= 20.0 kPa; xl= 1.0 fraction; Vell= 0.0 m/s; zl= 0.0 m; mdotl= “20000/3600” kg/s; } 
State-2: H20; 

Given: { p2= 20.0 kPa; x2= 0.0 fraction; Vel2= 0.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 
State-3: H20; 

Given: { p3= 100.0 kPa; T3= 25.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; } 

State-4: H20; 

Given: { p4= £C p3” kPa; T4= “T3+10” deg-C; Vel4= 0.0 m/s; z4= 0.0 m; mdot4= “mdot3” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1, State-3; e-State = State-2, State-4; Mixing: false; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 

# ****** DETAILED OUTPUT: 

# Evaluated States: 

# State- 1: H20 > Saturated Mixture; 

# Given: pl= 20.0 kPa; xl= 1.0 fraction; Vell= 0.0 m/s; 

# zl= 0.0 m; mdotl = “20000/3600” kg/s; 

# Calculated: Tl= 60.062 deg-C; yl= 1.0 fraction; vl= 7.6516 m A 3/kg; 

# ul= 2456.7214 kj/kg; hl= 2609.708 kj/kg; sl= 7.9086 kJ/kg.K; 

# el= 2456.7214 kj/kg; jl= 2609.708 kj/kg; Voldotl= 42.509 m A 3/s; 

# Al= 4250896.5 m A 2; MM1= 18.0 kg/kmol; 

# State-2: H20 > Saturated Mixture; 

# Given: p2= 20.0 kPa; x2= 0.0 fraction; Vel2= 0.0 m/s; 

# z2= 0.0 m; mdot2= “mdotl” kg/s; 

# Calculated: T2= 60.062 deg-C; y2= 0.0 fraction; v2= 0.001 m A 3/kg; 

# u2= 251.3691 kj/kg; h2= 251.3895 kj/kg; s2= 0.832 kJ/kg.K; 
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# e2= 251.3691 kj/kg; j2= 251.3895 kj/kg; Voldot2= 0.0056 m A 3/s; 

# A2= 565.0207 m A 2; MM2= 18.0 kg/kmol; 

# State-3: H20 > Subcooled Liquid; 

# Given: p3= 100.0 kPa; T3= 25.0 deg-C; Vel3= 0.0 m/s; 

# z3= 0.0 m; 

# Calculated: v3= 0.001 m A 3/kg; u3= 104.8785 kj/kg; h3= 104.9788 kj/kg; 

# s3= 0.3673 kJ/kg.K; e3= 104.8785 kj/kg; j3= 104.9788 kj/kg; 

# mdot3= 313.4869 kg/s; Voldot3= 0.3145 m A 3/s; A3= 31447.955 m A 2; 

# State-4: H20 > Subcooled Liquid; 

# Given: p4= “p3” kPa; T4= “T3+10” deg-C; Vel4= 0.0 m/s; 

# z4= 0.0 m; mdot4= “mdot3” kg/s; 

# Calculated: v4= 0.001 m A 3/kg; u4= 146.6718 kj/kg; h4= 146.7725 kj/kg; 

# s4= 0.5052 kJ/kg.K; e4= 146.6718 kj/kg; j4= 146.7725 kj/kg; 

# Voldot4= 0.3154 m A 3/s; A4= 31544.615 m A 2; 

# Property spreadsheet starts: 


# Property spreadsheet starts: 


# State 

p(kPa) 

T(K) 

X 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

20.0 

333.2 

1.0 

7.6516 

2456.72 

2609.71 

7.909 

# 02 

20.0 

333.2 

0.0 

0.001 

251.37 

251.39 

0.832 

#03 

100.0 

298.2 


0.001 

104.88 

104.98 

0.367 

# 04 

100.0 

308.2 


0.001 

146.67 

146.77 

0.505 


# Mass, Energy, and Entropy Analysis Results: 

# Device- A: i-State = State- 1, State-3; e-State = State-2, State-4; Mixing: false; 

# Given: Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; 

# Calculated: Sdot_gen= 3.9177132 kW/K; Jdot_net= 0.0 kW; Sdot_net= -3.9177132 kW/K; 


Prob.5.34. Air enters an adiabatic horizontal nozzle at 400 C with a velocity of 50 m/s. The inlet area is 
240 cm A 2. Temp of air at exit is 80 C. Given that the sp. vol. of air at the inlet and exit are respectively 0.2 
m A 3/kg and 1.02 m A 3/kg, find the area of cross-section of the nozzle at the exit. Assume that enthalpy 
of air is a function of temp only and that cp = 1.005 kJ/kg.K. [VTU-BTD-July 2006:] 
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TEST Solution: 

This problem is the same as Prob.5.13 which was solved with EES. 

1. Go to the Daemon tree and locate System Analysis - Open - Generic - Single Flow: 


r 

Basic Tools 


1 

Unit Desk Cal Tables 3 
Co n werfer Charts 


System Analysis 

* I , 


States & Pro perties 


Closed 


Open 


I 

Uniform 

System 


1 

Uniform 

Flow 


I 

Unsteady Process 


Steady Slate Steady State 
(cycles) I 


Unsteady Process 


Generic 

r 


Specific 

I 


Generic 


I 1 

Reciprocating Cycles HVAC/Psychrometry Comb. 3 Equilibrium 


I 

Uniform 

System 


I 

Non-Mixing 

Non-Uniform 


I 

Semi-Mixing 

Non-Uniform 


1 

Mixing 

Non-Uniform 



Single-Flo 


I 

Non-Mixing Multi-Flow 


I 

Mixing Multi-Flow 


Specific 


Vapor and Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 

Psychrometry 


Gas 

Dynamics 


Combustion 
3 Equilibrium 


2. Select the Perfect Gas (PG) Model (cp = const.) for Material model, since air is the 
working substance: 


<ST 

l K. 

PC Model 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

Working fluids such as H20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 
possibility of a phase transition. 


SL Model 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 

Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simple st gas model which obeyes the ideal gas equation (pv = RT) and assumes 
Launches the PG Model eats to be constant. In the IG (ideal gas) model, specific heats are assumed 

PG Mo e| Generic, Single-Flow, Open-Steady Daemon I ction 0 f temperature only. The RG (real gas) model uses generalized 


~pv = zRT 


I compressiblity charts and is useful for qases near the critical or super-critical 
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3. We get the following screen after clicking on PG model. Now, choose Air as the Working 
substance from the drop down menu. Then, enter known values of Tl, Veil, vl and A1 for 
State 1. Click on Calculate. We get: 


Generic, Open Steady, Single-Flow, Daemon: PG Model 


> Daemons > Systems > Open > Steady > Generic > SingleFlow > PG-Model 

^■Era^»l*rSB K, m 

vl =0.2 m A 3/kg {Specific volume] 
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4. Enter data i.e. T2, v2 and mdot2 = mdotl for State 2, hit Enter: 



5. Go to Device Panel, enter State 1 and State 2 for i-state and e-state respectively; 
enter Qdot = 0 and Wdot_ext = 0 for the nozzle and click on Calculate. We get: 


Generic, Open Steady, Single-Flow, Daemon: PG Model 


theimofluids.net > Daemons > Systems > Open > Steady > Generic > SingleFlow > PG-Model 

I mi m ar ra ■ *, m 

Move mouse over a variable to display its value with more precision. 



6. Now, click on SuperCalculate. Go to State Panel. We get: 


State 1: 
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And, State 2: 



Thus: A2 = 76.22 cm A 2 , Vel2 = 802.95 m/s.... Ans. 
Also, p2 = 99.36 kPa, mdotl = mdot2 = 6 kg/s ... Ans. 

7. Indicative T-s diagram for Plots tab: 

T, K 

740.46 


317.83 

6 35 s, kJ/kg.K 7 76 


( 6 . 19 , 726 . 7-46 ) 
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8. I/O panel gives the TEST code etc.: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { Tl= 400.0 deg-C; vl= 0.2 m A 3/kg; Vell= 50.0 m/s; zl= 0.0 m; Al= 240.0 cm A 2; } 
State-2: Air; 

Given: { T2= 80.0 deg-C; v2= 1.02 m A 3/kg; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# p; n d of TEST-code 
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and stimulating program that will help them reach their full potential. 
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# Property spreadsheet starts: # 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

965.92 

673.2 

0.2 

183.13 

376.31 

7.053 

# 

2 

99.36 

353.2 

1.02 

-46.16 

55.19 

7.058 


# 


Prob.5.35. Steam at 1 MPa and 250 C enters a nozzle with a velocity of 60 m/s and leaves at 10 kPa. 
Assuming the flow process to be isentropic and the mass flow rate to be 1 kg/s, determine: (i) the exit 
velocity (ii) the exit diameter. [VTU-BTD- Jan./Feb. 2005] 



TEST Solution: 

1. Choose the System analysis ... Single Flow daemon as shown below: 


Basic fools 


r 


1 


System Analysis 

, I , 


I 

States & Pr operties 


Unit DeskCat Tables & 
Co n veiTer C harts 


Open 


Unsteady Process 


Closed 

1 

Steady Slate Steady State 
(cycles) I 


I 

Uniform 

System 


1 

Uniform 

Flow 


Unsteady Process 


Generic 


Specific 
I 


Generic 


ir i i 

Reciprocating Cycles HVAC/Psychrometry Comb. Equilibrium 


I I I 1 

Uniform Non-Mixing Semi-Mixing Mixing 

System Non-Uniform Non-Uniform Non-Uniform 



[Single-Flo 


I 

Non-Mixing Multi-Flow 


I 

Mixing Multi-Flow 


Specific 


T 


T 


I 

Vapor and Gas Refrigeration HVAC Gas Combustion 

Po wer Cycles Cycles Psy chro me try Dynamics 3 Equilibrium 
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2. Choose the PC model for Material model: 





Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 




compressed-liquid sub-model, except for species with an asterisk (H20 3C as opposed to H20), which uses 
compressed liquid table for better accuracy. 



Model 

: low, Open-Steady Daemon 
Virsammiy ui a pnaae ua 


PC Mode 

Launcnes tne k<_ 
Generic, Single-F 

20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 
nsition. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


Pure Solid and Pure Liquid: Constant density and constant specific heats ( c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 

Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 

Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation (pv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 


3. Select H20 as working substance. Enter data, i.e. PI, Tl, Veil and mdotl for State 1. 
Hit Enter: 


c„ = const 
PG Model 


pv = zRT 
RG Model 


SL Model 


Generic, Open Steady, Single-Flow, Daemon: PC Model 
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4. Enter P2,s2 = si (since isentropic) and mdot2 = mdotl and hit Enter: 



Note in the above screen shot that immediately other parameters for State 2 are calculated. 
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5. Go to Device Panel, enter State 1 and Styate 2 for i-state and e-state; also, enter Qdot = 0 
and Wdot_ext = 0. Hit Enter: 



6. Now, click on SuperCalculate. Then, go to State Panel. 


See State 1: 



And, State 2: 
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Thus: 

exit velocity, Vel2 = 1225.31 m/s, 
exit area = 0.01002 m A 2 = 100.188 cm A 2 
Therefore, exit dia = d2 = 11.29 cm 

7. I/O panel gives TEST code etc: 

# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>PC-Model; v-10.bb06 

# start of TEST-code 

States { 

State- 1: H20; 

Given: { pl= 1000.0 kPa; Tl= 250.0 deg-C; Vell= 60.0 m/s; zl= 0.0 m; mdotl= 1.0 kg/s; } 
State-2: H20; 

Given: { p2= 10.0 kPa; s2= “si” kJ/kg.K; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; Wdot_ext= 0.0 kW; T_B= 298.15 K; } 

} 


# p n d of TEST-code 


# Property spreadsheet starts: 


# State 

p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

1000.0 

523.2 

0.2327 

2709.89 

2942.57 

6.925 

# 02 

10.0 

319.0 0.8 

12.2762 

2070.91 

2193.67 

6.925 


# ****** CALCULATE VARIABLES: 

sqrt(A2 * 4 /pi) = 0.11294417026465142 m 

# i.e. d2 = 11.29 cm ...Ans. 
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Prob. 5.36. Air flows steadily through a rotary compressor. At entry the air is at 20 C and 101 kPa. At 
exit it is at 200 C and 600 kPa. Assuming the flow to be adiabatic, (i) evaluate the work done per unit 
mass of air if the velocities at inlet and exit are negligible, (ii) what would be the increase in work input 
if the velocities at inlet and exit are 50 m/s and 110 m/s? [VTU-BTD-Ja./Feb. 2004] 


Out 



Fig.Prob.5.36 
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TEST Solution: 

1 . Choose the System analysis . . . Single Flow daemon as shown below: 


Basic fools 


System Analysis 

, I , 


States & Properties 


'ts 


Open 


Unsteady Process 


Closed 

1 

Steady State Steady State 
(cycles) I 


I 

Uniform 

System 


1 

Uniform 

Flow 


1 

Unsteady Process 


Generic 


r 


Specific 


Generic 


T 


1 


Reciprocating Cycles HVAC/Psycfro me try Comb. & Equilibrium 


I 

Uniform 

System 


I 

Non-Mixing 

Non-Uniform 


I 

Semi-Mixing 

Non-Uniform 


1 

Mixing 

Non-Uniform 



[Single- Flo 


I 

Non-Mixing Multi-Flow 


I 

Mixing Multi-Flow 


Specific 


Vapor and Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 

Psychrometry 


I 

Gas 

Dynamics 


Combustion 
& Equilibrium 


2. Choose the IG model for Material model: 


c„ - const 
PG Model 


pv=RT 

IG Model 


pv = zRT 
RG Model 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation {pv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 


Launches the IG Model 
Generic, Single-Flow, Open-Steady Daemon 


im expands steadily in a nozzle from an inlet-state to an exit-state with 
phase change. For specific examples, click on the help icon at the 


bottom margin of the daemon. 


A. When velocities are negligible: 
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3. Choose Air as working substance. Enter PI, Tl, mdotl for State 1. Veil = 0 by default. Hit 
Enter. We get: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 



^ Mixed SI ^ English < [ ©Case-0 v j > | Iv Help Messages C 
State Panel Device Panel 


Super-lterate 


Super-Caiculate 


IE3 


Exergy Panel 


jJ ©State-1 v 


No Yes 


I* pi I* Tl I - rhol I - " vl f“ ul 


1101.0 

■ i kPa 

v 20.0 

^jdeg*C 

v |l .20052 

kg/m^ 

V io.83297 

m^g 

v 1*89.13149 

kJ/kg 

- 

r hi 


r si 


W Veil 


w 

zl 


r el 



|-5.00139 

jlU/Vg 

v| (6.86686 

kJ/kg. K 

fc t'o o 

1 m,s 

v 0.0 


■ m 

V 1-89.13149 

U/kg 


r ji 


phil 


pstl 


w 

mdotl 


r Voldotl 



1*5.00139 

kJ/kg 

3r 

KJ/kg 


\kMg 

vj 1.0 


|jkg/s 

v||0.83297 

111*3/8 

v 


r ai 


MM1 


R1 


r c _pi 



(83297.13 

m*2 

V 28.97 

kg/kmol 

v ]0. 28699 

kJ/kg. K 

V 1 1.00347 

fcJ/kg.K 



4. Enter P2, T2, mdot2=mdotl for State 2; press Enter: 



5. Go to Device Panel, enter State 1 and State 2 for i-state and e-state respectively. 
Also, Qdot = 0, Wdot_ext = 0. Click Calculate, and SuperCalculate. We get: 



Thus: Work done on unit mass of air = -183.08 kW.... (negative sign indicates work input since it is 
a compressor). .Ans. 
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6. Indicative T-s diagram: 



6.16 


s, kJ/kg.K 


7.55 
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7. I/O panel gives TEST code etc: 


# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: Air; 

Given: { pl= 101.0 kPa; Tl= 20.0 deg-C; Vell= 0.0 m/s; zl= 0.0 m; mdotl= 1.0 kg/s; } 
State-2: Air; 

Given: { p2= 600.0 kPa; T2= 200.0 deg-C; Vel2= 0.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device-A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; T_B= 298.15 K; } 

} 


# p n( i of TEST-code 


# Property spreadsheet starts: 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

101.0 

293.2 

0.833 

-89.13 

-5.0 

6.867 

# 

2 

600.0 

473.2 

0.2263 

42.29 

178.08 

6.842 


B. When Veil = 50 m/s and Vel2 =110 m/s: 
The procedure is: 

i) Enter the Veil value , Calculate, and 

ii) enter Vel2 value, Calculate, and then 

iii) SuperCalculate. 
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We get: 
State 1: 


Generic, Open Steady, Single-Flow, Daemon: IG Model 


^ > Systems > Open > Steady > Generic > SinqleFlow > IG-Model 

I ee®i ra a? m ■ *, m 

Move mouse over a variable to display its value with more precision. 



State 2: 



Device Panel: 



Thus: Work done on unit mass of air = -187.88 kW...when inlet and exit velocities are 
considered. . . Ans. 
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# OUTPUT OF SUPER-CALCULATE # 

# Daemon Path: Systems>Open>SteadyState>Generic>SingleFlow>IG-Model; v-10.bb05 

# start of TEST-code 

States { 

State- 1: Air; 

Given: { pl= 101.0 kPa; Tl= 20.0 deg-C; Vell= 50.0 m/s; zl= 0.0 m; mdotl= 1.0 kg/s; } 
State-2: Air; 

Given: { p2= 600.0 kPa; T2= 200.0 deg-C; Vel2= 110.0 m/s; z2= 0.0 m; mdot2= “mdotl” kg/s; } 

} 

Analysis { 

Device- A: i-State = State- 1; e-State = State-2; 

Given: { Qdot= 0.0 kW; T_B= 298.15 K; } 

} 


# p n( i of TEST-code 
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# ****** DETAILED OUTPUT: 


# 

State 

p(kPa) 

T(K) 

v(m A 3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 

1 

101.0 

293.2 

0.833 

-89.13 

-5.0 

6.867 

# 

2 

600.0 

473.2 

0.2263 

42.29 

178.08 

6.842 


Prob. 5.37. Steam at a pressure of 1.4 MPa, 300 C is flowing in a pipe. Connected to this pipe through 
a valve is an evacuated tank. The valve is opened and the tank fills with steam until the pressure is 1.4 
MPa, and then the valve is closed. The process takes place adiabatically and K.E. and P.E. are negligible. 
Determine the final temp of steam. [Ref. 2] 


TEST Solution: 

This is a problem on charging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 

1. Select System Analysis - Open - Unsteady Process daemon as shown below: 


thermofIuids.net > Daemon Map 1 

i 

*- Home tff 
W TEST 


-f 



Basic Tools 
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| I | 

Closed Open 

1 l 
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(cycles) 
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1 

Mixing 

Non-Uniform 


I 

Single-Flow 


I 

Non-Mixing Multi- Flow 


I 

Mixing Multi-Flow 


Specific 
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Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 
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I 

Gas 

Dynamics 


Combustion 
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2. Hovering the mouse pointer on ‘Unsteady Process brings up the following message window: 


Click to go to page: TEST> Da emoos>Systems> Open > Unsteady Processes 


Open Processes: The open 
system is unsteady; moreover, 
only the beginning and final 
conditions are relevant so that 
the instantaneous changes can 
be integrated out. However, 
unlike a closed process, the open 
process equations also involve 
the inlet and/or exit state(s). 


Examples (chapters 4, 6) include 
charging and discharging of 
cylinders, inflating a tire, a 
pressure cooker discharging steam, etc. 


l * 

Sm&c ■ 3' 

jf 


0 


2m Icivi 

m k 

x = 1MH 



© 

mk\Ki 

m.iK 


3. Choose PC model for Material model since Steam is the working substance: 


Open-System Unsteady Process Daemons: Select a Material Model 



4. Select H20 as the substance and enter data, i.e. PI = 0, ml = 0 (since tank is evacuated), 
and Voll, for State 1. Press Enter: 


Open Process Daemon: Phase-Change (PC) Model 
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5. For State 2, enter P2, and m2 = m3, not known yet. Press Enter: 



6. Enter data for State 3, i.e. state of fluid in the pipe, P3, T3 and Vol3 = Voll. Press Enter: 



7. Go to Process Panel, enter i-state = State 3, e-state = Null, and State 1 and State 2 for b-state 
and f-state respectively, as shown. Also, Q = 0, W_ext = o. Press Enter: 



Process Panel 



Open Process - A 

Mass: = (m, - m e ) 

Energy [m f e f - ) = ( [mj ; - mj^ ) + Q 

A E -'net Q 

Entropy \(m f s f -m^ b )= (j m & - m e s g )+ — 

' v ' ' y ' A 



WinHip: 

Work in negative 
Heat in positive 


253 


Download free eBooks at bookboon.com 















Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


8. Click on SuperCalculate. Go to State Panel: 


State 2: 



Thus, T2 = 452.1 deg. C .... Ans. 

9. I/O panel gives TEST code etc: 


# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>Process>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: H20; 

Given: { pl= 0.0 kPa; Vell= 0.0 m/s; zl = 0.0 m; ml= 0.0 kg; Voll= 0.4 m A 3; } 

State-2: H20; 

Given: { p2= 1400.0 kPa; u2= a h3” kj/kg; Vel2= 0.0 m/s; z2= 0.0 m; m2= “m3” kg; } 
State-3: H20; 

Given: { p3= 1400.0 kPa; T3= 300.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; Vol3= “voll” m A 3; } 

} 

Analysis { 

Process- A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

Given: { Q= 0.0 kj; W_ext= 0.0 kj; T_B= 298.15 K; } 

} 


# p n( i of TEST-code 
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# Property spreadsheet: 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

0.0 

0.0 0.0 

Infinity 

0.0 

0.0 

0.0 

#02 

1400.0 

725.3 

0.2357 

3040.33 

3370.3 

7.459 

#03 

1400.0 

573.2 

0.1823 

2785.14 

3040.33 

6.953 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

# Given: Q= 0.0 kj; W_ext= 0.0 kj; T_B= 298.15 K; 

# Calculated: S_gen= 1.1094596 kJ/K; Delta_E= 6671.911 kj; Jdot_net= 6671.911 kj; 
Delta_S= 16.368185 kJ/K; 

# Sdot_net= 15.258725 kj; 


Prob.5.38. A 1 m A 3 tank contains ammonia at 150 kPa, 25 C. The tank is attached to a line flowing 
ammonia at 12300 kPa, 60 C. The valve is opened, and mass flows in until the tank is half full of liquid, 
by volume at 25 C. Calculate the heat transferred from the tank during this process. [Ref. 2] 
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TEST Solution: 

This is a problem on charging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 

1. Select System Analysis - Open - Unsteady Process daemon as shown below: 
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2. Choose PC model for Material model since NH3 is the working substance: 


Open-System Unsteady Process Daemons: Select a Material Model 


1 themnofluids.net > Daemons 

> Systems 

> Open > Process 1 

1 4t> Home of 1 

Ivteit 1 

m 

< 3 ? 

rail," hi 


Launch the open process eaemon by selecting a sutiable material model. 


PC Model 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

Launches the PC Model R-12, NH3 R_i 34 a N2 # C02, etc., should be treated as PC fluids if there is any 

Open-System Unsteady Process Daemon 



Examples: Saturated water vapor escapes from a pressure cooker. To determine the heat transfer during 
such an open process. For specific examples, click on the help icon at the bottom margin of the daemon. 
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I Law of Thermodynamics applied to Flow Processes 


3. Select Ammonia (NH3) as the substance and enter data, i.e. PI, Tland Voll for State 1. 
Press Enter: 

Open Process Daemon: Phase-Change (PC) Model 




4. Enter the data, viz. T2, Vol2 = Voll, and y2 = volume fraction = 0.5 for State 2; press Enter: 



5. Enter data for State 3 (i.e. fluid flowing in the line), i.e. P3, T3; press Enter: 
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I Law of Thermodynamics applied to Flow Processes 


6. Go to Process panel, enter State 3 for i-state, Null for e-state, and State 1 and State 2 for 
b-state and f-state respectively. (See the fig. below). Press Enter, and click SuperCalculate: 



Thus: Q = -379647.3 kj .... Ans... (negative sign indicates that heat is rejected). 
7. I/O panel gives TEST code etc: 


# OUTPUT OF SUPER-CALCULATE : 

# Daemon Path: Systems>Open>Process>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: Ammonia(NH3); 

Given: { pl= 150.0 kPa; Tl= 25.0 deg-C; Veil- 0.0 m/s; zl= 0.0 m; Voll= 1.0 m A 3; } 
State-2: Ammonia(NH3); 

Given: { T2= 25.0 deg-C; y2= 0.5 fraction; Vel2= 0.0 m/s; z2= 0.0 m; Vol2= “voll” m A 3; } 
State-3: Ammonia(NH3); 

Given: { p3= 1200.0 kPa; T3= 60.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; } 

} 

Analysis { 

Process- A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

Given: { W_ext= 0.0 kj; T_B= 298.15 K; } 

} 

# p n d of TEST-code 
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# Property spreadsheet starts: 


# State 

p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

150.0 

298.2 

0.9557 

1380.43 

1523.78 

6.113 

# 02 

1003.2 

298.2 0.0 

0.0033 

309.85 

313.14 

1.171 

#03 

1200.0 

333.2 

0.1238 

1404.74 

1553.27 

5.236 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

# Given: W_ext= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= -379647.3 kj; S_gen= 30.783335 kJ/K; Delta_E= 93206.055 kj; Jdot_net= 
472853.38 kj; 

# Delta_S= 351.2871 kJ/K; Sdot_net= 1593.847 kj; 
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I Law of Thermodynamics applied to Flow Processes 


Prob. 5.39. A 0.12 m A 3 rigid tank initially contains refrigerant R134a at 1 MPa and 100% quality The 
tank is connected by a valve to a supply line that carries R134a at 1.2 MPa and 36 C. Now the valve is 
opened and the refrigerant is allowed to enter the tank. The valve is closed when it is observed that the 
tank contains sat. liquid at 1.2 MPa. Determine (a) the mass of R134a that has entered the tank (b) the 
amount of heat transfer [Ref. 1] 


TEST Solution: 

This is a problem on charging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 

1. Select System Analysis - Open - Unsteady Process daemon as shown below: 


thermo-fluids.net > Daemon Map 
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I Law of Thermodynamics applied to Flow Processes 


2. Choose PC model for Material model since R134a is the working substance: 


Open-System Unsteady Process Daemons: Select a Material Model 


thermofluids.net > Daemons 

> Systems 

> Open > Process 1 


w 

\ZJMa y\ 


Launch the open process eaemon by selecting a sutiable material model. 


^ o'- 

■\K., 

PC Model 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

oZ C -s" tern U C n*e°ad!, Process Daemon R - 12 ' NH3 ' R - 134a ' N2 ' C02 ' etc - should be ***** as PC fluids d *** is ar "' 



Examples: Saturated water vapor escapes from a pressure cooker. To determine the heat transfer during 
such an open process. For specific examples, click on the help icon at the bottom margin of the daemon. 


3. Select R134a as the substance and enter data, i.e. PI, xland Voll for State 1. Press Enter: 


Open Process Daemon: Phase-Change (PC) Model 



Move mouse over a variable to display its value with more precision. 
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Observe that mass ml is immediately calculated. 


4. Enter P2, x2 and Vol2 = Voll for State 2. Press Enter. Immediately, mass m2 is calculated: 
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I Law of Thermodynamics applied to Flow Processes 


5. Enter data i.e. P3, T3 and m3 = (m2-ml) for the State 3, i.e. the R134a flowing in the line. 
Press Enter: 
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6. Go to Process panel, enter i-state = State 3, e-state = Null; and enter State 1 and State 2 for 
b- state and f- state respectively. Also, W_ext = 0, and press Calculate and SuperCalculate. 
We get: 



Thus: Q = 1031.51 kj (going in to the system), m3 = mi = (m2 - ml) = 128.25 kg .... Ans. 

7. Get TEST code etc from the I/O panel: 

# OUTPUT OF SUPER-CALCULATE: 


# Daemon Path: Systems>Open>Process>PC-Model; v-10.bb06 

# start of TEST-code 


States { 

State- 1: R-134a; 

Given: { pl= 1000.0 kPa; xl= 1.0 fraction; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.12 m A 3; } 
State-2: R-134a; 

Given: { p2= 1200.0 kPa; x2= 0.0 fraction; Vel2= 0.0 m/s; z2= 0.0 m; Vol2= “voll” m A 3; } 
State- 3: R-134a; 

Given: { p3= 1200.0 kPa; T3= 36.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; m3= “m2-ml” kg; } 

} 

Analysis { 

Process-A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

Given: { W_ext= 0.0 kj; T_B= 298.15 K; } 

} 

# p n d of TEST-code 
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# Property spreadsheet starts: 


# State 

p(kPa) 

T(K) 

X 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

1000.0 

312.5 

1.0 

0.0204 

250.17 

270.55 

0.913 

# 02 

1200.0 

319.4 

0.0 

9.0E-4 

116.05 

117.12 

0.421 

#03 

1200.0 

309.2 


9.0E-4 

100.82 

101.85 

0.373 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: ie-State = State-3, State-Null; bf-State = State- 1, State-2; 

# Given: W_ext= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 1031.5059 kj; S_gen= -0.12388586 kJ/K; Delta_E= 14094.179 kj; Jdot 
net= 13062.673 kj; 

# Delta_S= 51.161983 kJ/K; Sdot_net= 47.826183 kj; 


Prob.5.40. A 100-L rigid tank contains carbon dioxide gas at 1 MPa, 300 K. A valve is cracked open, and 
carbon dioxide escapes slowly until the tank pressure has dropped to 500 kPa. At this point the valve is 
closed. The gas remaining inside the tank may be assumed to have undergone a polytropic expansion, 
with polytropic exponent n =1.15. Find the final mass inside and the heat transferred to the tank during 
the process. [Ref:2] 

Note that this problem is the same as Prob.5.19, which was solved earlier with EES. 

Now, we shall solve it with TEST: 

TEST Solution: 

This is a problem on discharging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 
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1. Select System Analysis - Open - Unsteady Process daemon as shown below: 


thermo-fluids.net > Daemon Map 
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I Law of Thermodynamics applied to Flow Processes 


2. Choose PG model for Material model since C02 is the working substance. Select C02 as 
the working substance and enter data, i.e. PI, T1 and Voll for State 1. Press Enter: 


Open Process Daemon: Perfect Gas (PG) Model 

■ 4k Home of 
I & T E S T 


thermofluids.net > Daemons 


Move mouse over a variable to display its value with more precision. 


<* Mixed O si English 
State Panel 


iS 


^ c : * | 

I* Help Messages On Super-lterate [ 


Super-Calcuiate 


B g II 

F? pi 


©State-1 vj > 


□ 


■ " No * Yes 


r hi 


1000.0 


3 |o. 05667 
I* zl 


^ 1*8996.771 kJ/kg vj |-8 940.098 ki/tcfl 

r ei r n 


rrfS/kg 


r si 


I* veil 


[443043 kJ/kg.K V 0 


] 1*8996.771 [ kJ/kg |-8940.098 [ fcJ/l eg 

W Voll MM1 


I - ml 


' |o.1 V ; |44.01 kg/kmol 

kl 


kJ/kg 


R1 


C _p1 


[1.76449 
c vl 


|0 18891 kJ/kg.K 


kJ/kg.K v |o. 65476 kJ/kg.K 


I UnitLess v 


Note that ml is calculated as 1.76449 kg. 

3. Enter P2, T2, Vol2 = Voll for State 2; press Enter: 



Note that m2 is calculated as 0.96573 kg. 
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4. State 3 is the state of gas flowing out. Its enthalpy goes on changing during the ‘flowing out’ 
process. Let us take the enthalpy as the average of that at the beginning and end of flow, i.e. 
h3 is average of hi and h2. And m3 is equal to (ml - m2). Press Enter: 



5. Go to Process panel, enter i-state = Null, e-state = State 3; and enter State 1 and State 2 for 
b-state and f-state respectively. Also, W_ext = 0. Click on Calculate and SuperCalculate. 
We get: 


Move mouse over a variable to display its value with more precision. 
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Thus: Q = 20.13 kj (heat transferred in to the tank), m2 = 0.966 kg, T2 = 274.1 K .... Ans. 
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6. TEST code and other details can be seen in the I/O panel: 

# OUTPUT OF SUPER-CALCULATE: 

# Daemon Path: Systems>Open>Process>PG-Model; v-10.bb05 

# start of TEST-code 


States { 

State- 1: C02; 

Given: { pl= 1000.0 kPa; Tl= 300.0 K; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.1 m A 3; } 

State-2: C02; 

Given: { p2= 500.0 kPa; T2= “T1*(P2/P1) A ((1.15 - 1)/1.15)” K; Vel2= 0.0 m/s; z2= 0.0 m; 
Vol2= “Voir m A 3; } 

State-3: C02; 

Given: { h3= cc (hl+h2)/2” kj/kg; Vel3= 0.0 m/s; z3= 0.0 m; } 

} 


Analysis { 

Process- A: ie-State = State-Null, State-3; bf-State = State- 1, State-2; 
Given: { W_ext= 0.0 kj; T_B= 298.15 K; } 

} 

# p n d of TEST-code 
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# ****** DETAILED OUTPUT: 

# Evaluated States: 

# State- 1: C02 > PG-Model; 

# Given: pl= 1000.0 kPa; Tl= 300.0 K; Vell= 0.0 m/s; 

# zl= 0.0 m; Voll= 0.1 m A 3; 

# Calculated: vl= 0.0567 m A 3/kg; ul= -8996.771 kj/kg; hl= -8940.098 kj/kg; 

# sl= 4.4304 kJ/kg.K; el= -8996.771 kj/kg; jl = -8940.098 kj/kg; 

# ml= 1.7645 kg; MM1= 44.01 kg/kmol; Rl= 0.1889 kJ/kg.K; 

# c_pl= 0.8437 kJ/kg.K; c_vl= 0.6548 kJ/kg.K; kl= 1.2885 UnitLess; 

# State-2: C02 > PG-Model; 

# Given: p2= 500.0 kPa; T2= “T1*(P2/P1) A ((1.15 - 1)/1.15)” K; Vel2= 0.0 m/s; 

# z2= 0.0 m; Vol2= “Voll” m A 3; 

# Calculated: v2= 0.1036 m A 3/kg; u2= -9013.751 kj/kg; h2= -8961.977 kj/kg; 

# s2= 4.4851 kJ/kg.K; e2= -9013.751 kj/kg; j2= -8961.977 kj/kg; 

# m2= 0.9657 kg; MM2= 44.01 kg/kmol; R2= 0.1889 kJ/kg.K; 

# c_p2= 0.8437 kJ/kg.K; c_v2= 0.6548 kJ/kg.K; k2= 1.2885 UnitLess; 

# 

# State-3: C02 > PG-Model; 

# Given: h3= “(hl+h2)/2” kj/kg; Vel3= 0.0 m/s; z3= 0.0 m; 

# Calculated: T3= 287.0334 K; u3= -9005.262 kj/kg; e3= -9005.262 kj/kg; 

# j3= -8951.037 kj/kg; m3= 0.7988 kg; MM3= 44.01 kg/kmol; 

# R3= 0.1889 kJ/kg.K; c_p3= 0.8437 kJ/kg.K; c_v3= 0.6548 kJ/kg.K; 

# k3= 1.2885 UnitLess; 

# Property spreadsheet starts: 

# State p(kPa) T(K) v(m A 3/kg) u(kj/kg) h(kj/kg) s(kj/kg) 

# 1 1000.0 300.0 0.0567 -8996.77 -8940.1 4.43 

# 2 500.0 274.1 0.1035 -9013.75 -8961.98 4.485 

# 3 287.0 -9005.26 -8951.04 

# Property spreadsheet ends 

# Mass, Energy, and Entropy Analysis Results: 

# Process- A: ie-State = State-Null, State-3; bf-State = State- 1, State-2; 

# Given: W_ext= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 20.132648 kj; Delta_E= 7169.9126 kj; Jdot_net= 7149.78 kj; Delta_S= 
-3.4860783 kJ/K; 
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Prob.5.41. A rigid tank has a volume of 0.06 m A 3 and initially contains two phase liquid-vapour 
mixture of H20 at a pressure of 15 bar and a quality of 20%. As the tank contents are heated, a pressure 
regulating valve keeps the pressure constant in the tank by allowing sat. vap. to escape. Neglecting KE 
and PE changes (a) determine the total mass in the tank, in kg and the amount of heat transfer, in kj, 
if heating continues until the final quality is 0.5 (b) plot the total mass in the tank, and the amount of 
heat transfer versus final quality, x, ranging from x = 0.2 to 1. [Ref. 3] 

TEST Solution: 

This is a problem on discharging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 


1. Select System Analysis - Open - Unsteady Process daemon as shown below: 
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2. Choose PC model for Material model since H20 is the working substance. 





Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 

pv 



compressed-liquid sub-model, except for species with an asterisk (H20* as opposed to H20), which uses 
compressed liquid table for better accuracy. 

PC Mode] 

Launches the PC 
Generic, Single-F 

Model 

r !ow, Open-Steady Daemon 
pj33iL>mcyw-a pnaae ua 

20, R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 
nsition. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 

Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 

Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation (pv = P7) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 


c p - const 
PG Model 


pvy&T 
RG Model 



SL Model 
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3. Select H20 as the working substance and enter data, i.e. PI, xl and Voll for State 1. 
Press Enter: 


Open Process Daemon: Phase-Change (PC) Model 






UrKrilvl 

<3* 


P? tB 

I Move mouse over a variable to display its value with more precision. j 



4. Enter P2 = PI, x2 and Vol2=Voll for State 2, press Enter: 



5. For the fluid flowing out, it is State 3. Enter P3 = PI, x3 = 1 (since, by data, it is sat. vap.) 
and m3 = (ml -m2). Press Enter: 
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6. Go to Process panel. Enter i-state = Null, e-state = State 3 (i.e. fluid flowing out). Also, enter 
States 1 and 2 for b-state and f-state respectively. Enter W_ext = 0. Press Enter. We get: 



Note that Q is calculated as: Q = 2545.65 kj ...= heat supplied.... Ans. 

7. Click on SuperCalculate: Go to I/O panel to see TEST code etc: 

# OUTPUT OF SUPER-CALCULATE : 


# Daemon Path: Systems>Open>Process>PC-Model; v-lO.cbOl 

# start of TEST-code 


States { 

State- 1: H20; 

Given: { pl= 15.0 bar; xl= 0.2 fraction; Vell= 0.0 m/s; zl= 0.0 m; Voll= 0.06 m A 3; } 
State-2: H20; 

Given: { p2= “pi” bar; x2= 0.5 fraction; Vel2= 0.0 m/s; z2= 0.0 m; Vol2= “voll” m A 3; } 
State-3: H20; 

Given: { p3= ££ pl” bar; x3= 1.0 fraction; Vel3= 0.0 m/s; z3= 0.0 m; m3= “ml-m2” kg; } 

} 

Analysis { 

Process- A: ie-State = State-Null, State-3; bf-State = State- 1, State-2; 

Given: { W_ext= 0.0 kj; T_B= 298.15 K; } 

} 

# p n d of TEST-code 


273 


Download free eBooks at bookboon.com 










Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


# ****** DETAILED OUTPUT: 


# Evaluated States: 


# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 

# 


State- 1: H20 > Saturated Mixture; 

Given: pl= 15.0 bar; xl= 0.2 fraction; Vell= 0.0 m/s; 
zl= 0.0 m; Voll= 0.06 m A 3; 

Calculated: Tl= 198.3066 deg-C; yl= 0.9662 fraction; vl= 0.0273 m A 3/kg; 
ul= 1193.3645 kj/kg; hl= 1234.291 kj/kg; sl= 3.1409 kJ/kg.K; 
el= 1193.3645 kj/kg; jl= 1234.291 kj/kg; ml= 2.1979 kg; 

MM1= 18.0 kg/kmol; 

State-2: H20 > Saturated Mixture; 

Given: p2= “pi” bar; x2= 0.5 fraction; Vel2= 0.0 m/s; 
z2= 0.0 m; Vol2= “voll” m A 3; 

Calculated: T2= 198.3066 deg-C; y2= 0.9913 fraction; v2= 0.0665 m A 3/kg; 
u2= 1718.751 kj/kg; h2= 1818.4697 kj/kg; s2= 4.3799 kJ/kg.K; 
e2= 1718.751 kj/kg; j2= 1818.4697 kj/kg; m2= 0.902 kg; 

MM2= 18.0 kg/kmol; 
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# State-3: H20 > Saturated Mixture; 

# Given: p3= “pi” bar; x3= 1.0 fraction; Vel3= 0.0 m/s; 

# z3= 0.0 m; m3= “ml-m2” kg; 

# Calculated: T3= 198.3066 deg-C; y3= 1.0 fraction; v3= 0.1319 mA3/kg; 

# u3= 2594.395 kj/kg; h3= 2792.101 kj/kg; s3= 6.445 kJ/kg.K; 

# e3= 2594.395 kj/kg; j3= 2792.101 kj/kg; Vol3= 0.1709 m^3; 

# MM3= 18.0 kg/kmol; 

# Property spreadsheet starts: 


# Property spreadsheet starts: 


# State 

p(kPa) 

T(K) 

X 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

1500.0 

471.5 

0.2 

0.0273 

1193.36 

1234.29 

3.141 

# 02 

1500.0 

471.5 

0.5 

0.0665 

1718.75 

1818.47 

4.38 

# 03 

1500.0 

471.5 

1.0 

0.1319 

2594.4 

2792.1 

6.445 


# Mass, Energy, and Entropy Analysis Results: 

# Process- A: ie-State = State-Null, State-3; bf-State = State- 1, State-2; 

# Given: W_ext= 0.0 kj; T_B= 298.15 K; 

# Calculated: Q= 2545.6465 kj; S_gen= -3.1388645 kJ/K; Delta_E= -1072.5034 kj; Jdot 
net= -3618.15 kj; 

# Delta_S= -2.9524584 kJ/K; Sdot_net= -8.351734 kj; 


(b) Plot m2, Q against final quality, x2: 

The procedure is simple: Go to State 2, change x2 to desired value, press Calculate, and then SuperCalculate. 
All calculations are immediately up-dated. Read the values of m2 from State 2, and Q from the Process 
panel. Do this for all desired values of x2 and tabulate as shown below: 


x2 

m2 (kg) 

Q(kJ) 

0.2 

2.19789 

0 

0.3 

1.486 

1398.07 

0.4 

1.123 

2112.21 

0.5 

0.902 

2545.65 

0.6 

0.754 

2836.7 

0.7 

0.648 

3045.63 

0.8 

0.567 

3202.9 

0.9 

0.505 

3325.56 

1 

0.455 

3423.9 
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Now, plot the results in EXCEL: 
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Prob.5.42. A well insulated rigid tank of volume 10 m A 3 is connected to a large steam line through 
which steam flows at 15 bar and 280 C. The tank is initially evacuated. Steam is allowed to flow into the 
tank until the pressure inside is P. (a) Determine the amount of mass in the tank, and the temp in the 
tank , when P = 15 bar (b) Plot the quantities in part (a) versus P ranging from 0.1 to 15 bar. [Ref. 3] 

TEST Solution: 

This is a problem on charging a tank. i.e. Uniform State, Uniform Flow type of problem. See eqn. 
5.14 at the beginning of this chapter. 
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1. Select System Analysis - Open - Unsteady Process daemon as shown below: 


thermofIyids.net > Daemon Map 1 

i 

Home of 

^ TEST 


m 



I 

Basic Tools 


System Analysis 
, I , 


I 

States 8t Pr operties 


Closed 


Open 


I 

Uniform 

System 


1 

Uniform 

Flow 


I 

Unsteady Process 


Steady Stale Steady State 
(cycles) I 


Unsteady Process I 


Generic , 


r 


Specific 

I 


Generic - 


T 


Reciprocating Cycles HVAC/Psychrometry 


I 

Uniform 

System 


I 

Non-Mixing 

Non-Uniform 


T 


1 

Comb . 3 Equilibrium 
1 


Semi-Mixing 

Non-Uniform 


Mixing 

Non-Uniform 


1 

Single-Flow $ 


I 

Non- Mixing Multi-Flow 


I 

Mixing Multi-Flow 


I 

Specific 


1 

Vapor am) Gas 
Power Cycles 


I 

Refrigeration 

Cycles 


I 

HVAC 

Psychrometry 


I 

Gas 

Dynamics 


Combustion 
3 Equilibrium 


2. Choose PC model for Material model since H20 is the working substance. 


\K. 


Pure Phase-Transition Fluid: The phase-change (PC) model can be used to determine states of sub-cooled 
(compressed) liquid, super-heated vapor, and saturated mixture of liquid and vapor phases. Based on the 
saturation and super-heated tables, the model is quite accurate. Sub-cooled liquid is modeled with the 
compressed-liquid sub-model, except for species with an asterisk (H20 se as opposed to H20), which uses 
compressed liquid table for better accuracy. 


PC Mode ^ aunches the p C Model 20 , R-12, NH3, R-134a, N2, C02, etc., should be treated as PC fluids if there is any 

mout Generic, Single-Flow, Open-Steady Daemon 

puasiLMiiLy ui a pi lose li e/nsition. 


Examples: Analyze a steady steam turbine with one inlet and one exit. For specific examples, click on the help 
icon at the bottom margin of the daemon. 



SL Model 


Pure Solid and Pure Liquid: Constant density and constant specific heats (c p = c v = c) characterize the 
solid/liquid (SL) model. Beside a wide selection to choose from, a new solid or liquid can be created by 
assigning custom material properties. 

Working substances such as steel, iron, copper, aluminum, wood, water, oil, etc., which can be assumed to 
maintain their condensed (solid or liquid) phase when a system undergoes other changes, can be analyzed 
with the SL model. 


Examples: Liquid water is pumped steadily from a given inlet-state to a given exit-state with no possibility of 
a phase change. For specific examples, click on the help icon at the bottom margin of the daemon. 



Cf, '= const 


PG Model 



pv = zRT 
RG Model 


Pure Gas: A pure gas has a fixed chemical composition across space and time. 
Oxygen, nitrogen, and air are examples of a pure gas. The PG (perfect gas) model is 
the simplest gas model which obeyes the ideal gas equation {pv = RT) and assumes 
specific heats to be constant. In the IG (ideal gas) model, specific heats are assumed 
to be function of temperature only. The RG (real gas) model uses generalized 
compressiblity charts and is useful for gases near the critical or super-critical 
conditions for which PC-model data are not avaiable. 
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3. Select H20 as the working substance and enter data, i.e. Pl= 0, ml= 0 (since evacuated 
tank), and Voll = 10 m A 3 for State 1. Press Enter: 


Open Process Daemon: Phase-Change (PC) Model 


1 1, * rmofluidishet > Daemons > Systems > Open > Process > PC-Model 

m ttvrm m 

Move mouse over a variable to display its value with more precision. 



4. For State 2, enter P2, Vol2, and u2 = h3 for filling an evacuated tank (see eqn. 5.14 at the 
beginning of this chapter). Press Enter: 



5. State 3 refers to the fluid in the line. Enter P3, T3 and m3 = m2. Press Enter: 
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6. Now, go to Process Panel. Enter i-state = State 3, e-state = Null, and enter States 1 and 2 for 
b-state and f-state respectively. Also, enter Q = o (since the tank is insulated) and W_ext = o 
(since no external work). Press Calculate: 



7. Now, click on SuperCalculate. Go to States panel. We get: 


State 1: 



State 2: 
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We see that: T2 = 423.99 C, m2 = 46.46 kg ... Ans. 

8. I/O panel gives TEST code etc: 

# OUTPUT OF SUPER-CALCULATE: 

Daemon Path: Systems>Open>Process>PC-Model; v-lO.cbOl 

# start of TEST-code 

States { 

State- 1: H20; 

Given: { pl= 0.0 bar; Vell= 0.0 m/s; zl= 0.0 m; ml= 0.0 kg; Voll= 10.0 m A 3; } 

State-2: H20; 

Given: { p2= 15.0 bar; u2= “h3” kj/kg; Vel2= 0.0 m/s; z2= 0.0 m; Vol2= “voll” m A 3; } 
State-3: H20; 

Given: { p3= 15.0 bar; T3= 280.0 deg-C; Vel3= 0.0 m/s; z3= 0.0 m; m3= “m2” kg; } 

} 

# p n d of TEST-code 

# Property spreadsheet starts: 


# State p(kPa) 

T(K) x 

v(m3/kg) 

u(kj/kg) 

h(kj/kg) 

s(kj/kg) 

# 01 

0.0 

0.0 0.0 

Infinity 

0.0 

0.0 

0.0 

# 02 

1500.0 

697.1 

0.2107 

2991.81 

3307.88 

7.341 

#03 

1500.0 

553.2 

0.1626 

2747.96 

2991.81 

6.834 


281 


Download free eBooks at bookboon.com 


Basic Thermodynamics: Software Solutions 
Part II 


I Law of Thermodynamics applied to Flow Processes 


(b) Plot m2, T2 against final pressure P2: 

Procedure: Go o State 2, change P2 to desired value, and click on Calculate, and SuperCalculate. 
Immediately, all calculations are updated. Read the values of T2 and m2. Repeat this procedure for all 
desired values of P2. Tabulate the results as shown below: 


P2(bar) 

m2 (kg) 

T2 (deg.C) 

0.1 

0.315 

414.05 

0.5 

1.577 

414.31 

1 

3.154 

414.65 

2 

6.31 

415.32 

4 

12.63 

416.67 

6 

18.95 

418.01 

8 

25.27 

419.35 

10 

31.60 

420.68 

12 

37.94 

422.0 

13 

41.11 

422.67 

15 

47.46 

423.99 
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Now, plot the results in EXCEL: 
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